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Introductory Remarks ' 


A conference on ‘Rheo-optics of Polymers” was held on Saturday, 
August 24 at the University of Massachusetts. This conference, jointly 
sponsored by the Division of High Polymer Physics of the American Physi- 
cal Society and the Polymer Research Institute.of the University of Massa- 
chusetts, was attended by over 100 representatives of universities and 
industry from this country and abroad. It was jointly organized by 
Professor Rodney D. Andrews of the Massachusetts Institute of Technology, 
Dr. Donald G. LeGrand of the General Electric Research Laboratories, 
and Professor Richard 8. Stein, Director of the Polymer Research Institute 
of the University. 

“Rheo-optics” is a new term designating the use of optical methods to 
study flow. This conference was concerned with the use of these methods 
for the study of the flow and deformation of polymers. Methods used 
included birefringence, x-ray diffraction, dichroism, and light scattering. 

The association of these optical methods with standard rheological 
methods aids in the elucidation of mechanisms responsible for their be- 
havior. The study of the time-dependence of the optical behavior often 
permits a determination of the rates of the contributing mechanisms. 

Samples studied in the papers presented at this symposium included 
polymers in the form of solutions, melts, and amorphous and crystalline 
solids. These were studied under conditions of flow, rotational and 
vibrational strain, steady state deformation, relaxation, and creep. 
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Streaming Birefringence of Polymer Solutions 


W. PHILIPPOFF, Products Research Division, Esso Research and 
Engineering Co., Linden, New Jersey 


Synopsis 


Streaming or flow-birefringence is the optical anisotropy induced by laminar flow in 
liquids. The optical anisotropy can be described by an optical tensor which is coaxial 
with the stress tensor in an elastic material In a tension experiment this is obvious. 
In laminar flow this has been proved by investigating the distribution of the bire- 
fringence in three directions of space. The amount of anisotropy, usually described by 
the degree of birefringence An, is proportional to the difference in principal stresses of the 
mechanical stress tensor in a very large range of variables, therefore giving the possi- 
bility of using optical measurements for the determination of mechanical stresses. 
However, the relationship is valid either for concentrated solutions or for conditions 
where the solvent and the polymer have essentially the same index of refraction n—a 
condition usually termed matching solvent. The proportionality constant between the 
difference in principal stresses and An—the stress-optical coefficient C—is, as also follows 
from theory, independent of concentration, shear stress temperature, molecular weight 
and its distribution, only determined by the structure of the monomer unit and its n. 
Another quantity in flow birefringence is the so-called extinction angle which is measured 
between the direction of the principal optical axis and the direction of flow. A reascnable 
amount of experimental data has been obtained to show that this angle can be quan- 
titatively correlated with the so-called ‘‘recoverable shear’’ that is the shear component 
of the elastic strain tensor in laminar flow. Therefore, the experimental determination 
of the angle gives a measure of the elastic strain. A survey of the experimental results 
in the field shows the generalities of this approach in a variety of polymer and solvent 


systems. 


The birefringence induced by mechanical stresses in transparent solids 
has been used for more than half a century as a tool for determining stress 
distributions in models of structures under steady loading.'~* The basis 
of such studies lies in the observance that the amount of birefringence An 
is exactly proportional to the tensile stress in the sample or more exactly 
An is proportional to the difference in principal stresses, and is oriented in 
the direction of the principal tensile stress. This relationship is called the 
stress-optical law, the proportionality constant is called the stress-optical 
coefficient C, and is measured in Brewsters (10~'* egs). A similar bire- 
fringence occurs during the flow of liquids and is usually termed streaming 
or flow birefringence. It has been widely ascribed to the orientation and 
deformation of the coiled molecules in a polymer solution. In this review, 
a molecular origin of the birefringence will not be discussed, but the con- 
tinuum-mechanical aspects of the problem will be described in order to 
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Fig. 1. Flow-birefringence experiments (An and x) as well as flow curve (7) for a 1 
vol. % solution of polystyrene (M = 9 X 10°) in bromoform. Flow curve recalculated 


for concentric cylinder device from capillary data. re: = 88; right curve An versus Ap. 


emphasize the possibility of applying stress-optical analysis to liquids, 
whereas up to now it has predominantly been used for solids. 

The main difficulty is that the stressing of a liquid is a three dimensional 
phenomenon, distinct from the case of solids, in which the stress is in one 
direction only. In order to overcome the difficulties the stressing is usually 
performed experimentally in two dimensions (comparable to the method 
used in solid models) with the third dimension of the sample in the 
direction of the beam of light being very large compared to the 
others. The description of the optical phenomenon has generally been 
made using the tensor of polarizability of the material. The difference of 
the principal axes of this tensor determines a quantity called the degree of 
birefringence, An and the orientation of the principal axis of the polariza- 
bility tensor towards the direction of flow is called the extinction angle x. 
With solids the direction of the principal axis of stress is given directly by 
the direction of the applied stress, but with the flowing solution, this is not 
generally the case. To obtain the direction of the principal tensile stress 
in flowing viscoelastic liquids Lodge‘ suggested the application of the 
stress-optical law, used hitherto only for solids. The direction of the prin- 
cipal stress is therefore coincident with that of the polarizability tensor. 
Both are inclined to the direction of flow .at the extinction angle x which 
The major axis of the tensor of the rate of 


changes with the rate of shear. 
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(a) CD) 
Fig. 2. (a) Flow-birefringence experiments and flow curve for 3% solution of Vistanex 
B-100 polyisobutylene in matching solvents: 6-Decalin 10-White oil of 150 ep; (6) plot 
of An versus Ap for experiments in Figure 2a. 


deformation in liquids, however, is always inclined to the direction of flow 
by an angle of 45°. There exists therefore between the axes of these tensors 
a finite angle, as was emphasized by Weissenberg® some years ago. A 
derivation showing that the direction of the principal axis of deformation 
is inclined under the same angle x has been formulated using the theory of 
finite elastic deformations. This, of course, requires that the streaming 
solutions acquire elastic properties. The fact that a polymer solution does 
indeed have an elastic energy in flow is not widely accepted even though it 
was suggested many years ago. A physical explanation for its origin has 
been suggested by Mooney: relaxation of a viscoelastic material occurs at 
a rate dependent on the material, and when it is stressed at a rate deter- 
mined by the instrument this leads to a condition of elastic stress. The 
resultant elastic state of stress can then be described by the use of the 
theory of finite elastic deformations of elastic bodies. 

The investigation of polymer solutions using this general idea leads to (a) 
the determination of the relation between An and the difference in principal 
stresses Ap and (b) the correlation of the extinction angle x with the degree 
of elastic strain. 

The application of the stress-optical law requires a proportionality of An 
and Ap for a flowing solution, Ap being equal to 27/sin 2x, where 7 is the 
shear stress applied. ‘This has been proved in very numerous experiments, 
as shown in Figures 1 and 2 for polystyrene and polyisobutylene® solutions. 
It is quite important to note that in the case of polymer solutions 7, x, and 
An can all be nonlinear functions of the rate of shear D; but in all cases the 
proportionality of An and Ap, (i.e., a constant stress optical coefficient C) 
has been obtained. This proportionality is valid over 3-4 decades of Ap 
or An. It must be emphasized that the D — An curves, dependent on 
concentration, molecular weight, and solvent for one and the same polymer, 
may be linear, rising over-proportionally, or saturating to a limit in the 
range of investigation. This statement is important, for conclusions as to 
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TABLE I 


Stress-Optical Coefficients C, in Brewsters 


Polymer Solvent C 
Polyisobutylene Decalin + 1600 
Polyvinylacetate Chloroform + 1000 
Polymethylmethacrylate Benzene + 1000 

_ +2200 


Polyethylene 
Natural rubber 
Liquid polyisobutylene 


— +2000 
— + 1000-1500 


Polybutadiene Toluene +3500 
Isotactic polystyrene Bromoform — 10300 
Atactie polystyrene Bromoform — 6800 
Poly p-chlorostyrene Bromoform — 12000 
Poly 3-4 dichlorostyrene Bromoform — 12900 
Poly 2-5 dichlorostyrene Bromoform — 6900 
Poly 2-5 dimethylstyrene Bromoform —8100 





the softness or rigidity of dissolved molecules have been made from these 
curves. Sometimes, as seen in Figure 1, two influences (for example x 
and 7) compensate each other so that An is proportional to D; when the 
non-Newtonian behavior (r — D curve) is very pronounced a saturation 
is observed; when the solution is nearly Newtonian a change in x causes an 
over-proportional rise of An with D (see 3% Vistanex in 10 Fig. 2a). 

The equation of Kuhn and Kuhn? for the molecular origin of the flow 
birefringence predicts, and experiments have confirmed, that for a given 
combination of polymer and solvent, C is independent of molecular weight, 
concentration, and rate of shear within the wide range cited above, and 
the only important property is the structure of the monomer unit. How- 
ever, a restriction is imposed in that these simple conditions are valid only 
when the index of refraction of the polymer is essentially equal to that of the 
solvent, a condition usually termed matching solvent. Values for C of a 
rariety of polymers are shown in Table I. / 

The second step, namely the correlation of the extinction angle x with 
the degree of elastic strain, has been achieved by the equation s = 2 
cot 2x where s is the so-called recoverable shear.4 The deformation of 
a spherical volume element of an elastic body in simple shear (similar to 
laminar flow) to an ellipsoid in a state of elastic strain (affine transform) 
leads to a difference of the principal axis of strain: s = A — (1/A), or 
A= 8/2+ V 32/4 + 1 = cot x where J is the ratio of the major axis of 
the ellipsoid to the diameter of the initial sphere. It is obvious that as 
soon as the extinction angle yx differs from 45°—the value observed for 
normal low-molecular-weight liquids—an elastic energy appears. The 
proof that this is indeed the case is shown in Figure 3' where the recover- 
able shear s for two solutions is measured with three entirely different 
methods giving it an unquestionable physical reality: (1) the flow bire- 
fringence—using the extinction angle; (2) the Weissenberg Rheogoniome- 
ter, using the relations P;,; — Px». = rs where Py, — Ps. is the difference of 
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Fig. 3. Recoverable shear s = f (r) for a 1% solution of nitrocellulose in butylacetate and 
for 6% natural rubber in toluene measured by three different methods. 


normal stresses along the flow lines and perpendicular to them; (3) the 
direct measurement of s by the so-called recoil in a rotational viscometer. 
It is seen that s = f(r) is a nonlinear function beginning as a linear relation 
and saturating at higher values of r. Many more experimental cases exist 
where the first two methods were employed in wide ranges of polymers and 
concentrations, rates of shear, etc., with excellent agreement between the 
two. 

The investigation of a stress distribution in flow using liquids gives, in 
certain respects, more experimental possibilities than the investigation of 
solids using the stress-optical law. 

The occurrence of a polarizability tensor in the flowing solution, whose 
main axis is inclined by the angle x toward the direction of flow leads 
theoretically to a birefringence to be observed in the three planes:® in the 
1-2 plane (the one usually observed), the birefringence Any: and the angle x12 
are measured; in the 1-3 plane Any = Any cos 2x2 and x13 = 0; in 
the 2-3, Ano; and x23 are both zero. Experiments made in observing the 
birefringence between two transparent rotating parallel plates indeed 
showed, '* see Figure 4, that not only does the expected effect (Anj;) occur 
but its value is equal to the value precalculated by means of the results ob- 
tained in the 1-2 plane, i.e., from Any: and x1». 

In the theory of the mechanical behavior of polymers the symmetry 
present in laminar flow leads to an equation P22. = P;;'' meaning that the 
stresses normal to the direction of the flow lines are equal. An experimen- 
tal proof of this relation has been carried out in a fairly complicated instru- 
ment by Roberts.'* Experiments have also been carried out in a capil- 
lary,'* where the flow birefringence looking axially down the capillary in 
the direction of flow An.; was observed. The mentioned equality of normal 
stresses must lead according to the “stress optical law” to an optically 
isotropic material, i.e., Ans; = 0. This was indeed the case in all solutions 
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Fig. 4. Flow-birefringence experiments in the 1-2 and 1-3 planes; 30% polystyrene in 
toluene. 


used, whereas the birefringence in the other two planes was considerable, 
meaning that Ps. — P33 is in any case less than 0.1% of Pi. — P22». In other 
cases this method gives indeed observable results: solutions or suspensions 
of solid particles gave a very noticeable optical anisotropy which is ob- 
servable as a cross of isoclines, giving a qualitative distinction between soft 
and hard particles. These last two types of experiments have definitely 
shown that the usual flow-birefringence An, adequately describes the 
polarizability tensor in laminar flow and therefore can alone be considered. 

All these relations have been derived for a solution in which the solvent 
does not contribute to the viscosity or the flow birefringence. When one 
goes to more dilute solutions, in which the viscosity is commensurable with 
that of the solvent and the contribution of the birefringence of the solvent 
is not negligible, the conditions become considerably more complicated but 
can be analyzed (as has been done in a series of publications’). 

Reviewing the whole field one sees the flow birefringence of solutions 
of a viscoelastic material (such as polymer solutions) can be described by 
the following tensors: (a) the rate-of-shear tensor, whose principal axis is 
45° to the flow direction; (6) the stress tensor, whose shear component 
is directed in the direction of flow and whose principal tensile 
axis is directed at the angle x to the direction of flow; (c) the tensor of 
elastic response or elastic deformation, whose shear component (the re- 
coverable shear s) is the direction of flow and whose principal axis is in- 
clined at the same angle x as the stress tensor; (d) the tensor of polar- 
izability, whose principal axis is also inclined by the angle x to the direction 
of flow. In the theory of the flow birefringence of coiled molecules, be- 
ginning with Kuhn and Griin'‘ the shape of the molecule, deformed in 
laminar flow, is used to calculate the polarizability tensor of the solution, 
introducing as a constant the difference of polarizabilities of the statistical 
element. However, the ratios of axes of these tensors vary widely. The 
polarizability tensor is nearly spherical because the birefringence is very 
small as compared to the mean index of refraction or the change in polar- 
izability is negligible as compared to the mean polarizability. The elastic 
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deformation tensor, as well as the shape of the deformed sphere it describes, 
deviate strongly from the spherical shape. The ratio of the principal axes 
of the tensor is: p = A* = cot? x, or for large values of s: p + s%. The 
values of p can be quite large, as the following calculation shows: for s 
= 1, x = 31°43’ and p = 2.62; fors = 2, x = 22°30’, and p = 5.82. 
In Figure 1, the smallest x measured was 3°, which leads to s = 19 and 
p = 360. Even under such extremely large deformations, as can be seen 
in Figure 1, the stress-optical law remains valid. From Figure 2 it is seen 
that Apmax ~ 10° dynes/cm.? has been experimentally achieved, which 
corresponds to a tensile stress in the polymer (3%) of 3.3 kg./em.”. For 
rubbers (as quoted in Treloar’s book) stresses up to 50 and 85 kg./cm.? 
have been used, still complying to the linear An versus Ap relation. These 
facts emphasize the generality of this law. The mentioned equality of the 
recoverable shear s as measured by three different methods allows one to 
measure the normal stress P;, and Px» using flow birefringence alone.® 
If the C is known, as is the case, and the relation P22. = P33 is used, as has 
been experimentally proved above, then: 


(Pu —_ P 9) C = Any = Any COs 2x 


When C is not known, a measurement of the viscosity at the same rate of 
shear gives 7 and: 


Py — Px» = rs = 2r cot 2x 


Further the analysis of laminar flow has been attempted in still another way. 
Gill® suggested the use of a plastic tube in an elliptical guide. The tube 
and its liquid contents are rotated and the optical phenomena are ob- 
served by axially transmitting polarized light. This approximates the 
stressing with rotating axes and gives rise to a cross of isoclines determining 
a birefringence An and the extinction angle x. A calibration with New- 
tonian liquids'® (liquid polyisobutylenes) gave an exact quantitative cor- 
relation with the usually measured flow-birefringence. In viscoelastic 
liquids the x diminished with rate as predicted, but due to experimental 
difficulties, no exact correlation was as yet possible. Such a stressing is 
similar to a dynamic experiment with w = D. Then x = 1/2 6, when 
6 is the loss angle and tan 6 = @G”/G’ calculated from the complex modulus 
7* = G’ + jG". The stress optical behavior of polymer solutions is, 
as far as could be ascertained, exactly parallei—even quantitatively so— 
to the one of rubberlike polymers or elastomers. For these materials, as 
described, An is determined by the stress. Therefore, in a dynamic ex- 
periment the phase of An and 7 should be expected to be identical—this is 
indeed the case as subsequent results presented at this symposium have 
shown. Concerning the suspensions of rigid particles—such as the micellar 
solution of milling yellow—and crystalline polymers no such prediction is 
possible—the experimental results are the only answer. 
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Résumé 


La biréfringence d’écoulement consiste dans l’anisotropie optique induite par un flux 
laminaire dans les liquides. On décrit l’anisotropie optique par un tenseur optique qui 
est coaxial avec le tenseur d’élongation dans un matériel élastique. Ceci est évident 
dans une expérience de tension. En flux laminaire, ceci a été prouvé en étudiant la 
distribution de la biréfringence dans trois directions de l’espace. La quantité d’aniso- 
tropie, habituellement décrite par le degré de biréfringence An est proportionnelle a la 
diffusion de tension principale du tenseur d’élongation mécanique dans un trés grand 
domaine de variables, donnant par 14 la possibilité d’utiliser les mesures optiques pour 
la détermination des tensions mécaniques. Cependant, la relation est yalable pour les 
solutions concentrées et pour des conditions ot le solvant et le polymére ont essentielle- 
ment le méme indice de réfraction n—, condition habituellement appelée solvant de 
comparison. La constante de préportionnalité entre la différence en tension principale 
et An—le coefficient de tension optique C—est, comme prévu par la théorie, indépend- 
ante des concentrations, température de tension de cisaillement, poids moléculaire et sa 
distribution, et uniquement déterminée par la structure de l’unité monom-rique et son n. 
Une autre quantité en biréfringence d’écoulement est l’angle d’extinction qui est mesuré 
entre la direction de l’axe optique principal et la direction du flux. On a obtenu un 
nombre suffisant de résultats expérimentaux pour montrer qu’un angle peut étre quanti- 
tativement mis en corrélation avec le cisaillement récupérable, qui est le composant de 
cisaillement du tenseur d’élongation élastique en flux laminaire. De plus, la détermina- 
tion expérimentale de |’angle donne une mesure de la tension élastique. Une étude de 
résultats expérimentaux obtenus dans ce domaine montre les généralités oe cette ap- 
proximation dans une variété de systemes de polymére et solvants. 
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Zusammenfassung 


Strémungs- oder Fliessdoppelbrechung ist die optische Anisotropie, die durch laminare 
Strémung in Fliissigkeiten verursacht wird. Die optische Anisotropie kann durch einen 
optischen Tensor, der zum Verformungstensor in elastischem Material koaxial ist, be- 
schrieben werden. Dies wird in einem Spannungsexperiment offensichtlich. Bei Lam- 
inarstr6mung wurde dies durch Untersuchung der Doppelbrechungsverteilung in den 
drei Richtungen das Raumes bewiesen. Der Anteil der Anisotropie, der gewéhnlich 
durch den Grad der Doppelbrechung An ausgedriickt wird, ist der Differenz der Haupt- 
spannungen des mechanischen Spannungstensors in einem sehr grossen Variablenbereich 
proportional, und gibt daher die Méglichkeit, optische Messungen zur Bestimmung 
mechanischer Spannungen zu verwenden. Die Beziehung ist entweder fiir konzentrierte 
Lésungen oder fiir Bedingungen, bei denen das Lisungsmittel und das Polymere im 
wesentlichen denselben Brechungsindex n haben, eine Bedingung, die gewshnlich 
“Lésungsmittel-Angleichung’ genannt wird, giiltig, Der Proportionalitatsfaktor 
zwischen der Hauptspannungsdifferenz und An—der spannungs-optische Koeffizient C— 
ist, wie auch aus der Theorie folgt, unabhangig von der Konzentration, Schubspannung, 
Temperatur, Molekulargewicht und seiner Verteilung, und ist nur durch die Struktur der 
Monomereinheit und ihres n bestimmt. Eine weiter Grisse bei der Strémungsdoppel- 
brechung ist der sogenannte Ausléschungswinkel, der zwischen der Richtung der op- 
tischen Hauptachse und der Fliessrichtung gemessen wird. Eine betrichtliche Anzahl 
experimenteller Daten zeigt, dass dieser Winkel quantitativ zu der sogenannten ‘‘re- 
versiblen Scherung,”’ nimlich der Schubkomponente des elastischen Verformungstensor 
bei laminarem Fliessen in Beziehung gesetzt werden kann. Daher gibt die experi- 
mentelle Bestimmung des Winkels ein Mass fiir die elastische Verformung. Ein Uber- 
blick iiber die experimentellen Ergebnisse zeigt die allgemeinen Merkmale dieser Methode 
bei einer Reihe von Polymer- und Lésungsmittelsystemen. 
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Streaming Birefringence as a Rheological 
Research Tool 


HAROLD WAYLAND, Division of Engineering and Applied Science, 
California Institute of Technology, Pasadena, California 


Synopsis 
This paper reviews the theoretical and experimental evidence currently available to es- 
tablish criteria to aid the experimenter who wishes to use streaming birefringence as a tool 
for studying fluid flow. Instrumental techniques are reviewed, and the possibilities for 
both qualitative and quantitative flow visualization by streaming birefringence in pure 


Newtonian liquids, colloidal suspensions and dilute and concentrated polymer solutions 


are discussed. Methods of reconstruction of the velocity and stress. fields from 


birefringence data are outlined. 


I. Introduction 


If a pencil of plane polarized light is passed through a flowing liquid 
which exhibits streaming birefringence (SBR), two pieces of experimental 
information can, in general, be obtained: (/) the orientation (with 
respect to some preassigned set of axes) of the isoclinic axes, and (2) the 
amount of birefringence (i.e., the relative shift in phase between a beam 
polarized along one of the isoclinic axes to a beam polarized along the 
other. This retardation is usually measured as the fraction of a wave- 
length retardation per centimeter of path length for the particular wave- 
length used, although theoretically it is more apt to be expressed as the 
ratio of the index of refraction difference for light polarized along these 
two axes to the mean index of refraction of the medium for this wave- 
length.) The relationship between the direction of the isoclinic axes and 
the streamlines of the flow will depend on the mechanism which gives rise 
to the double refraction. It is the purpose of this paper to review the 
theoretical and experimental evidence currently available to try to es- 
tablish useful criteria to aid the experimenter who wishes to use SBR as a 
tool for studying fluid flow. In particular it is important to ascertain 
whether or not quantitative interpretation is possible and to set criteria 
to permit the choice of the method of interpretation applicable to the 
system under consideration. 

The fact that, in shear flow, certain liquids exhibit different indices of 
refraction for light polarized in different directions was first observed by 
Clerk Maxwell' in 1866. The first published suggestion of its use as a 
flow visualization tool followed a demonstration by Humphrey? at the 

il 
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London Physical Society in 1923. It is interesting to note that at this 
demonstration one of the audience raised the warning that it would be 
unwise to use Dr. Coker’s method of analysis of photoelastic systems for 
interpreting flowing systems. This warning has too often passed un- 
heeded! In spite of early interest, SBR became a useful quaiitative tool 
for flow studies only during World War II, and a solid basis for quantitative 
studies—both theoretical and instrumental—has been laid only during 
the past eight or ten years. 

Except for the early work of Sadron* and Alcock and Sadron,‘ which will 
be discussed in detail in Sections II and III, attempts at systematic in- 
terpretation of flow patterns by means of SBR seem to have first been 
published by Hauser and Dewey.’ They used bentonite solutions in 
concentrations of the order of 1%. Dewey went into some detail on 
possible quantitative interpretation in his doctoral thesis,’ but lacking a 
suitable theoretical basis his results were inconclusive. Hauser and 
Dewey’s work on use of SBR for studies of flow separation and transition 
to turbulence did lead to some useful results. They made both still and 
motion pictures of the flow around two-dimensional models of automobiles, 
although this appears to have been done as a demonstration of the method 
rather than for design purposes. Taking advantage of personal advice 
from Dewey, Knapp and his co-workers*:’ at the California Institute of 
Technology developed a polarized light flume early in World War II. 
This used a 0.2% solution of white Hector bentonite for the optically 
active medium. It was used to study flow separation at the nose and tail 
of underwater weapons, such as torpedoes and depth charges. The 
models were three dimensional, the fluid was not strictly Newtonian, and 
essentially nothing was known about how to make a quantitative inter- 
pretation even of two-dimensional data. But the method was extremely 
fruitful in improving body shapes on a cut and try basis. The work, 
unfortunately, received little publicity outside of military circles since it 
carried security classification at the time. 

In 1945 Leaf" used the SBR of a bentonite solution to study flow sep- 
aration in a firebox and was successful in improving the design with 
the aid of this method. In the same year Binnie" published his first 
paper on studies of the onset of turbulence in pipe flow using a solution 
of the dye, benzopurpurin. The use of SBR techniques has contributed 
considerably to the studies of the onset of turbulence in pipe flow in the 
hands of Binnie,'' Binnie and Fowler,'? and Lindgren,'* using visual, 
photographic, and photocell observation of the presence and magnitude 
of turbulence and turbulent flashes. Lindgren found that a concentration 
of bentonite as small as 0.062% gave a measurable influence on the velocity 
of a turbulent slug. Although it is fair to say that the usefulness of SBR 
of colloidal solutions for studying flow separation and the onset of tur- 
bulence is well established, Lindgren’s work demonstrates that great care 
must be exercised to avoid effects due to the presence of the colloid when 
fundamental fow mechanisms are being studied. 
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As far as I know, SBR techniques have not been used with visco- 
elastic liquids to study qualitative differences in their behavior vis a vis 
Newtonian liquids. Vor example, such techniques should be useful in 
observing the type of secondary flow predicted by Ericksen'*:” in the 
flow of viscoelastic fluids in elliptic cylinders. Tor this type of observa- 
tion a precise knowledge of the mechanism giving rise to the SBR is not 
necessary. 

lor quantitative interpretation three types of systems will be considered: 
pure liquids with Newtonian viscosity behavior, colloidal solutions and 
solutions of polymers. The physical mechanism giving rise to SBR in 
ach case is quite different, and it is only since adequate theories have been 
developed and experimentally verified that it has been possible to put 
flow visualization by SBR on a sound footing. The state of the theoretical 
knowledge on SBR as of about 1958 is summarized in the excellent review 
article by Jerrard.'® 

It should be made clear in the beginning that quantitative studies will 
be limited to two-dimensional flow situations. It is theoretically possible 
to predict the amount of SBR associatea with a given shear rate or stress 
pattern (depending on the appropriate mechanism) for complicated 
flows. However, the tensor nature of the optical properties of a bire- 
fringent system makes the inverse problem extremely difficult, if not 
impossible. It is not clear that such a solution is unique even for cylin- 
drically symmetric flows such as Poiseuille flow. 


II. Instrumental Techniques 
1. Photographic Methods 


The simplest method of collecting data on the optical properties of a 
two-dimensional flow field is to use the well established photographic 
techniques employed in viscoelastic studies in solids. These are well 
documented in the literature, as in the classic work by Frocht,” and will 
only be briefly sketched here. The conventional physical arrangement 
is shown in Figure 1. l’irst let us consider the case in which monochromatic 
light from the source is rendered parallel by the lens, plane polarized by 
the polarizer, passed through the flow field and viewed through the an- 
alyzer crossed to the polarizer (with the quarter wave plates out of the 
system). In this case there will be no light transmitted through those 
portions of the flow (a) in which the retardation between the fast and slow 
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Fig. 1. Flow polarimeter. 
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Fig. 2. Types of circular polariscopes: (a) crossed; (b) parallel; (¢) mixed. 


axes is an integral multiple of the wavelength of light used (called csochro- 
matics) and (6) in which the plane of polarization of the incoming light 
is parallel to either of the isoclinic directions (called zsoclinics). If the 
polarizer and analyzer are rotated together, coupled in the crossed position, 
the isochromatic pattern will remain fixed, while the isoclinie will, in 
general, vary. A series of isoclinic maps can be taken by systematic 


rotation of the polarizer-analyzer combination. 
Although the isochromatics can be obtained from the configuration just 
discussed, they are more usually obtained by passing circularly polarized 
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light through the test region. This is obtained by introducing a quarter 
wave plate just after the polarizer, with axes at 45° to the plane of polariza- 
tion. A second quarter wave plate follows the test region. If its fast 
axis is parallel to the slow axis of the first \/4 plate and the polarizer is 
crossed with respect to the analyzer we have the crosscd circular po- 
lariscope (Vig. 2a). If the fast axes of the \/4 plates are parallel and the 
polarizer and analyzer are parallel, we have the parallel circular polariscope 
(Fig. 2b). In either of these cases extinction will occur wherever the retarda- 
tion is an integral number of wavelengths. Ifthe quarter wave plates remain 
parallel, but the analyzer is crossed with respect to the polarizer, we have 
the mixed circular polariscope (Vig. 2c). In this case extinction occurs 
when the retardation is an odd number of half wavelengths. 

These methods, using monochromatic light, are primarily useful when 
the maximum retardation in the field is several wavelengths so that several 
fringes are available to use in analyzing details of the flow pattern. Frac- 
tional fringes can be more easily studied if white light is used, but the 
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Fig. 3. Sadron’s spectographic flow polarimeter. 


circular polariscope loses some of its quantitative character since a quarter 
wave plate can give strictly circularly polarized light at only one 
wavelength. 

Sadron’ developed a spectroscopic technique to permit analysis of small 
differences of retardation. He was dealing with a system (plane Poiseuille 
flow in sesame oil) in which the isoclinies were at 45° to the streamlines. 
A beam of plane polarized light with the plane of polarization parallel 
to the streamlines was passed through a full wave retardation plate, cut 
to give one wavelength retardation at 5800 A., with optic axes at 45° to 
the plane of polarization (Fig. 3). The analyzer was crossed with respect 
to the polarizer. The light from the analyzer was photographed through 
a spectroscope. With no flow, extinction would occur at \)» = 5800 A. 
across the field. A retardation introduced by the flow would shift the 
extinction position, and the amount of retardation was correlated with 
position across the flow. This method can be made extremely sensitive, 
but requires knowledge of the isoclinie directions. Except for the early 
work of Sadron* and Alcock and Sadron,‘ this method does not seem to 
have been exploited, 
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2. Point by Point Measurements 


Point by point scanning of the flow field is particularly desirable when 
the amount of birefringence is small. This has been developed by Way- 
land and Sutera and Wayland’ to permit obtaining optical data to 
within 0.1 mm. of a wall. To obtain a narrow pencil of light the optical 
arrangement of Bjérnstahl” is the most suitable. The pinhole nearest 
the light source (Fig. 4) is imaged at the lower window of the flow channel 
and the second pinhole is imaged at the upper window of the flow channel. 
The diameter of the pinholes can be so chosen that their images have 
equal diameters. It has proved possible to work with images such that the 
diameter of the first dark diffraction zone is 0.2 mm. Since monochro- 
matic light is used, only a simple lens is necessary. The beam can be 
translated parallel to itself by means of a plane parallel sheet of glass 
(the deflector plate in Fig. 4) which can be rotated about an axis per- 
pendicular to the axis of the light beam. The thickness and index of 
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Fig. 4. Optical system for point-by-point scanning 
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refraction of the glass can be chosen to give the required amount of transla- 
tion, although the system is practical only if the total translation is limited 
to a few millimeters. 

Various methods are available to analyze the effect of the birefringent 
medium on this polarized pencil of light. Since this point by point method 
seems particularly suitable for the study of systems with small retardation, 
a brief resume will be given of some useful methods for this case. The 
most difficult measurement is that of the angle of isocline. With crossed 
polarizer and analyzer the intensity of light passing the analyzer for a 
retardation 6 when one of the axes of the birefringent medium makes an 
angle of 6 with plane of polarization (see Fig. 5a) is 


I = sin? 28(1 —cos 8) (1) 


If the retardation is small, and 8 is approaching the extinction angle, 
we have 


I = B75? (2 
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Fig. 5. Arrangements for locating the angle of isocline: (a) location with crossed polarizer 


and analyzer; (b) location with modulated optical rotation. 
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It is seen that in this case the null is proportional to the square of both 
the angular position and the amount of birefringence. A variety of tech- 
niques has been developed to increase the sensitivity of location of the 
angle of isocline. Of these, only the method using a modulated optical 
rotation will be discussed in detail. Certain visual half-shadow methods 
give linear dependence of intensity with angular position and amount of 
birefringence, but lack the advantage of ac photoelectric detection. The 
reiationship of these methods to the modulated methods is discussed by 
Wayland and Badoz.”! 

Referring to Figure 5b, monochromatic light polarized parallel to 
the x, axis passes through the birefringent medium M of retardation 6 
with one optic axis at the angle 6 to the x; axis. The light then passes 
through a quarter wave plate Q oriented as shown, and then through an 
optical rotator 7 which rotates the plane of polarization sinusoidally 
through +7. This modulation can readily be accomplished by passing 
the light through a Faraday cell in which the magnetic field is actuated by 
an alternating current. The light then passes through the analyzer A, 
crossed with respect to the polarizer. If the light falling on the polarizer 
is of unit intensity, the intensity of light passing the analyzer will be given 
by the expression (when the rotator gives the fixed rotation y) 


I = '/,[1 — cos 2y(cos? 28 + sin? 26 cos 6) —sin 2y sin 28 siné] (3) 


If y = yo sin wt and yo is small, this expression consists of a continuous 
term, a term at the driving frequency and a term at double the driving 
frequency. If the amplifying system on the photomultiplier current is 
tuned to amplify the fundamental and to reject de and the second harmonic, 
the detected signal will be given by the expression 


S = '/;sin 2y sin 26 sind = 6Byo sin wt (4) 


which is linear in both the retardation 6 and the angular displacement from 


the nuli point 8. 

To determine the amount of birefringence with this arrangement, 
the entire assembly of polarizer, quarter-wave plate and analyzer is rotated 
through 45° from the isocline 8 of the medium. If the analyzer is now 
uncoupled from the system it can be rotated independently to seek a new 
null position. If it is rotated an angle a from being crossed with respect to 
the polarizer, the intensity of light transmitted by the analyzer is 


I = '/, [1 — cos 2a(cos 2y cos 6 — sin 27 sin 8) 
— sin 2a(sin 2y cos 6 + cos 2y sin 6)] (5) 
Modulating the rotation and detecting only the fundamental we get 
I = yo sin wtl[eos 2a sin 6 — sin 2a cos 6] = yo sin wt sin (6 — 2a) (6) 


This vanishes when a = (6/2) and the change in signal is linear in angular 
unbalance when near the null point. 
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In practical application the angle of isocline can readily be determined 
to 1° for a relative retardation of 5 X 10~-*° cm. using the mercury green 
line. The least count of the measurement of retardation in the apparatus 
constructed at the California Institute of Technology was 0.3 A. For 
extremely small birefringences the precision was increased by measuring 
the retardation as a function of angular position over about 200° of rotation 
and obtaining both the amount of birefringence and the isoclinic positions 
from a smooth curve drawn through these data points. This limitation 
was not due to the lack of sensitivity and stability of the detecting system, 
but to limitations in measuring the angular position of the analyzer 
(+0.01°) and in extraneous birefringence in the windows. A detailed 
discussion of the experimental arrangements and the correction of errors 
will be found in the doctoral thesis of Intaglietta,?* and will be published 
in detail elsewhere. 


lil. Pure Newtonian Liquids 


Certain oils and organic liquids show SBR well within the range of 
Newtonian viscosity behavior. As early as 1928 Raman and Krishnan** 
developed a theory of SBR for pure liquids. This was based upon two 
experimental facts: (/) an extensive series of measurements on organic 
liquids by Vorlinder and Walter?‘ appeared to show that all liquids 
showing measurable SBR possessed highly asymmetric molecules, and 
(2) measurements in simple shear flow showed that for these pure liquids 
the optic axes were at 45° to the streamlines, i.e., were parallel to both 
the principal stress and strain rate axes (since these are coincident for a 
Newtonian fluid). Raman and Krishnan assumed that such asymmetric 
molecules would tend to be oriented by the stresses in the fluid in such a 
manner that the long axis of the molecule would be along the axis of 
principal tension. Thermal agitation would give rise to Brownian motion, 
tending to disorient the molecules. Classical statistical mechanics was 
applied to calculate the resulting equilibrium distribution for a given flow 
condition, giving the relationship between the local stress field and the 
resulting birefringence. 

Sadron*® later threw considerable doubt on the validity of the Raman- 
Krishnan theory, and the Vorlinder-Walter data were proved to be 
virtually useless since many of the measurements were made in non- 
laminar flow which Jerrard® and Wayland" have shown can lead to com- 
pletely misleading results. 

Peterlin and Stuart”” have developed a theory of SBR of pure liquids in 
which the flow orientation of the asymmetric molecules is included as well 
as certain internal field corrections due to molecular interactions. This 
predicts that the angle of isocline will deviate from coincidence with the 
principal strain rate axes for high enough velocity gradients, but for most 
organic liquids with essentially Newtonian viscosity this requires such 
high gradients that thermal effects make measurements extremely dif- 
ficult. The definitive experiment distinguishing between these two theories 
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is reported by Tolstoi.** The birefringence induced in purified transformer 
oil by a combination of shear flow and an electrostatic field was compared 
with theoretical treatments based on the statistical mechanical theory 
and the flow orientation theory. The results gave strong support to the 
flow orientation theory of Peterlin and Stuart. 

Both theories predict that at moderate shear rates in simple shear 
flow the optic axes will remain parallel to the principal stress and strain 
rate axes. At very high shear rates the Peterlin-Stuart theory predicts 
a departure from coincidence of the optic and strain rate axes, but in the 
pure liquids studied to date this departure does not show up until the 
shear rate exceeds 10‘ sec.—'. No critical work has been done on the re- 
lationship of the SBR to the stress or strain rate in flow other than of the 
Couette type. Since the relaxation time for the angular rotation of the 
molecules of these liquids is very short, it is reasonable to expect that the 
optic axes will coincide with the principal stress and strain rate axes. 
This means that if the liquid is calibrated in Couette flow, a measurement 
of the direction of the isoclinic axes and the amount of birefringence at a 
network of points in the flow will permit a determination of the principal 
stresses and strain rates from the optical measurements. [rom these 
data the flow field can be reconstructed if it is known along some boundary 
or streamline, as will be discussed in Section IV. 

The earliest reported attempt at quantitative flow measurements by 
SBR was made with a Newtonian liquid, sesame oil. This was the work 
of Sadron* and Alcock and Sadron‘ using the spectroscopic method dis- 
cussed in Section II.1. Their measurements on flow between parallel 
plates agreed with the calculated flow within about 4%, with most of the 
error occurring close to the walls where strains in the windows near the 
clamped edge may be large and spurious birefringence due to reflections 
can introduce significant errors. These latter errors can be reduced by 
‘are in maintaining the parallel nature of the light beam (now possible 
with a laser light source) and by covering the wall with a dielectric of 
index of refraction less than that of the flowing liquid so that light reflected 
from the walls at grazing incidence will be totally reflected.2® In the case of 
diverging flow the agreement was not so good. Their published data 
are not sufficiently detailed to ascertain whether this lack of agreement is 
due to extraneous birefringence, failure to take into account the variation 
in the isoclinie direction as the flow is traversed, or something more fun- 
damental. Critical experiments are still needed to clarify the validity 
of the assumptions concerning the relationship between the SBR of these 
Newtonian liquids and the stress and strain rate fields in converging and 
diverging flow situations. The specific birefringence of such liquids is 


generally quite small—ethy] cinnamate, one of the more active substances, 


has a specific birefringence of 5.04 K 10~-* (A./em.) sec. A careful set of 
experiments in a complicated, but calculable, flow using point by point 
measurements would be helpful in clarifying the potential role of using the 
SBR of pure Newtonian liquids as a flow visualization tool. 
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Point by point studies have been made with ethyl cinnamate in cy- 
lindrical Couette flow above the Taylor transition velocity.’* Experi- 
mentally it was found that the isoclinics remained at 45° to the mean 
streamline direction even in the turbulent portion of the flow. Assuming 
that the amount of birefringence was proportional to the mean velocity 
gradient, and using the calibration made in laminar flow, the velocity 
profile was reconstructed by numerical integration from the zero velocity 
at the fixed outer wall. With the velocity of the moving wall 168 cm./ 
sec. integration gave 180 cm./sec. or 7% high; for the wall at 252 cm./ 
sec. integration gave 249 cm./sec. or 3% low; and for the wall at 492 cm. 
sec., integration gave 476 cm./sec. or 3% low. This result is, at best, 
empirical, since there is no adequate theoretical treatment for SBR in 
a turbulent flow. It is, however, an excellent illustration of the pos- 
sibilities in the point by point method: for studying boundary layers. 
The measurements were made in a 1 mm. gap, and the basic data for the 
highest rotative speed are shown in Figure 6a. Squire® compared these 
data with his theoretical treatment of turbulent cylindrical Couette 
flow. Introducing the dimensionless velocity V, related to the linear 
velocity of the wall of the inner cylinder, he plotted the velocity gradient 
deduced from the birefringence data against the reciprocal of the dimension- 
less velocity from the fixed wall to the center of the gap and against 1/ 
(1-V) from the moving wall to the center of the gap (Fig. 6b). The hori- 
zontal dotted portion of the curve corresponds to a laminar sublayer 
(constant velocity gradient) near the walls, and the sloping linear portion 
is consistent with the behavior predicted from his theoretical treatment in 
turbulent flow. 


IV. Suspensions of Rigid Particles 


Most engineering work on flow visualization with SBR has been done 
with colloidal solutions. Although a considerable body of theory is 
available on which to base flow interpretation from the optical data, all 
of this theory applies only to situations in which interactions among 
particles are neglected. There is a great temptation to increase particle 
concentration to such a point as to give large optical effects, but in every 
case which I have reviewed in which the particle concentration exceeds 
a few tenths of a per cent there is evidence that the medium no longer 
follows the laws used in interpreting the data. The theoretical treat- 
ment on which flow interpretations are based will be summarized and the 
current state of the experimental art will be discussed in terms of this 


background. 

Peterlin and Stuart*! laid the basis for a quantitative theory of the 
SBR of a dilute solution of rigid ellipsoids in Couette flow, which Way- 
land*®* has extended to general two-dimensional laminar flow. These 
treatments are based on the following assumptions: 

(1) The particles are rigid ellipsoids of revolution. 

(2) There are no hydrodynamic, chemical, or electrical interactions 
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(3) The number of particles along the length of the light beam is 
large. 

(4) There is a finite number of discrete particle sizes (although the 
extension to a continuous distribution is straightforward). 

(5) The unperturbed flow field throughout distances comparable to 
the largest particle dimension can be calculated from a knowledge of the 
velocity and the velocity gradient at any point, i.e., second and higher 
space derivatives of the velocity can be neglected over one particle length. 

(6) The linearized Navier-Stokes equations can be used in describing 
the motion of the individual particles between Brownian fluctuations, 
hence inertial terms can be neglected in describing the rotation of the 
particles. 

(7) The particles are subject to Brownian motion, and the diffusion 
law is valid. 

(8) The relaxation time for rotation of the individual particles is 
short enough that the particles can be assumed to be in equilibrium with 
the flow at every point. 

(9) The particles are short with respect to the wavelength of light 
used. 

On the basis of these assumptions the position of the isoclinic x from 
the first principal strain rate axis, and the amount of birefringence, An, 
for a monodisperse system will be given by 


E sin 2Ao E Ek? ( in? DAe + =) £0 (*) 
= — ——_— _ sin? 2/ - 
7 6D ame; ot Se D 
+, 1 (v2 . 240? :) é 
= — 1 — k2)+... (7) 
12D | 27D? ( 4° 35 se 


nes (ee Va KE? («i 294. 4 s5) 10 (=) 

a na /\D 18D? pare 35 D4] | 
se (AV FG Eis 4 6b? r*) 4 | (8) 
~ A\B/\ a D 18D? 4 3 ree 


where EF is the magnitude of the principal strain rate; Ay is the angle be- 
tween the principal strain rate axis and the streamline; 6 = (a,? —a2’) 
(a;2 + as”) with a; = major axis and a, = minor axis of ellipsoid of rev- 
olution; y, = 2H sin 2Ao, the shear rate along the streamline; D is the 
rotary diffusion coefficient for the ellipsoid, c is the volume concentration ; 
G = gi —ge, the optical anisotropy for the particle; and # is the mean 
index of refraction for the solution. 

In Couette flow Ay = 45°, hence 2h = y,. Down the centerline of a 
symmetric converging or diverging channel flow, Ao = 0, hence x = 0 
and the quadratic term in An is considerably different from that for 
Couette flow. In the case of sufficiently elongated particles, b = 1, and 
we have the following expression for An: 
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an = (“fr —o200 (8) 
oe ae 7 


On the centerline of converging or diverging channel flow 


An = 4 @t 0.038 G) | (10 
ya — V.U5 
D) D ee 


If there is a sufficiently strong shear component that Ay ~ 45° one can 
interrelate / to An as measured in Couette flow even outside the range in 
which An is linear with shear rate. lor more general flows, however, 
this correlation cannot be made without a knowledge of the angle of 
isocline and of the axial ratio parameter b. In the case of a polydisperse 
system, the expressions are considerably more complicated and, to work 
in the nonlinear range, would require an elaborate calibration. 

It seems probable that no work on flow visualization using colloidal 
solutions has met all of the assumptions listed above. Tobacco mosaic 
virus in the concentrations used by Sutera and Wayland'® probably comes 
the closest. to meeting all of these conditions, but in their concentrations 
(less than 0.39%) the effects are very small. Furthermore, the particles 
are not short compared to the wavelength of light used (particle length 
3000 A., A = 5460 A.). By increasing the viscosity approximately 80 fold 
by addition of glycerine, it was possible to increase the relaxation time 
so that the birefringence measured at a point reflected the previous flow 
history rather than giving information on the flow field at this point. 

Bentonite, one of the most popular materials, probably does not consist 
of rigid particles, and shows only a small specific birefringence at con- 
centrations low enough that interaction effects are not significant as shown 


Couette flow 


by Lindgren. '* 

The substance most recently to achieve popularity for flow birefringence 
work is the dye milling yellow. As yet I have seen no accurate char- 
acterization of its size and shape but the phenomenological data reported 
in the literature lead me to raise the following questions concerning inter- 
pretation of flow data from SBR of milling yellow in the concentrations 
used in the experiments so far reported (1'/,;-1'/2%). In the first place, 
the amount of birefringence shows a far from linear dependence on con- 
centration, which suggests strong interaction phenomena, and immediately 
vitiates application of the theory for rigid, isolated particles. Another 
indication of interaction effects is the viscometric data reported by Honey- 
eutt and Peebles** which shows a rather strong shear dependence of vis- 
cosity, stronger than seems explicable on the basis of orientation effects 
alone, and sufficient to question the use of a constant viscosity in cal- 


culating the flow fields studied to date. 
Even assuming that the milling yellow solutions obey the laws for a 
monodisperse, non-interacting colloidal solution of rigid ellipsoids, it is 
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Order of isochromatic observed in channel 





Rate of deformation, sec. ai 
Fig. 7. Data of Prados and Peebles for m illing yellow solution T-7: (——) data in 
Couette flow; (---) calculated data for centerline of converging flow. 


interesting to look at the effect of the observed nonlinearity of the An 
versus velocity gradient curve on the possible interpretation of the local 
value of #. The curve of Prados and Peebles*‘ for their solution T-7 
at 25.20°C. is given as a solid line in Figure 7. Since this was taken in 
Couette flow, the coefficients in eq. (9) were empirically determined to 
match this curve. On the centerline of converging channel flow, where 
Ao = 0, the quadratic dependence is reduced by a factor of 6.85, giving 
the dotted curve in Figure 7. It is true that the total strain rate is not 
apt to be large on the centerline of converging channel flow, but certainly 
caution must be exercised when working outside the region in which An 
is linear with strain rate. 

I feel that we can honestly raise the question as to whether the dis- 
crepancies reported by Bogue and Peebles® between optically derived 
velocity profiles and the theoretical profiles for channel flow are due to 
entrance conditions, as they suggest, or to something more fundamental, 
such as non-Newtonian character of the milling yellow solution, or the 
need to take the nonlinear behavior of An with strain rate into considera- 
tion. An even more fundamental difficulty may arise: the basic re- 
lationship between An and strain rate may not even be of the form of 
eq. (8), and hence An may depend on F in a somewhat different fashion 
when the flow is not Couette flow. This should not be surprising: the 
interactions which lead to an amount of SBR considerably larger than that 
expected from orientation of individual particles may introduce a mech- 
anism somewhat akin to network deformation such as appears to be re- 
sponsible for the SBR of certain polymer solutions (to be discussed in the 
next section) superposed on a flow orientation component. In this case 
the SBR would be related to a complicated combination of stress and 
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(a) (b) 


Fig. 8. Coordinate systems for flow analysis: (a) local Cartesian coordinates; (b) local 
streamline cocrdinates. 


strain rate—an area in which little satisfactory theoretical work has been 
done. The fact that Prados and Peebles,*4 and Bogue and Peebles® get 
reasonably good agreement with expected flow profiles seems to indicate 
that these effects are not extremely large in the situations studied, but 
they may well represent a fundamental limitation on the use of milling 
yellow, or any other colloidal suspension, in sufficient concentrations for 
photographic techniques to be applicable. 

If one is working.in the range of validity of eqs. (7) and (8) it is possible 
to reconstruct the flow field from a knowledge of the birefringence field and 
a knowledge of the flow along some streamline—usually an application of 
the no-slip condition to a boundary. Prados and Peebles*‘ give step by 
step procedures for various flow situations based on an assumed relation- 
ship between the SBR and the flow parameters. The following treatment!® 
is more directly based on the theory of SBR in generalized two-dimensional 
flow. If we are in the linear range, eqs. (7) and (8) reduce to 


xX = c; E sin 2Ay (11) 
An = ck (12) 


l'rom a knowledge of An at each point in the flow we can obtain F provided 
we obtain c. by calibration in Couette flow. Suppose we start from a 
point on a fixed boundary where we know the streamline direction and 
have measured An and the directions of the isoclinics. From An we 
calculate 2 from the Couette calibration and the deviation of the isoclinics 
from the principal axes, X, = ¢#. We know that one of the optic axes 
in the unknown flow will make an angle of Ag — xX with the streamline, 
and this we know from the position of the isoclinics. Consequently we 
solve for Ap from the transcendental equation 

(Ag —X) = Ao —X, sin 2Ao (13) 
If we take a local set of Cartesian coordinates with the x-axis in the positive 
streamline direction (Fig. 8a) we have the following relationship for an 
incompressible fluid 
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; ; 1] Ou | ov cores 
dey = bye = 3 & ok = = E sin 2Ao (15) 


Since x has been chosen to be tangent to the streamline, (0v/dx) = 0 at 
that point and eq. (15) can be rewritten 


ou 
oy 


From eqs. (14 and 16) we compute the velocity increment to a neighboring 
point (0, Ay) obtaining 


u(O, Ay) = uw + (2E sin 2Ao) Ay 
v(0, Ay) = —(E cos 2Ao) Ay (17) 


= 2K sin 2Ao (16) 


This gives enough information to determine the total velocity and the 
streamline direction at the new point, from which the procedure can be 
continued. Starting from a series of points on the known streamline, 
the streamline pattern can be developed. 

In the case of a monodisperse population in which b is known it is pos- 
sible to set up an iterative procedure for determining the principal strain 
rate directions from the optical data in the quadratic region, but pre- 
liminary attempts make it seem doubtful if it is worth the effort. For 
polydisperse systems a look at the equations** is enough to discourage 
working outside the linear range except in the case in which the shear 
component dominates. If we are dealing with a flow in which the maxi- 
mum deviation of sin 2A9 from unity can be neglected, the analysis of the 
flow is considerably simplified. If we need a knowledge of the angle of 
isocline, as required in the above method, we must limit the analysis to 
the case in which Ay ~ 45°. The application of the above method is 
readily extended to this case, obtaining the appropriate value of EH from 
the plot of An against EF for Couette flow, even if not linear. The applica- 
tion to pure liquids, where An is proportional to # and the angle of isocline 
is always coincident with the principal strain rate axis is obvious. 

Since variation in Ao has considerably less influence on An than on x, 
the method of Bogue and Prados®* using only An appears very attractive. 
In their published example of converging channel flow, they use cylindrical 
coordinates which are streamline coordinates for their case. They also 
introduce a stream function, but it does not appear to be fundamental to 
their method, and the physical picture is clearer without it. Let us 
first introduce general two-dimensional orthogonal curvilinear coordinates 
a, 8 (Fig. 8b) with the line element of the form 


ds? = h,? da? + h,? dp? (18) 
The strain rates will be 
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If we now take the curves 8 = const. to be the streamlines of the flow and 
a = const. their orthogonal trajectories we will have ug = Oand uz = V. 


(19) 


Keeping in mind that 


1 Oh, l 1 Ohg ] . 

= 5 = — (20) 
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where r is the radius of curvature of the streamline and r’ that of its 
orthogonal trajectory, we write 


(21) 
a = 4 
EAR 7ecce 


The expression for 47 becomes 


ma! (ory OV av ‘ (2 1 n) v:| is 
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A network of values of £ is available from the optical measurements. 
Starting from a wall, which we know to be a streamline, and along which 
V = 0 from the no-slip condition we can calculate (OV)/(On) from eq. 
(22). (We know r from the geometry of the wall, and must estimate 
r’.) With a value for (OV)/(On) we can calculate the velocity increment 
along the orthogonal to the streamline, and hence carry the velocity 
field on into the flow. An iterative process can be used to improve the 
estimation of r’ as the streamline pattern is developed and the orthogonal 
trajectories can be more precisely drawn. In cases in which the general 
comportment of the streamlines can be estimated in advance the calcula- 
tions will be less difficult. 


V. Solutions of Polymers 


There is a considerable body of experimental data, much of it due to 
Philippoff and his coworkers, on various polymeric systems which follow 
a stress-optical law, at least in the simple flow situations so far reported, 
of the type used so effectively in photoelastic studies. Caution must be 
used, however, in concluding that all polymer solutions have these same 
properties. A brief discussion of some of the theoretical background will 
help keep perspective on the use of SBR for flow visualization studies in 


polymeric systems. 
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Philippoff* has given a rather extensive treatment of the interrelation- 
ship to be expected between elastic stresses and the optical properties of 
viscoelastic liquids if it is assumed that there is a simple correlation be- 
tween the elastic deformation tensor and the optical tensor. Lodge*’ 
appears to be the first to develop a rational theoretical model explaining 
such relationships. He pointed out that, in certain flowing liquids one 
might expect the same mechanism, such as network deformation, to give 
rise both to birefringence and stress. If these two phenomena arise 
from the same source, it is reasonable to suppose that the optic and stress 
ellipsoids will be coaxial, and the differences of the principal axes of either 
ellipsoid will be proportional to the differences of the corresponding axes 
of the other ellipsoid. He developed a model compatible with this hy- 
pothesis based on the method of Green and Tobolsky in their theory of 
rubberlike elasticity. The molecular interactions between the molecules 
of the polymer are assumed to be sufficiently localized and long-lived to be 
treated as junctions in a deforming network, and he considers that the 
stress, apart from an additive isotropic component, arises entirely from 
entropy changes consequent on the deformation of the network. Some 
parts of the network will be older than others and will be more deformed, 
but it is assumed that a network element of a given age gives its own 
contribution to the total stress, independent of the rest of the network, 
of an amount given by the statistical theory of rubberlike elasticity. 
The total stress is the sum of the contributions. The birefringence is 
assumed to be due to orientation by network deformation of optically 
anisotropic chain links. Hydrodynamic flow forces are neglected, as are 
contributions from orientation of free chains by solvent flow and from 
double refraction of entire macromolecules. The flow gradient is also 
taken to be so small as not to affect the mean lifetime of junction points. 
On this model he concludes that, for laminar flow in the x; direction, with 
a gradient I in the x2 direction, the angle of isocline will be given by 

2821 


tan 2x = — 
811 — See 


i.e., the isoclinie axes coincide with the principal stress axes; and that 


An sin 2x = (=) ay( =) amet (24) 
nt o 3n k1 


where # is the mean index of refraction; a, a are principal polarizabilities 
of a chain; 7 is the viscosity of solution; and 7 is the viseosity of solvent. 
With this model he also finds 


821 
“= * const. (25) 


This lack of dependence of 7 on T is due to his assumption of a quasi- 
tatic process: it is assumed that a chain, on leaving the network, is ables 
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to retract fully before joining the network again. The concentration and 
shear rate range are, therefore, limited. In particular, the solution must 
be sufficiently concentrated to possess network structure. The upper 
limit of concentration is also influenced by the assumption that the mean 
distance between successive chains is large enough for a Gaussian ap- 
proximation to be valid. 

Philippoff and his coworkers*** have studied the relation of bire- 
fringence and shear stress for solutions of polyisobutylene (PIB) in dec- 
alin from concentrations of 0.1-9% as well as 100% PIB over a shear 
range from 10—' to 10‘ sec.—'. His results for the various solutions are 
shown in Figure 9a. Using viscosity data obtained with various types of 
viscometers, depending on the range, these are replotted in Figure 9b 
in which An, the amount of birefringence, is plotted against half the dif- 
ference of principal stresses 7/sin 2x. This plot uses the whole shearing 
stress (i.e., the viscosity of the solvent has not been subtracted out as 
suggested by Lodge). Subtracting out the solvent viscosity would pull 
the points for the lowest concentrations to the left just about to the limit 
of scatter of the data, and would have less and less influence at the higher 
concentration. The linearity of the data as plotted may be an argument 
for using the viscosity of the solution, as suggested by Philippoff, or, if the 
corrected data must be used, may represent a departure of the behavior 
of the solution from the Lodge theory. There is some real doubt if the 
Lodge theory should be expected to be valid at the lower concentrations, 
since Janeschitz-Kriegl*® has shown that form birefringence becomes 
important in such cases: but if this is true, the stress-optical law will 
not be valid, either. 

Philippoff has also made direct normal stress measurements using a modi- 
fied Weissenberg rheogoniometer. These data are plotted in Figure 10. 
The dotted lines represent the same data as deduced from the birefringence 
measurements, correlating with the total shear stress and not correcting 
for the viscosity of the solvent. We see that the fit is pretty good, but far 
from perfect, especially considering the logarithmic scale. Unfortunately 
there are no direct mechanical measurements at the lower concentrations. 

From these results it appears that, for PIB in decalin, the stress-optical 
law is well enough demonstrated to be usefully applicable over a wide 
range of stresses and a wide range of concentrations. For the higher 
concentrations, it is not particularly important whether the SBR depends 
on (» — m) or on 7. The optical measurements are, in general, more 
sensitive than mechanical measurements for normal stresses provided we 
van be sure that the deviation of the angle of isocline from the principal 
strain rate direction is due to stress effects and not to orientation effects. 
The fact that the data for 100% PIB also fell in line with the solution 
data raises further theoretical questions since Lodge’s theory definitely 
should not hold for this ease. 

Concentrated polymer solutions have been ‘considered by Janeschitz- 
Krieg] in another quasi-static model. He concludes that the amount of 
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Fig. 10. Birefringence as a function of difference in principal stress for polyisobutylene in 
Decalin. 


birefringence depends on I'(n — mo) provided the velocity gradient is suf- 
ficiently small. He also calculates the amount of form birefringence for 
isolated polymer molecules in dilute solutions when the index of refraction 
of the polymer and that of the solvent are different, and finds a distinct 
and measurable effect. Cerf‘! has looked at the case of dilute solutions 
of deformable molecules and has shown that for very low velocity gradients 
they should behave like rigid particles, for very high gradients like de- 
formable systems, with all intermediates in between. This type of be- 
havior has been experimentally observed by Leray.*? 

We can summarize the theoretical findings by saying that, if we can 
describe the mechanism producing SBR by a quasi-static statistical process 
in which form birefringence is small and overall orientation of asymmetric 
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macromolecules is not important, the stress-optical law will be valid. 
If, on the other hand, the form of the flow field is varying too rapidly from 
point to point for equilibrium to be approximated, we should not expect 
to have a simple stress optical law. The question naturally arises: how 
fast is too fast? This will depend on the relaxation time spectrum of the 
elements giving rise to the SBR. Just because we appear to get a simple 
stress optical law for a substantially stagnant flow such as Couette flow does 
not mean a priori, that the stress optical law will still be valid in a more 
general situation: the rate of change of hydrodynamic ‘history’? must 
be slow enough that the molecular arrangement reflects the local conditions. 
If we must take the prior history into account we have a rather cumber- 
some tool, at best. For the systems, such as polyisobutylene, for which 
a stress-optical law is valid in Couette and capillary flow there still seems 
to be too little evidence to judge how rapidly the optical properties will 
follow a change in flow conditions, either for transient changes or com- 
plicated steady flow patterns. 

Assuming we have a fluid which strictly follows the stress-optical law, 
what of a quantitative nature can we deduce from the optical: measure- 
ments? First we must confine our flow studies to two-dimensional situa- 
tions. If we measure both the angle of isocline and the amount of bire- 
fringence we can deduce the direction of the two axes of the deviatoric 
stress ellipsoid lying in the plane of the flow, and the difference of the prin- 
cipal stresses, and no more. For two-dimensional incompressible flow we 
have two equations of motion, the equation of continuity and the con- 
stitutive equation with which to work. We have six unknown quantities 
in the plane: sy, S12, 822, u, v, and p. We know sy and sy — 8, so pre- 
sumably we have enough information to reconstruct the flow. Practi- 
cally this is not simple. We probably do not know enough about the 
system to use the constitutive equation effectively, and the equations of 
motion are, in general, nonlinear. The reconstruction of the deviatoric 
stress field can be carried along lines equivalent to those used for flow 
analysis of colloidal solutions given in Section IV, with the simplification 
that the optic axes are coaxial with the principal stress axes. Bogue and 
Peebles® suggest using the Airy stress function as a means of construct- 
ing the stress pattern from the SBR data, similar to the method of re- 
constructing the flow pattern in Newtonian systems from the stream func- 
tion. This will require a measurement of the stress distribution (including 
the normal stresses) over the boundaries. The stress distribution will, 
however, not give the flow pattern unless we know the relationship between 
the stress tensor and the strain-rate tensor for the polymer being studied. 

One of the most ingenious and important applications of the stress- 
optical law has been to check the Weissenberg hypothesis concerning the 
equality of se. and s;; in shear flow. If the magnitude of the SBR is, 
indeed, proportional to the difference of the principal stresses, then if 
observations are made along the streamlines in shear flow no birefringence 
should be detected. Philippoff‘* made observations along the axis of a 
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capillary in a variety of solutions of polyisobutylene, rubber, polysiloxane 
and nitrocellulose, finding no birefringence, hence deducing that the Weis- 
senberg hypothesis is valid for these systems. Systems of rigid particles, 
such as vanadium pentoxide and tobacco mosaic virus showed marked 
birefringence, as would be expected for orientable systems. It would be 
interesting to extend this type of experiment to the type of system dis- 
cussed by Cerf‘! in which orientation effects should predominate at low 
shear and deformation effects at high shear. 


VI. Summary and Conclusions 


1. Qualitative Studies 


The use of streaming birefringence for qualitative observation of the 
character of flow separation, onset of turbulence, existence of secondary 
flows, etc., has already proved extremely effective. Great care must be 
exercised, however, that the rheological properties of the birefringent 
liquid used are fully understood, and, indeed are those of the system of 
interest. Even small departures from Newtonian behavior may introduce 
marked qualitative differences in the nature of the flow. SBR should be 
particularly valuable in qualitative studies of the flow of viscoelastic 
systems to detect the existence or absence of secondary flows. 


2. Quantitative Studies 


(a) Pure Newtonian Liquids. Pure liquids with Newtonian viscosity 
behavior can probably be used for quantitative analysis of complex two- 
dimensional flows using SBR if the measuring techniques are sufficiently 
sensitive. The case in which the fluid elements change hydrodynamic en- 
vironment with great rapidity has not been checked experimentally al- 
though the short relaxation time for the small molecules of these substances 
will probably insure that the optical properties are reflecting the local flow 
situation rather than the past history. 

(b) Colloidal Solutions. Flow analysis from the SBR data for colloidal 
solutions can be made with assurance only if all portions of the flow are 
within the range in which both the retardation An and the variation of the 
angle of isocline x are linear with velocity gradient in Couette flow, the 
change in birefringence with concentration is linear, and the relaxation 
time of the colloidal particles is short enough that one can be sure that the 
local optical effect is due to the local flow field and does not depend on the 
upstream history of the particles. This is, indeed, a very stringent 
limitation. For many engineering purposes one can work considerably 
outside this range to obtain “adequately” quantitative data but one must 
recognize that the discrepancies with theoretically predicted results are 
probably due to systematic errors and not to random errors. 

(c) Polymeric Solutions. Certain polymeric solutions, particularly poly- 
isobutylene in decalin have been sufficiently studied to assure that a stress 
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optical law is valid and the Weissenberg hypothesis is also valid. Fora study 
of the stress patterns in flow of such viscoelastic systems, SBR should be 
most useful. Care must be taken that time-history effects are not impor- 
tant (unless it is these one wants to study). If a new polymeric system is 
to be used, it would be essential to ascertain that the stress-optical law is 
valid; that the Weissenberg hypothesis is valid; and that the relaxation 
time is short enough that the optical properties correctly reflect the local 
stress situation. 

The construction of the flow pattern for these systems will require 
measurements other than with SBR, or the knowledge of the constitutive 
equation of the system. 


3. Instrumental 


Photographic techniques are the easiest to use if the amount of re- 
tardation is sufficient to give adequate resolution. It seems probable 
that even boundary layer studies might be possible by using a light beam 
which has virtually no converging or diverging rays, as can be obtained 
with a laser. This should be particularly valuable with polymer solutions, 
where large effects are often available. It seems questionable if any 
colloidal system will show enough birefringence for the photographic 
method to be applicable for quantitative investigation, or even for qualita- 
tive investigations unless non-Newtonian effects are being studied. 

Point-by-point measuring techniques can be used for spatial resolution 
of the order of 0.1 mm. and can be made extremely sensitive and relatively 
easy to use with photoelectric detection. 
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Résumé 


Cet article rassemble les données théoriques et expérimentales généralament valables 
pour établir des critéres permettant 4 l’expérimentateur qui le désire d’utiliser la biré- 
fringence d’écoulement comme un instrument pour |’étude de I’écoulement fluide. On 
y décrit et on discute les techniques expérimentales et les possibilités pour rendre visible 
l’écoulement qualitatif et quantitatif par la biréfringence d’écoulement dans les liquides 
Newtoniens puis les suspensions colloidales et des solutions polymériques diluées ou 
concentrées. On signale des méthodes de reconstruction de la vitesse et des champs de 
tension a partir des données de biréfringence. 


Zusammenfassung 


In der vorliegenden Mitteilung werden die gegenwirtig vorhandenen theoretischen 
und experimentellen Unterlagen zur Aufstellung von Kriterien, die dem Experimentator 
bei der Beniitzung der Str6mungsdoppelbrechung als Werkzeug zur Untersuchung des 
Fliissigkeitsfliessens behilflich sein kénnen, besprochen. Ein Uberblick tiber die Ver- 
suchmethodik wird gegeben und die Méglichkeit fiir eine qualitative und quantitative 
Darstellung des Fliessens durch Strémungsdoppelbrechung in rein Newton’schen 
Fliissigkeiten, killoiden Suspensionen sowie verdiinnten und konzentrierten Polymerlé- 
sungen wird diskutiert. Methoden zur Rekonstruktion des Geschwindigkeits- und 
Spannungsfeldes aus Doppelbrechungsdaten werden angegeben. 
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Synopsis 


Aqueous suspensions of Milling Yellow dye have been found to be useful for engineering 
flow studies by employing established birefringent techniques. In the present work, the 
basic birefringent and rheologic properties of Milling Yellow suspensions were measured 
for a range of deformation rates important to the engineering applications. A concentric- 
cylinder polariscope was used to measure the amount of birefringence and the angle of 
extinction for ten Milling Yellow suspensions at 25°C. varying in concentration from 
1.25 to 1.46% dye by weight. The optical properties exhibited by these suspensions 
are similar to those of a number of colloidal suspensions as reported by past investigators. 
Rheologic data for Milling Yellow suspensions in the same concentration range were 
measured at 25°C. for deformation rates of 50-10‘ sec.~! by means of a capillary vis- 
cometer. Over the complete range of deformation rates the rheologic data are accu- 
rately described with the Powell-Eyring relationship for pseudo-plastic materials. 


Introduction 


The engineering application of optical techniques based on flow bire- 
fringence for determination of velocity profiles and stress distribution in 
laminar flow, detection of regions of flow separation and boundary layer 
reattachment, and onset of turbulence is well-known and employed widely. 
In 1953, some of the possibilities for visualization of flow phenomena with 
aqueous suspensions of Milling Yellow dye were discussed," and since then 
it has been shown that quantitative measurements of velocity and stress 
distributions can be made for two-dimensional laminar flows.?'* Other 
work with Milling Yellow dye suspensions includes the dynamic bire- 
fringence studies of Thurston, '4 the work on laminar-turbulent transition by 
Gilbrech,® and a large number of flow studies not reported in the literature. 

It is the purpose of this paper to present experimental data on the basic 
birefringent and rheologic properties of Milling Yellow suspensions. Such 
data were obtained at 25°C. for a range of dye concentrations of interest for 
flow studies. Birefringence and extinction angle measurements were made 
with a concentric cylinder polariscope using well-established techniques.* 
The rheologic measurements were made with a capillary viscometer. 
Viscosity data were obtained from the pseudo-shear diagrams by standard 
methods.” 
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Optical Theory 


The use of plane polarized light to measure the amount of birefringence 
produced by liquid motion is well-known, and the quantitative features of 
the phenomenon have been described by many writers.*® In summary, if a 
plane-polarized light wave moves through a birefringent liquid and then 
through a second polarizer, the vibration of the transmitted light wave can 


be represented by 


; . 10 
u = asin 2¢ sin x cos (ct — x — 6/2) (1) 


where u is the displacement of the vibration from the axis at time t; a is 
the amplitude of the vibration; A is the wavelength of the light vibration; 
« is the arbitrary phase retardation indicating the displacement is not 
necessarily zero at t = 0; 6 is the relative retardation due to passage 
through the birefringent liquid; and ¢ is the angle between optic axis and 
plane of vibration of the incident wave. 

Since the intensity of transmitted light is proportional to the square of 
its amplitude of vibration, eq. (1) indicates that zero intensity results when 


@ = N(x/2) andé = NA (2) 


where N is the fringe order (= 0, 1, 2, 3, .. .). 
These are the well-known conditions for the extinction angle and the 


isochromatic fringes, respectively. 
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Isocline 
Fig. 1. Top view of annulus of the concentric cylinder polariscope showing the cross 
of isocline and the optical orientation of the birefringent liquid. Parallel lines indicate 
plane of vibration of light leaving the polarizer (direction OX in Fig. 1). 
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The extinction angle X is defined by Edsall* as the angle made by the 
optic axis of the birefringent liquid with the streamlines. In the concentric 
cylinder polariscope used in this work, four dark areas, forming the well- 
known cross-of-isocline, appear in the annulus of the flow cell at the points 
where ¢ = 0, 7/2, x, and 3(m/2), as indicated by eq. (2). This is shown 
pictorially in Figure 1. The small angle made by the arms of this cross 
with the axes of polarization OX and OY is the extinction angle, x. 

The isochromatic fringes which are indicative of the amount of bire- 
fringence appear as dark bands concentric with the axis of the flow cell. 
These bands move outward across the annulus and the higher order fringes 
move into the annular field in succession as the speed of the outer cylinder 
is increased. The amount of birefringence can be measured by counting 
the isochromatic fringes as the speed of the rotating cylinder is increased 
from zero, thereby determining N. (N = 0 corresponds to the zero retarda- 
tion occurring when the liquid is at rest.) The amount of birefringence ex- 
pressed as the difference in the indices of refraction for the extraordinary 
and the ordinary rays passing through the test liquid contained in the flow 
cell can be calculated from 

Nx 
|ne — mol = (3) 


Y 
A 


where |n, — mo| is the amount of birefringence and S is the length of light 
path in the liquid. 

{quation (3) makes use of the condition for the existence of the isochro- 
matic fringes as given by eq. (2) and the calculation of the relative retarda- 
tion of the two light waves as determined by the length of light path and 


light speed for the two rays. 





Experimental 


Test liquids were prepared from commercial Milling Yellow NGS dye 
(National Aniline and Dye Co.) by mixing the powder with distilled water 
to give about 1% dye by weight. The mixture was then heated to the boil- 
ing point to dissolve the dye, and evaporated to the desired concentration. 
For the test liquids investigated the dye concentration was determined by 
carefully evaporating a weighed sample of the liquid to dryness at 100°C., 
and then weighing the residue. Determinations of the dye concentration 
by this means could be reproduced within +0.002% dye by weight. 

Birefringent and rheological properties of Milling Yellow dye suspensions 
in the concentration range of 1.25-1.45% by weight were determined at 
25.0°C. Following is a brief description of the concentric cylinder polari- 
scope, capillary viscometer, and experimental procedures employed in the 
measurements. 


Concentric Cylinder Polariscope 


The measurements of extinction angle and amount of birefringence for the 
Milling Yellow dye suspensions were carried out in a concentric cylinder 
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Fig. 2. Schematic diagram showing the components of the concentric cylinder polariscope. 


polariscope manufactured by the Rao Instrument Company, New York. 
An instrument almost identical in every respect with the one used in the 
present work was described by Edsall,* and their article may be consulted 
for any details not given here. The assembly of major components of this 
unit is indicated schematically in Figure 2. All components except the 
monochromatic filter are attached to an optical bench for rigid, precise 
alignment. 

The light source is a tungsten arc lamp (Sylvania type C-25) which gives 
an intense light, emanating from a small area, thereby approximating a 
point source. It is mounted at the focal point of a collimating lens, which 
focuses parallel rays of light upward through the bottom of the flow cell. 
A monochromatic light filter (green, wavelength 5100 A.) is placed above 
the lens for measurements of the birefringence, but is removed for extinction 
angle determinations. The polarizer is located below the flow cell in a 
friction mounting, which allows it to be adjusted and removed for cleaning. 
The analyzer is mounted above, on a support which may be removed com- 
pletely from the optical bench to allow the flow cell to be filled con- 
veniently. The analyzer support assembly also contains a quarter wave 
plate which may be positioned below the analyzer for use as a compensator. 
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Fig. 3. Cross sectional assembly drawing of concentric cylinder flow cell. 


Its utilization would enable one to measure amounts of double refraction 
corresponding to a retardation of less than one wavelength of light. How- 
ever, the Milling Yellow suspensions employed in the present investigation 
produced relatively large retardation at low rates of shear, and the com- 
pensator was not needed. Not indicated in Figure 2 is a set of movable 90° 
cross hairs, mounted between the polarizer and the bottom of the flow cell. 
These were used to measure the extinction angle. 

A drawing of the flow cell is shown in Figure 3. The annular flow region 
is formed by the space between the fixed inner cylinder and the rotating 
outer cylinder. The inner cylinder is mounted from above by a supporting 
plate, which holds it free from the outer cylinder and assures that its 
position is concentric to that of the outher cylinder. The outer cylinder 
assembly is made up of two elements; a rotor, mounted on sealed ball 
bearings and equipped with a drive pulley, and a cup which fits in the rotor 
and forms the outer surface of the flow cell. The bottom of the cup is of 
glass, and a glass top window prevents the formation of a meniscus in the 
liquid, which would act as a lens and distort the parallel light rays. Except 
for the glass top and bottom windows, all cylinder parts are of stainless 
steel. Each glass piece is mounted between two neoprene washers, and is 
held in place by a threaded fitting. The annular dimensions can be varied 
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by interchanging the three inner cylinders provided. This permits measure- 
ments to be made over a wide range of shear rates, with the same range of 
instrument rotation speeds. The radial dimensions of the cylinders are as 
follows: outer cylinder, 0.500 in.; large inner cylinder, 0.480 in.; medium 
inner cylinder, 0.460 in.; small inner cylinder, 0.349 in. 

The liquid depth in the cell is determined by the distance from glass 
bottom to top glass ring of the outer cylinder cup. For the cylinders 
furnished by the manufacturer, this distance was 1.768 in. However, in 
order to increase the clarity of the interference bands in the liquid, four 
new cylinders were fabricated, one outer, and three inner, with radical 
dimensions identical to those of the original cylinders but shorter in 
length, so as to provide a liquid depth of 0.570 in. This modification was 
made to improve fringe definition at high orders. For fringe orders above 
nine to ten, the observed optical phenomena in the polariscope become in- 
distinct. This results from overlapping of light waves of slightly different 
wavelengths, because of incomplete filtering action by the monochromatic 
light filter. 

The outer cylinder of the flow cell is connected by a pulley and an elastic 
belt to the drive motor. The drive unit itself consists of a variable-speed, 
direct current operated motor and a regulated power supply. The motor is 
provided with six pulleys, three operating on a low speed shaft driven by a 
gear train, and three directly on the armature shaft. This arrangement pro- 
vides a continuous speed range for the outer cylinder of the flow cell, from 
approximately 0.8 to 4000 rpm. This range was more than adequate for all 
the measurements made in this investigation. 

The speed of rotation of the cell was measured by counting revolutions 
of the outer cylinder for timed intervals, up to about 120 rpm. At higher 
speeds, a flashing stroboscope (Type BL-631 Strobotac, manufactured by 
the General Radio Corporation) was used. 


Birefringence and Extinction Angle Measurements 


The techniques used in the measurement of the extinction angle and the 
birefringence, while simply in principle, required considerable practice and 
refinement before reproducible results were obtained. Two factors, in 
particular, worked to disrupt the accuracy of the results: (a) temperature 
variations and (6) evaporation of water from the solution in the flow cell. 
Since the birefringent properties of Milling Yellow suspensions are highly 
dependent on temperature and dye concentration, both of these sources of 
error were important. A temperature rise of 0.2-0.3°C. was sufficient to 
lower noticeably the amount of birefringence at a given rate of shear. There- 
fore, it was necessary to empty the cell and allow it to cool at frequent 
intervals. The effect was most marked in tests on the more viscous solu- 
tions, and even with the dilute solutions at high rates of shear. The effects 
of evaporation were noticeable if the fluid was left in the cell for longer 
than fifteen or twenty minutes. This was manifested in an increased 
amount of birefringence (due to increased dye concentration) for a given rate 
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of shear. By making measurements rapidly and refilling the cell at frequent 
intervals, it was possible to keep the effects of evaporation to a minimum. 

To prepare the concentric cylinder polariscope for measurements, first 
the optical system was adjusted. The collimating lens was set to give 
parallel rays of light by placing the inner and outer cylinders in position, 
and projecting the image of the annulus on a piece of white cardboard, held 
about a foot above the cell, in a darkened room. The lens was moved until 
a clear sharp image of the annulus, without halos, appeared on the card- 
board. The plane of polarization of the polarizer was oriented with the 
cross hair scale by use of a 3° bifield polarizer. A 3° bifield polarizer is a 
plate of polarizing material really consisting of two polarizers whose planes 
of polarization are 3° apart. A line scribed on the plate bisects these 
planes. To adjust the instrument polarizer, the bifield was set so that this 
line was perpendicular to the desired plane of polarization (aligned with a 
cross-hair at the 90° position). Then the instrument polarizer was turned 
until both halves of the bifield appeared equally dark. A very accurate 
adjustment could be made in this manner. The analyzer support assembly 
was mounted on the optical bench, and the analyzer turned until it produced 
maximum extinction of the light coming from the polarizer. These adjust- 
ments for the polarizer and analyzer were checked at intervals, but rarely 
needed correction. 

When the optical system had been adjusted, the inner and outer cylinders 
were assembled as shown in Figure 3, and were slowly filled with a long- 
needled syringe. Care was taken to prevent the entrapment of air bubbles 
in the cell. Once the cell had been filled, the optical measurements were 
begun. 

The measurements of the amount of birefringence were carried out as 
follows. The motor was started and its speed slowly increased until a dark 
band was observed to move outward across the annulus. When the band 
was at the center of the annulus (this could be estimated accurately by 
eye, since the annular gap was quite narrow) the speed of rotation of the 
outer cylinder was measured. This was done by timing a number of revolu- 
tions with a stopwatch at low speeds (less than 120 rpm) or with the 
Strobotac at higher speeds. The speed of the cylinder was increased and 
recorded as each band moved into the center of the annulus. At intervals, 
the cylinder was slowed down so that previously observed bands appeared, 
and second speed measurements were made. This served as a check as to 
whether or not the liquid was being heated up enough to change its bire- 
fringent properties. If a point could not be reproduced, the motor was 
stopped; the cylinder was cleaned and allowed to cool, then refilled with 
fresh liquid and the irreproducible portion of the run repeated. The 
cylinder was never run for more than fifteen or twenty minutes at a time 
without refilling, in order to prevent the evaporation effects from becoming 
noticeable. 

Much trial and error indicated that better results could be obtained by 
making separate runs to measure the extinction angle, rather than trying to 
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measure it along with the amount of birefringence. Attempts to measure 
the two simultaneously, slowed down the run to such an extent that only a 
few measurements could be made before it became necessary to refill the 
cell to prevent evaporation effects. The extinction angle measurements 
were made by aligning one of the instrument cross-hairs to coincide with an 
arm of the cross of isocline, as viewed through the analyzer. Four such 
measurements were made at a given speed of rotation, one for each arm of 
the cross, and the arithmetic average of the readings was taken. This pro- 
vided somewhat greater accuracy ii the extinction angle values, which were 
difficult to determine, even at low speeds of rotation, to much closer than 
+0.5°. At high speeds, variation in rate of shear across the annulus, due 
to non-Newtonian flow behavior, caused the cross to spread out and be- 
come quite diffuse. This reduced the accuracy of the extinction angle 
readings to about +2.0°. Since the extinction angle values were less ac- 
curate and did not seem to be as sensitive to changes in temperature and 
concentration as for the amount of birefringence, it was not necessary to 
empty and refill the cell as frequently during the extinction angle measure- 
ments as it was during measurements of the birefringence. 

The amount of birefringence was calculated from the observed fringe 


order means of eq. (3). 
The shear rate at the center of the annular space of the flow cell was 


calculated from 


Gf = eee (4) 
(7) + r2)?(r2? — 11”) 


where G@ is the shear rate; is the speed of rotation of outer cylinder; 1; is 
the radius of inner cylinder; and rz is the radius of outer cylinder. 


Capillary Viscometer and Procedure 


The capillary viscometer used for rheologic measurements is a simple 
modification of the Ostwald-type. The important features are illustrated 
in Figure 4. The capillary section is approximately 0.04 em. i.d. by 20 cm. 
long, and the sample bulbs are nominally 5 cc. internal volume. The di- 
mensions of the capillary and the volumes of the sample bulbs were deter- 
mined from very careful physical measurements and viscometric measure- 
ments with standard oils. In general it was found the viscosity of standard 
oils determined from the dimensions of the capillary agreed with the 
National Bureau of Standards’ values within 0.1%. 

Figure 5 shows schematically the viscometer with the necessary auxiliary 
equipment for pressure drop measurement and temperature control. Runs 
were made by timing with a stopwatch the interval for the volume of test 
liquid in the sample bulb to flow through the capillary. Data for a wide 
range of shear rates were obtained by using a range of drive-fluid pressures 


from a few millimeters to about 100 em. of mercury. 
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Fig. 4. Assembled capillary and flow control unit 
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Fig. 5. Schematic diagram of the capillary viscometer facility. 
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Viscosity values were calculated from the apparent values by means of 
the relationship 
] ] Tw d(1 ‘Ma) 
hee nee (5) 
a" 4 dty 
where yu is the viscosity at shear stress, 7»; ua is the apparent viscosity; and 
Yw is the wall shear stress. 

The apparent viscosity was calculated from the Hagen-Poiseville equation 
for steady laminar flow in a capillary. Values for the derivative in eq. (5) 
were obtained by numerical differentiation of the experimental relationship 
between apparent viscosity and wall shear stress. The shear rate corre- 
sponding to the viscosity was calculated from the equation, 


G=- (6) 
a 
Results 


Experimental results for the amount of birefringence and the extinction 
angle for nine Milling Yellow suspensions are given in Figures 6 through 10. 
The amount of birefringence is shown plotted as a function of shear rate in 
Figures 6 arid 7. For the liquids tested it can be noted that at low 
shear rates the birefringence is proportional to shear rate as suggested by 
the theoretical work of Boeder,'! Peterlin and Stuart,!! and Wayland." 
In the case of increasing shear rates the relation between Ine _ no| and shear 
rate becomes nonlinear, and for very high shear rates the data suggest a 
second range of linearity. The results clearly reveal the marked de- 
pendence on dye concentration in the range of 1.2-1.5% dye by weight. 
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Note: Curve designations are percentage 
by weight of dye in solution 
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Fig. 6. Birefringence as a function of rate of shear for liquids tested. 
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Note: Curve designations are percentage 
by weight of dye 
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lig. 7. Birefringence as a function of rate of shear for dilute solutions (high speed range). 
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Fig. 8. Extinction angle as a function of rate of shear for four representative solutions 
in a moderate shear range. 


Below about 1.2% dye the Milling Yellow suspensions do not exhibit 
birefringence, and at the higher concentrations where intense birefringent 
behavior is found, the dye suspensions are gels and show birefringence even 
in the absence of flow. 

Values of the extinction angle are plotted in Figures 8, 9, and 10, a func- 
tion of shear rate. In all instances, the extinction angle is 45° at zero shear 





48 F. N. PEEBLES, J. W. PRADOS, AND E. H. HONEYCUTT, JR. 


rate and then decreases to saturation values near 20° at very high shear 
rates. The shear rate range for saturation extinction angles decreases with 
increasing Milling Yellow concentration. 

The birefringent properties of Milling Yellow suspensions, as reflected 
by the high values of birefringence (10-4 compared to 10~'? for liquids 
showing molecular birefringence) and the shear rate dependence, suggest 
that these liquids can be considered as suspensions of rigid rods. If one 
accepts this model then the experimental values of extinction angle can be 
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Fig. 9. Extinction angle of three concentrated solutions at low rates of shear. 
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Fig. 10. Extinction angle of three dilute solutions at high rates of shear. 


used with theoretical results on particle orientation! to show that the effec- 
tive particle diameter is in the Milling Yellow suspensions about 0.5 u. 
Details of these calculations are given.!” 

In general, Milling Yellow suspensions exhibit non-Newtonian flow 
characteristics. Figures 11 and 12 show typical rheologic curves plotted as 
shear diagrams. The data for liquids in the range of concentrations from 
about 1.2 to 1.5% dye by weight show linearity at low and high shear rates, 
and as in the case of the optical data show a marked dependence on dye 
concentration. Viscosity values range from near 1.5 ep at shear rates ap- 
proaching 10* reciprocal sec.-400 ep at very low shear rates, depending on 


the dye concentration, 





—— 





a 


SE Ee ee 





MILLING YELLOW SUSPENSIONS 19 


As cited in the previous section, the shear rates were obtained by calcu- 
lation from the apparent viscosity values. The basis of this established 
method of calculation is that the shear stress is a unique function of shear 
rate. One test of the validity of this assumption is that the apparent vis- 
cosity of a given liquid is a function of wall shear stress and independent of 
capillary diameter. Figure 13 shows data on the apparent viscosity of a 
Milling Yellow suspension containing 1.442% dye as measured in three 
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Fig. 11. Wall shear stress as a function of shear rate for the solutions tested. 
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Fig. 12. Wall shear stress as a function of shear rate for the solutions tested. 
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vapillaries. The narrow band of dispersion of data points about the mean 
curve drawn through the data suggest the validity of the assumption re- 
garding uniqueness of the shear stress—shear rate relationship. 
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Fig. 13. Apparent viscosity as a function of wall shear stress for the 1.442% solution 
using three different capillaries. 
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Fig. 14. Viscosity data for 1.32% Milling Yellow suspension. 
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Fig. 15. Powell-Eyring parameters for Milling Yellow suspensions. 


The viscosity data for the Milling Yellow suspensions are well repre- 
sented by the Powell-Eyring equation for pseudo-plastic materials, ex- 
pressed as 


a A(uo — wo) sinh—!(G/A) 7 
p=, + Sere 7) 
where jo, He are Viscosity at zero and infinite shear rates; A is a constant; 
and G is the shear rate (= —(du/dr) for tube flow). 

The constants of eq. (7) were obtained from experimental data for a given 
test liquid by the familiar method-of-averages,’ i.e., the constants were 
calculated from the requirement that 


M etic inh-! (G,/2 
2 4 ~ E + Ae eu = 0 (8) 
+=0 e G; 
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Figure 14 shows a typical plot of viscosity data and comparison to the 
Powell-Eyring equation. In general, the agreement between experimental 
points and the correlation equation was within +1%, with an average 
deviation of about +0.5%. Figure 15 shows values of the Powell-Eyring 
parameters for a range of Milling Yellow concentrations. 

Use of the Powell-Eyring equation to describe the Milling Yellow vis- 
cosity data does not imply that other viscosity models are not applicable. 
Work is under way to examine other representations that will allow a unified 
treatment of the rheologic and birefringence data for Milling Yellow 


suspensions. 


Conclusions 


1. Experimental measurements of the amount of birefringence and ex- 
tinction angle for Milling Yellow dye suspensions over a wide range of shear 
rates show a behavior similar to that observed with other colloidal suspen- 
sions. The shear rate dependence of these optical properties suggest that 
the suspensions may be considered as suspensions of rigid, ellipsoidal 
particles. 

2. Viscosity data for Milling Yellow suspensions reveal non-Newtonian 
behavior in general, with Newtonian ranges at the low and high shear rate 
extremes. Accurate description of the rheologic data for the complete 
range of shear rates can be made with the Powell-Eyring equation. 


The authors wish to acknowledge indebtedness to the late Robert Marsh Boarts, late 
Professor and Head of the Chemical Engineering Department at the University of 
Tennessee, who suggested to the senior author that flow birefringence techniques should 
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from the Office of Naval Research, Mechanics Branch. 
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Résumé 


Des suspensions aqueuses de colorant Milling Yellow ont été trouvées utiles pour des 
études d’écoulement en employant des techniques bien établies de biréfringence. Dans 
le présent travail, les propriétés de birefringence fondamentales et les propriétés rhéo- 
logiques de suspensions de Milling Yellow ont été mesurées dans un domaine de vitesses 
de déformation importante pour les applications industrielles. Un polariscope 4 cylindre 
concentrique a été utilisé pour mesurer le degré de biréfringence et l’angle d’extinction 
pour dix suspensions de Milling Yellow 4 25°C, en variant le concentration de 1.25 a 
1.46% de colorant en poids. Les propriétés optiques de ces suspensions sont semblables 
a celles d’un nombre de suspensions colloidales comme cela a été raporté par des études 
précédentes. Les données rhéologiques pour des suspensions de Milling Yellow dans le 
méme domaine de concentration ont été mesurées 4 25°C pour des vitesses de déforma- 
tion de 50-10‘ sec.~! au moyen d’un viscométre capillaire. Dans le domaine complet 
des vitesses de déformation, les données rhéologiques sont décrites parfaitement au 
moyen de la relation de Powell-Eyring pour des matériaux pseudoplastiques. 


Zusammenfassung 


Wiassrige Suspension en des Farbstoffes Milling Yellow erwiesen sich fiir technische 
Fliessuntersuchungen mit bekannten Doppelbrechungsmethoden als brauchbar. In 
der vorliegenden Arbeit wurden die grundlegenden Doppelbrechungs- und Fliesseigen- 
schaften von Milling-Yellow-Suspension in dem fiir die technische Anwendung wichtigen 
Bereich an Deformationsgeschwindigkeiten emessen. Zur Messung des Betrages der 
Doppelbrechung und des Ausliéschungswinkels an zehn Milling-Yellow-Suspension 
bei 25°C im Konzentrationsbereich von 1,25 bis 1,46 Gewichtsprozent Farbstoff wurde 
ein Konzentrisch-Zylinder-Polariskop verwendet. Die optischen Eigenschaften dieser 
Suspensionen weisen eine Ahnlichkeit mit den von friiheren Autoren fiir eine Anzahl von 
kolloiden Suspensionen mitgeteilten auf. Rheologische Daten wurden fiir Milling- 
Yellow-Suspensionen bei 25°C und Deformationsgeschwindigkeiten von 50-10 sec™ in 
einem Kapillarviskosimeter bestimmt. Im ganzen Bereich an Deformationsgeschwindig- 
keiten kénnen die rheologischen Daten durch die Powell-Eyring-Beziehung fiir pseudo- 
plastische Stoffe genau wiedergegeben werden. 
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Particle Size Distribution in Dilute Liquid/ 
Liquid Dispersions by Light Scattering 


K. EK. LINDSEY, University of Massachusetts, D. C. CHAPPELEAR, 

Monsanto Chemical Company, D. M. SULLIVAN, Shawinigan Resins 

Corporation, V. A. AUGSTKALNS,* University of Massachusetts, 
Amherst, Massachusetts 


Summary 


1. Equal energy per unit mass (or volume) proved to be a valid scaling 
criterion for particle size distribution as well as mean particle size for 
geometrically similar vessels with a volume range of 4 to 1. 

2. Binodal particle size distributions were observed. Although the 
coarse fraction varied as the —6/5 power of the agitation rate, the fine 
size was relatively independent of agitation rate. 

3. The system with the greater drop viscosity gave a smaller fine-frac- 
tion and a larger coarse-fraction than the other system (IB/W). 

4. Both systems approached equilibrium size slowly (many minutes). 
Equilibrium was approached more rapidly at high agitation rates. Dif- 
ferences were small between different tank sizes at equal energy dissipation 
rates per unit mass. 

5. The average size of the fine fraction remains about the same through 
the run (although the amount increases). The size of the coarse fraction 
decreases with time. More experiments are necessary to elucidate the 
mechanism behind this. 


Introduction 


We became interested in this problem because of practical considerations 
in emulsion and suspension polymerization. Often it is important for the 
product to have the right particle size or distribution of sizes. Yet the 
mechanism of breakup and the factors affecting size are little understood. 
Very little is known about the affect of physical properties, e.g., viscosity, 
interfacial tension, density, etc. Even less is known about the cause 
of variation from one piece of apparatus to another. 

This is essentially a progress report of two phases of our work: 


1. Equilibrium dispersions in which agitation is carried on long enough so 
that the particle sizes have reached some distribution which does not vary 
further with time. These dispersions were produced in geometrically simi- 


* Now in service, U. 8. Army. 
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lar tanks but of different agitator speeds selected to give equal power input 
per unit mass. 


2. Transient dispersions in which the agitation was started with the 
stratified phases, and the change in particle size distribution with time was 
observed. The same apparatus was used. 


Particle Size-Measurement 


The use of light scattering seemed to be the most advantageous method of 
measuring drop sizes and size distributions of unstabilized dispersions hav- 
ing drop radii of the order of 0.2-50 uw. The technique has been described 
by Sloan, Wortz, and Arrington.? 

There is a relationship between the intensity of light / scattered at any 
angle, @ from the incident beam and the size of the particle that scatters it. 
Sloan has defined conditions under which a plot of /@? versus @ can be used 
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to obtain information about the average particle size and the distribution 
of particle sizes. 

The method requires that certain conditions be met: (1) There is no 
angular dependent absorption of the incident light; (2) Secondary scatter- 
ing is not significant; (3) The particle radius is in the region of 0.1—-100 yu; 
(4) The ratio of refractive indices of continuous and drop phases should 
be not too far from unity. 
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Fig. 2. Agitation tank. 


Sloan has defined conditions under which a plot of /@? versus logarithm 
6 for a monodisperse system has the appearance of curve A in Figure 1. 
The location of the peak of the curve along the abscissa is related to the 
diameter d by the relation d@ = J, where K = a constant, which according 
to Arrington should be 19.4 for the conditions of this work and for spherical 
particles. Curves B, C, D, and F indicate the progressive change of the 
light scattering curves as the dispersion becomes more polydisperse. A 
change in the average size would be accompanied by a shift of the peak 
of the curve along the abscissa. The occurrence of two peaks in a single 
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curve, as in F, would indicate that the dispersion contained two maxima of 
rather narrow particle size in its distribution curve. 

The Sloan method of plotting gives, then, curves which change in a 
manner similar to distribution curves, but they are not true size distribu- 
tion curves. The location of the peaks gives a reliable indication of average 
size. The flatness of the curves gives only semiquantitative indication of 
the distribution of sizes, but such information was adequate for this work. 

To obtain a more adequate idea of what was happening in the transient 
system, the Wortz differential analysis of the J versus @ curve was per- 
formed. Wortz has used this method to examine five distributions pre- 
viously examined by sedimentation methods and obtained very good 


agreement. 


Apparatus 


The experimental apparatus consisted of an agitation system, a contin- 
uous circulating sample system, and a light-scattering apparatus. 

Three sizes of geometrically similar agitation vessels of the conventional 
design often used in industry were used to observe a small range of scale-up 
effects. The tanks were of the baffled cylindrical type with diameters of 
7.5, 9.0, and 12 in. The agitators were turbines consisting of six vertical, 
flat, rectangular blades. (lig. 2.) The drive mechanism for the agitators 
had a variable-speed transmission and was mounted so that reaction 
torque could be measured. Power input to the fluid could be determined. 

The circulation system consisted of a sampling probe, tubes, flowmeter, 
and return tube. A bellows pump was used downstream from the cell in 
the equilibrium tests. In an effort to get still lower shear, a peristaltic- 
action rubber-tube pump was used in the transient studies. 
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Fig. 3. Projected and partially cut-a-way views of light scattering cell: (A) along 
light beam (8) horizontally across light beam; (D) exaggerated length of wet length 
path. 











| 
| 





The light scattering apparatus used in the equilibrium tests was from 
Monsanto Chemical Company‘ and that in the transient tests was a com- 


mercial Scattermaster. Figure 3 shows the details of construction of the 


scattering cell. 


Both apparatus gave essentially the same diameter on a latex of known size 
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Fig. 4. Light-scattering curves: (A) polystyrene latex (radius 0.6 uw); (8B) mulberry 
pollen (radius 6.5 w); (C) mixture. 


(1.171 + 0.013 uw) curve A, Figure 4 and a mulberry pollen (15 u) curve B, 
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The length of the light path through the cell was adjustable. 
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Figure 4. Curve C is for a mixture. 


The sample transit time between tank and cell is about twelve seconds. 
It was concluded from several kinds of tests that the sample in the cell is 
identical with that in the tank. Some of the larger drop sizes in the cell 
are broken up in recirculating pump before going back to the tank, but 


that this bias is uniform and small. 
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Fig. 5. Light-scattering curves at low speeds: (A) IB/W, 7.5 in. tank I.D.; (B) 
IB/W, 9.0 in. tank 1.D.; (C) IB/W, 12.0 in. tank I.D.; -(D) CH/W, 7.5 in. tank I.D.; 
(E) CH/W, 9.0 in. tank I.D.; (F) CH/W, 12 in. tank [.D. 


Materials 


Two different immiscible liquid systems were used. The first (IB/W) 
consisted of 2-butanol and water. The second (CH/W) consisted of cyclo- 
hexanol and water. These are both partly miscible and were selected to 
minimize interfacial tensions and phase density differences. This was 
considered to be helpful in minimizing drop coalescence in the circulating 
system. 

In both systems, the organic-rich phase was the dispersed phase (2.5% 
volume in equilibrium study and 0.5% in transient study) and the water- 


rich was the continuous phase. 


Equilibrium Dispersion Results 


With the IB/W system in the smallest tank, a low agitation speed was 
selected which was slightly above the minimum necessary to maintain 
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Fig. 6. Light-scattering curves for IB/W system (tank 7.5 in. I.D.); (A) theoretical 
monodisperse system with same average size as (B); (B) low agitator speed; (C) high 
agitator speed. 


dispersion. When the dispersion had equilibrated at this speed, light- 
scattering data were obtained. The agitator speed was then doubled and, 
after re-equilibration, a new light scattering curve was obtained. [or the 
two other tanks, the high and low speeds were selected to give the same 
energy dissipation rates per unit mass. These same experiments were 
repeated with the IB/W system using the same agitator speeds. 

lor illustration, results obtained with the three tanks at a single energy 
dissipation rate are compared in Figure 5. Similar results were obtained 
at the higher rates. A typical difference between high and low energy dis- 
sipation rates is illustrated in Figure 6. Most curves show two maxima 
corresponding to a coarse and a fine size range. Tor comparison, a Sloan 
curve for a monodisperse system is shown. 

A relationship reported by Vermeulen® and others* indicates that the 
diameters of sufficiently large drops vary as the —6/5 power of the agitator 
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Fig. 7. Typical light-scattering data during transient period. 


speed. Consequently, doubling the speed should decrease the diameter 
56.5%. The largest decrease in the mean diameter of the fine fraction 
was only from 3.64 to 2.86 u or only about 21%. In IB/W system where 
the coarse maxima are clearly evident, the corresponding shift of the coarse 
fraction is from about 33 to 15 uw or 54%, which is certainly better agree- 
ment. 

The microscale of turbulence 7 varied from 30 to 70 u in these experi- 
ments.* The —6/5 power should hold for drop diameters greater than the 


* The microscale of turbulence is a measure of the size below which viscous forces be- 
come greater than inertial forces. It is calculated from the kinematic viscosity v and the 





energy dissipation rate per unit mass ¢ according to the equation: 7 = av, */e. 
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Fig. 8. Cross plot of scattering data of constant time. (IB/W) 12” tank. Intermediate 
agitation. 


microscale of turbulence. Since this is far from the case with the fine 
fraction, this may explain the difference in behavior. 


Transient Results 


The period before the attainment of equilibrium particle size was studied 
not only for its practical value but also to provide, hopefully, some insight 
into the mechanism of drop breakup and coalescence. 

Runs were made with the IB/W system at three energy dissipation rates 
per unit mass: about 1.0, 3.2, and 7.2 ft.?/sec.* (1 ft.?/sec.* = 0.929 X 10° 
ergs/g. sec.). Runs were also made with the CH/W system at the lowest 
and highest energy rates in the two larger tanks. Reynolds numbers 
varied between 23,000 and 82,000, and the microscale of turbulence varied 
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Fig. 9. 16? versus @ curve for IB/W system in 12” tank. Intermediate agitation. 


between 37 and 70 uw. It would seem that conditions were favorable for 
turbulence which was locally isotropic? but not necessarily homogeneous. 

Since large changes occurred in the first minutes, rapid and accurate 
scanning was required. The intensity was read at only certain critical 
angles, so that a scan could be made in about two minutes. A plot of time 
against corrected intensity was made for each angle (lig. 7) and then cross 
plotted to give curves of intensity versus scattering angle at specific times 
(Fig. 8). This plot then formed the basis for the Sloan and Wortz analyses 
(Figs. 9 and 10). 

Figure 11 shows a typical Wortz analysis of the data. Drop size as 
ordinate is plotted against time from the beginning of the dispersion 
processing, with the cumulative per cent below a given size as a parameter. 
lor example, after four minutes, 90% by weight of the dispersed phase had 
a drop size of 28 u or less. 

The data when interpreted by the Sloan method reveal two peaks indi- 
cating a binodal size distribution (Fig. 9). Unfortunately, the coarse peak 
is more difficult to identify in the later stages « a run because it tends 
to be obscured by scattering from increasing amounts of fine material. 





PARTICLE SIZE DISTRIBUTION 


However, the plot indicates that the average size of the fine material re- 
mains about constant while that of the coarse decreases. 

Comparison of the CH/W with the IB/W system indicates that the effect 
of increased drop-phase viscosity (CH/W) is to decrease the average drop 
size of the fine fraction and increase the average size of the coarse fraction 
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Fig. 10. Size distribution during agitation. Log-probability plot. 12” tank. 
Intermediate agitation. 


at equilibrium. The viscosities of the aqueous continuous phases were ap- 
proximately the same in both systems. Differences in interfacial tensions 
and densities were relatively smaller than the drop viscosity difference, 
(25.5 eycles ‘sec. for CH/W versus 2.39 eycles/sec. for IB/W). 

Comparison of results in one tank at different speeds indicates that time 
to reach equilibrium size decreases as speed increases and that the total 
energy required to produce a given distribution may be constant and inde- 
pendent of speed. 





LINDSEY, CHAPPELEAR, SULLIVAN, AND AUGSTKALNS 


” 
z 
°o 
lo 
oO 
= 


MINUTES 


Fig. 11. Change in cumulative undersize distribution during agitation. IB/W. 12” 
tank. High agitation. 
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On Phenomenological Rheo-Optic Constitutive 


Relations 


KLLIS HAROLD DILL, Department of Aeronautics and Astronautics, 
University of Washington, Seattle, Washington 


Synopsis 
The birefringence exhibited by high polymers can be explained by considering light 
as an electromagnetic wave propagating through an anisotropic dielectric. The ob- 
served dependence of the birefringence on the deformation state is described by ex- 
pressing the dielectric tensor as a functional of the strain history. For small strains, the 
functional can be approximated by an integral so that the birefringence and strain are 
related by a law of the Boltzmann superposition integral form. 


Introduction 


The birefringence exhibited by high polymers can be explained by con- 
sidering light as an electromagnetic wave propagating through an ideal 
anisotropic dielectric. The nature of the wave is determined by the rela- 
tion* between the free charge potential (electric displacement) D and the 
electric field E. For a material at rest 


D; = 60K 4 Lf; (1) 


Here « is the dielectric coefficient in free space and K ,; will be called the di- 
electric tensor. 

Suppose that the z;-axis is everywhere a principal axis of the dielectric 
tensor. The electromagnetic field equations show! that a plane wave 
propagated along the z; axis is divided into two linearly polarized waves. 
The direction of vibration of each wave coincides with one of the remaining 
principal axes of the dielectric tensor and has a refractive index equal to the 
square root of the corresponding principal value of the dielectric tensor. 

This suggests the introduction of a new tensor N,; which has the same 
principal axis as the dielectric tensor but which has principal values N; 
equal to the square root of the principal values of the dielectric tensor: 


= me ‘ 
Ni; = | (2 
The birefringence of a wave propagated along the z; axis is 
A = N, = Nz (3) 


* All quantities are referred to a fixed rectangular cartesian system z;. All indices 
have the range 1, 2, 3 and a repeated index indicates summation over this range. 
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The N,,; will be called the refraction tensor since the principal values of 
N,, are the refractive indices of the propagated wave. The refraction ten- 
sor is regarded, here, as the fundamental optical quantity in the formula- 
tion of constitutive relations. Other tensors which have been used in the 
formulation of constitutive relations are the velocity tensor.? 


Vy = cK,” (4) 

















and the Fresnel tensor® 


(5) 







Ny = Kym 





The birefringence has also been used directly‘ in the formulation of constitu- 
tive relations but, since the birefringence is not a tensor quantity, this proce- 
dure is unsatisfactory for a general treatment. 











Constitutive Relations 





The birefringence of high polymers is observed to depend upon the 
history of mechanical states. Since the refraction tensor N ,;, completely de- 
termines the birefringence, one might express this dependence mathe- 
matically by « relation between N,,; and the history of mechanical states. 
Consider an isothermal process at temperature %. Let ¢ denote the pres- 
ent time and 7 denote a past instant. The history of deformation is com- 
pletely described® by the strain tensor e;;(7),0 S 7 St. If N,,(0 is uniquely 
determined by the history of deformation, then N,;(é) must be a functional 
of e;;(7) and a function of e;;(t): 
N,(6) = Pasles'(r); €:;(t) ] (6) 


r= 











tj 






The functional may be expanded in a power series :*!4 


t 
Fyle(r); eO)] = file) + f Jijrs(€(L), t, r)e;.(r)dr 
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The functions may also be expanded in a power series: 


fisle) = Niuj® + Giyrs€rs + Vijremn€rs8mn 7A (8) 


7 tj 












ete. 
If only small strains are considered, then all series can be truncated so 





that eq. (6) becomes linear in the strains: 


t 
N,(t) = NifP + agzce;s(t) + { A jirs(tyr)e,.(r)dr (9) 
J ¢ 





Since the constitutive relation should be invariant with respect to a shift 









of the time origin, the functions A ;, must be of the form A; ,(t — 7). 
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Therefore, the constitutive relation for a anisotropic material experiencing 


small, isothermal deformations is 


nt 
Ni() = Ny? + J Bie(t — 1)(d/dee-(1)dr (10) 
0 


where 
B ijx(0) = Qijer 
(dj dt) Biyr(t) = A ijpr(t) 
Since E is symmetric, there is no loss in generality if B is symmetric in the 
last two indices. Assume that N,, is symmetric then B is also symmetric 


in the first two indices. 
The material will be called isetropic if B is an isotropic tensor :? 
Bijer —_ Ab bxy + Bb 46 jr + 6:76 jr) (11) 
where 6;; is the Kronecker data. Since N,,; = N°6;;, the constitutive rela- 
tion for an isotropic material experiencing small, isothermal deformations is 


t 
Ny = N% + fia (t — r)6,,(d/dr) [e,,(r) ] 
0 
+ B(t — r)2(d/dr)[e,(r)|}dr (12) 


Equation (12) has the usual form of the convolution integral. All of the 
well-known results which have been of use in viscoelasticity® can be imme- 
diately applied. In particular, the functions A(¢) and B(¢) can be repre- 
sented in terms of their spectrum. In the case of a discrete spectrum, the 
constitutive relation can be related to a ‘‘spring-dash pot”’ model. 

Denote the deviatoric part of the refraction tensor by 

T | T 
Ry = Ny — = Norby 
3 
and the deviatoric part of the strain tensor by 
l 


egy = 


Then eq. (12) yields 


t 
R(t) = f2Be — r)(d/dr) [e;;(7) |dr (14) 
( 


) 


A relaxation process is defined as the response to a step change in strain. 
If ¢,; has the step change ¢,,° at ¢ = 0, then eq. (14) gives 


R,(t) = 2B(e,) . (15) 


Therefore B(t) will be called the optical relaxation function. 
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Suppose that it is found experimentally that (dB/dt) varies monotonically 
and that B(t) approaches a limit B.. when t— ©. Then define a function 
¢(t) such that 


Bit) = Bo + (B. — Boo) (16) 


Then ¢(t) increases from zero to unity with monotonic slope (Fig. 1) and 
will be called the reduced relaxation function. 


o,f 
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Fig. 1. Reduced relaxation function. 


Consider isothermal relaxation at a different temperature 6. Denote the 
optical relaxation function by F(¢) and define a reduced relaxation function 
f(®) analogous to eq. (16). Then there exists a function &(7), which may 
be called the reduced time (Fig. 1), such that 

I(r) = (&(7)) (17) 

In the case of stress relaxation of a linear viscoelastic material, it is found 

that 

E(r) = a(6)r (18) 
The function a(@) is termed the temperature shift factor and depends upon 
the choice of reference temperature %. Assume that eq. (18) also holds 
for isothermal optical relaxation but perhaps with a shift factor different 
from the mechanical case. 

Consider now relaxation at varying temperature 6(/)._ Again denote the 
optical relaxation by F(t) and introduce a function f(t) defined analogous 
to eq. (16). Then again there exists a reduced time &(t) such that eq. (17) 
holds. Therefore 


(df/dt) = (dé/dt)(do/dé) (19) 
In the isothermal case 
(df/dt) = a(6)(do/dé) 
This suggests the hypothesis 


(dé/dt) = a(6(t)) . (20) 
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and 
t 
é(t) = f aocs)ar (21) 
0 


in the case of varying temperature. 

Now suppose that the temperature varies in the same way, @(t), but the 
strain increment occurs at time 7. _ Denote the resulting creep function by 
F(t,r). In this case the fundamental constitutive relation, eq. (14), must 
be written 


t 
Ry = fe F(t,r) (d/dr) [€:;(7) |dr (22) 
0 


that is the kernel can not be expressed as a function of the difference ¢ — r. 

Analogous to the isothermal case, let F(t,7) = Fo and assume that 
[OF (t,r) ]/(Ot) varies monotonically with ¢ while F(t,r) > F. ast —~> o~. 
Then the function f(t,7) may again be defined as in eq. (16): 


F(t,r) = Fo + (F.. — Fo)f(t,7) (23) 
In this case there exists a reduced time p(/,7) such that 
S(t,7) = 6(e(t,7)) (24) 
and 
(Of/Ot) = (Op/dt)(0d/Op) (25) 
The similarity to eq. (19) suggests that 
(Op/dt) = a(@(t)) (26) 


Since f(t,7) = o(0) = 0, p(t,7) = 0. Therefore 
t 
p(i,r) = a(O(n))dn (27) 


= &(t) — &(r) 
Therefore the creep function at varying temperature has the form 
Fit,r) = Fo + (F.. — Fo o(E® — E(r)) (28) 
If Fy and F., are independent of temperature, then | 
F(t,r) = BE — &7)) (29) 


In this latter case, the fundamental constitutive relation for varying tem- 
perature, eq. (22), takes the form 


& 
Ry*(®) = f OB(E — 2) (d/dn) [en,*(n) In (30) 
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~! 
to 


where 
Ri* (EO) a R(t) 
€4* (E(0) 7 €:;(t)) 

Similar rezsoning may be applied to the trace of N4,. 

The eqs. (10), (12), and (30) are the fundamental constitutive equations 
in the various special instances. Similar results can be derived if one prefers 
to relate the birefringence to stress. Moreover, since the stress and 
strain histories are uniquely related for most materials, the fundamental 
equations here can usually be transformed into stress-birefringence rela- 


tions. 


Linear Viscoelasticity 


The constitutive equations for materials such that the stress depends on 
the deformation history may also be expressed in terms of functionals.®' 
For small isothermal deformations,’ 


t 
oi(t) = f 2G(t — 1r)(d/dr) [€;;(7) |dr (31) 
( 


) 
where 6 is the deviatoric part of the stress tensor. A similar expression® 
can be derived for the mean stress. 
Equations (14) and (31) imply 


t 
Ry = f C(t — 1r)(d/dr) [o4;(7) |dr (32) 
0 


Therefore a unique relation between birefringence and strain history implies 
a unique relation between birefringence and stress history in the case of 
linear viscoelastic materials. 

In the ease of varying temperature, hypotheses analogous to those used 
above concerning the temperature shift factor and the reduced time can be 
used'! to derive a constitutive relation for stress. In the notation of eq. 
(30), 


€ 
o*,(£) = f 20(— — n)(d ‘dn) [€:;*() |dn (33) 
0 


If the mechanical shift factor is the same as the optical shift factor, as 
supposed here, then eq. (30) and eq. (33) yield 


€ 
R,*() = f C(é — 9)(d/dn) [oy*(n) len (34) 
C 


) 


Plane Polariscope 


Now consider an idealization of the situation in a plane polariscope. A 
slice of the dielectric is bounded by the planes z; = 0,h. The z;-axis is 
everywhere a principal axis of K,, and therefore of Ny. By means of a 
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polarizer, a linearly polarized plane electromagnetic wave propagating 
along the z;-axis is obtained. Every such wave can be divided into the sum 
of monochromatic waves and each monochromatic wave can be decomposed 
into components such that the direction of vibration coincides with one of 
the principal axis of N,;. 

Consider one such wave. When this wave is incident on the dielectric 
it is partly reflected and partly transmitted. A similar reflection and trans- 
mission occurs at the exit surface. By applying the jump conditions of the 
electromagnetic field equations the amplitude of the exited wave can be 
related to the amplitude of the incident wave. If the exited wave is now 
viewed through a polarizer which has its axis at right angles to that of the 
first polarizer, the amplitude of the transmitted light is approximately!” 
proportional to the factor 


I = sin 2a sin 2nr (35) 


where a(t) is the angle between the axis of the polarizer and the first princi- 
pal axis of N,;, and 


n = (hw/2mc)A (36) 


is the fringe order. Here, h is the model thickness, w is the circular fre- 
quency, and c is the velocity of light in free space. 

Thus, if the model is subjected to a non-homogeneous deformation field, 
there will appear lines of zero light intensity (isoclinics) connecting the 
points where a(t) = 0, (2/2). There will appear lines of relative minimum 
intensity where the fringe order is an integer. Therefore by observing the 
isoclinie and isochromatic pattern, the functions a(t) and A(?) can be de- 
termined. 

Using the inverse of eq. (32), 


1¢ 
oi(t) = =f Y(t — 1r)(d/dr)[Rjj(7) |dr (37) 
2J0 


it follows that 


1 rt 
out) — a(t) a. y(t — 1r)(d/dr)[A(r) cos 2a(r) |dr 
2 Jo 


(38) 
1 t 
o(t) = =f v(t — 1r)(d/dr)[A(r) sin 2a(r) |dr 
“ 0 
The principal axes of stress and the difference in principal stresses can 


then be calculated. Similar results apply to the strain tensor. 
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Résumé 


La biréfringence présentée par les hauts polyméres peut s’expliquer en considérant 
la lumiére comme une onde électromagnétique se propageant 4 travers un diélectrique 
anisotrope. La dépendance de la biréfringence vis-d-vis de l'état de déformation est 
décrite en exprimant le tenseur diélectrique en fonction de ]’évolution de Ja tension 
interne. Pour de faibles tensions, on peut décrire approximativement la fonction par 
une intégrale de telle sorte que la biréfringence et la tension sont reliées par une loi du 
type intégrale de superposition de Boltzmann. 


Zusammenfassung 


Die Doppelbrechung von Hochpolymeren kann durch Betrachtung des Lichtes als 
eine elektromagnetische, sich in einem anisotropen Dielektrikum ausbreitende Welle 
verstanden werden. Die beobachtete Abhingigkeit der Doppelbrechung vom Deformie- 
rungszustand wird durch Darstellung des dielektrischen Tensors als ein Funktional der 
Verformungsvorgeschichte beschrieben. Fiir geringe Verformung kann das Funktional 
durch ein Integral angeninert werden, sodass Doppelbrechung und Verformung durch 
ein Gesetzt von der Form des Boltzmannschen Superpositionsintegrals verbunden sind. 
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Stress-Strain-Time-Birefringence Relations 


in Photoelastic Plastics with Creep 


CHINTAKINDI L. AMBA-RAO, School of Aeronautical and Engineering 
Sciences, Purdue University, West Lafayette, Indiana 


Synopsis 


The above is an exploratory study started in September 1961 with a view to solving 
the time dependent stress and strain distributions in structural problems through optical 
model analysis. The present analysis, however, is restricted to monotonically increasing 
uniaxial loading. Finally, “optical creep parallels mechanical creep’’ is the basis on 
which this analysis is built. 


Introduction 
Purpose 


The aim of the author in deriving the stress-strain-optice law for plastics 
in the presence of creep is to extend the conventional photoelastic method 
of stress analysis to the solutions of static structural problems in the pres- 
ence of creep. Earlier attempts made in this area are similar in nature; al- 
though the problems are different, they are closely allied. There is a classi- 
cal paper by Mindlin' and a recent paper by Theocaris and Mylonas;? the 
title of the former is rather misleading. Both papers deal with ‘‘the condi- 
tions under which models or birefringent coatings made of viscoelastic ma- 
terials may be employed to solve boundary value problems in the linear 
theory of elasticity”’ through photoelastic principles. 


Various Methods of Specifying Viscoelastic Properties 


According to Alfrey,* there are no less than seven distinct (but all mathe- 
matically equivalent) methods of specifying the properties of linear visco- 
elastic materials, divided into two classes. 


Class I 


(A) The Voigt model, consisting of a set of Voigt elements in series. 
(B) The Maxwell model, consisting of a set of Maxwell elements in 
parallel. 
(C) The operator equation, P(o) = Q(e) where P and Q are linear dif- 
ferential operators. 
(D) The mechanical impedance function Z(w). 
75 
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Class IT 


(£2) The creep curve—strain as a function of time at constant stress. 

(Ff) The relaxation curve—stress as a function of time at constant 
strain. 

(G) The dynamic modulus curve—elastic modulus as a function of fre- 
quency. 

The class IT methods are ideally suited for the development of an analysis 
which basically rests upon experimental observations, while class I methods 
are suitable for obtaining solutions of theoretical problems; and each 
method has its own limitations. 

Since birefringence is a consequence of the material being stressed (or 
strained), the author is inclined to believe that optical anisotropy is a con- 
sequence of mechanical anisotropy. Hence, the approach to the field of 
optical viscoelasticity may be made on lines similar to those of mechanical 


viscoelasticity. 
Survey of Literature 


The following lines from the doctoral dissertation (University of Illinois, 
1957) of Dr. Daniel Post* seem to be appropriate at this stage. 

“The stress optic law is well established but even today, the validity of 
the law rests largely upon experimental verification.4 There has been no 
completely satisfactory explanation of the stress-optical phenomenon as 
exhibited by amorphous solids in terms of the electromagnetic character 
of light and the electric and magnetic properties of the stressed solid.’’>* 

The author’s approach to the present problem is largely based on experi- 
mentally observed facts in an engineering laboratory. He is of the opinion 
that there is no unique way in which the stress-strain-time-birefringence re- 
lation can be derived in the presence of creep as can be seen from the 
published work. ‘Besides, statistical mechanics applied to the optical 
properties of the whole network on the lines of Kuhn and Griin’? and work 
on the plastics at the transition temperature may also yield fruitful results. 


The Optic Law in Elasticity 


The optic law in classical photoelasticity takes the following simple 


form: 
(n /d) = C,( 01 = do) 
(1) 
= ce. — €) 
where 
c.=c,E/1+» (2) 


* Dr. Post is a noted authority on optics applied to photoelasticity and is widely known 
for his fringe sharpening and fringe multiplication techniques. 
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~] 
n_ 


and where n = relative retardation, d = thickness of specimen, c, = stress 
optical coefficient, c, = strain-optical coefficient, 4; — o2 = principal stress 
difference, ¢, — €¢ = principal strain difference, / = Young’s modulus, and 
v = Poisson’s ratio. 

Within the elastic range, the stress is proportional to strain; hence, it is 
immaterial whether one considers the optical birefringence as being due to 
principal stress difference or strain difference. However, when the medium 
creeps (call it viscoelastic), the optic law takes a much more complicated 
form. 


The Optic Law for a Medium Which Creeps 
A historical review of the work includes the following: 
(n/d) = c.(o1 — o2) + ¢(€e1 — €2) (3) 


Coker and Filon.* c, varied from 5.4 to 9.6 Brewsters. c, varied from 0.4 
to 1.58 Brewsters. 


(n/d) = c,(t)(o1 — a2) (4) 
Durelli and Lake.® 
(n/d) = c,(t)(o1 — a2) (5) 
Theocaris and Mylonas.* 
1 
K, = =K,+at+ blogt (6) 
Co 


where k,, k,, a, b are empirical constants. 
(n/d) = Cq@(o1 — 92) ip) + Ci(o1 — 92) ¢) (7) 


Where c,») and c@) refer to the optical coefficients of the plastic phase and 
the elastic phase (Kuske™). ‘Though chemists question the multiphase 
nature of plastics,’’ Kuske experimentally verified his law on many ma- 
terials. 

Recently, Dill’! has put forward a hypothesis based on electromagnetic 
theory and Boltzmann’s superposition principle. Aithough it has some 
very interesting characteristics, it has a restrictive assumption that in the 
form presented it is valid only for linear viscoelastic materials. Professor 
Drucker of Brown University in a general lecture at the Fourth U. S. 
National Congress of Applied Mechanics pointed out that linear visco- 
elastic problems are far from reality. It is known that most of the metals 
at high temperature and plastics at high stress or strain levels do not be- 
have as linear viscoelastic materials. Crystallizable polymers and incom- 
pletely polymerized resins change their properties during creep and do not 
satisfy the superposition principle. This is an observation made by the 
author, following the work of Leaderman!? on Boltzmann Principle 


- 
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The Analysis 








Theoretical Foundation 





Guided by the existing knowledge, the present analysis is based strictly 
on experimentally observed facts. In the present analysis, the author con- 
centrated on class II methods suggested by Alfrey. Relaxation and creep 
tests are suitable for viscoelastic studies in the time scale of 1 to 10° sec., 
while dynamic modulus measurements are suitable for the time scale of 1 
to 10-* sec. Since the author’s main interests lie in static problems, the 
latter approach is ruled out. Main concentration is made on creep tests as 
no Instron Testers are available to perform relaxation tests. 














Creep Tests under Constant Load 





Creep tests under constant load are assumed to be the same as the creep 
tests under constant stress as there is negligible change in the cross-sectional 
area. The following analysis lays a theoretical foundation for the study of 
some model materials. A series of assumptions are made to start with, and 
the ultimate justification for the assumptions depends upon the final re- 
sults. 

Let it be assumed that 


€-(t) 












= filo) = A(o)g) 
= a7g(t) 












p = 1, if it has linear mechanical creep. 







e.(t) = e(t) — «, and ¢,(t) = € 


Dot denotes differentiation with respect to time. Similarly 







n(t) = n()) — n 


and 







and 





(9) 


where subscripts ¢ and e refer to creep and elastic components. In the 
absence of any specific information whether birefringence is due to stress or 
strain, it is reasonable to assume that birefringence is a function of both 







n,(t) = So(o,¢,l) 









stress and strain. 
Substituting eq. (8) in eq. (9), one obtains 


n. = f3(o,t) = Bio) fd) 
ov (t) 


q = 1, if it has linear optical creep. 
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It has been experimentally verified that g(t) and f(t) are similar in nature 
for many model materials. Equation (8) is valid for many materials for the 
primary and secondary creep stages of the creep curve. Rabotnov, a well- 
known authority on creep says in a recent paper,'* “when we take into 
account the great discrepancies in experimental data on creep, we see that it 
is rather pointless to strive for precision in such predictions. Therefore, 
improvements in the equations of creep are often held to be meaningless, 
and it is frequently suggested that one apply the simplest theories.”’ (The 
emphasis is made by the author.) 

The author fully concurs with Rabotnov’s remarks and believes in simple 
expressions for creep curves. 


Optical Creep Parallels Mechanical Creep 
Experiments by Maxwell" indicate that the relationship between change 
of index of refraction and stresses is the same as that between strain and 
stresses for an elastic medium. 
eS No = (0; +- Co( a2 +- 3) 
ne — No = cCro2 + C2(03 + 01) (11) 
ns; — No = Cio3 + C2(o1 + 2) 
where 7, 72, 23 are the three principal indices of refraction along the 1, 2, 3, 
directions and N,j is the radius of the index sphere before stress is applied. 
v 
€1 E (a2 + 93) 


v 


E (a3 + 01) 


€3 7 (o1 + 2) 

The creep behavior is considered to be a small perturbation from the 
elastic behavior at ¢ = 0. Hence, it is reasonable to assume that the 
mechanical creep strain parallels the optical creep birefringence.* This is 
the idea underlying the analysis developed in this report. From the litera- 
ture one could collect the following evidences. 

On Catalin, up to 10 min. after loading, Mylonas" obtained the relations: 


n(t) = o(A + Bt”) (13) 
e(t) = ofA’ + B't”) (14) 
where A, B, A’, B’ are constants. 
Coker and Filon® fitted the following equation (up to an hour after load- 
ing) for xylonite (celluloid) : © 
n-(t) = pt’? + ut (15) 
e(t) = al’ + bi (16) 
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where a and p are proportional to o? approximately; and b and q “‘are very 
small quantities not accurately determined.” (The emphasis is made by 
the author.) 


For Xylonite, eqs. (15) and (16) approximately reduce to 


n(t) = Ky,ot'” (15a) 
i>0 
e(t) = Kyot'” (16a) 
lor CR-39, the author found that 
n(t) = af®-0308 (17) 
t>0 
€(t) = af.0356 (18) 


Photoviscoelastic Law 


Choose a model material which creeps. Prepare optical and mechanical 
creep curves in uniaxial tension. In order to enable master curves to be 
prepared, eqs. (8) and (10) may be written as follows: 


le-(t)/A(o)] = g (19) 
[n-(t)/B(o)] = fd) (20) 
Linear viscoelastic model material is a particular case, where 


A(c) = Bic) = o 


One can easily see from eqs. (19) and (20) that, in the most general case, 


m(l)\ _ e,(L) ei , 
(ss) a (< e) + Constant (21) 


The justification for eq. (21) is the fact that for Xylonite and Catalin, in- 
vestigated upon and reported in the literature, and, from the experimental 
observations of the author on other materials," “the curve between the opti- 
cal creep strain per unit stress and time is of the same (or similar) type as 


120 140 160 
t 





MIN 
Fig. 1. «(t)/o and n(t)o versus time for Araldite 502 at 74°F. 
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40 60 80 100 120 149 160 180 
MIN t 
Fig. 2. ¢(t)/o and n(t)/o versus time for CR-39 at 72°F. 


the curve between the mechanical creep strain per unit stress and time.”’ 
This can be seen from the expressions for the optical and mechanical creep 
curves and the experimental results plotted in Figures 1 and 2. In Figures 1 
and 2, (n(t)/c) and (e(t)/o) are plotted against ¢ for CR-39 and Araldite 
502. They are found to be coincident within +3.5% from the mean, basing 
the deviation on the creep components only. Similar coincidence has been 
obtained even for cases in which the stress is arbitrary, but uniaxial. 


Effect of Rate of Loading and Rate of straining 


Experiments were also conducted to find the dependence of mechanical 
and optical creep curves on different rates of loading. The effect of strain 
rate on stress-strain and stress-birefringence curves was also studied. The 
conclusions are as follows: For CR-39, the maximum différence in bire- 
fringence and strain between a constant load creep test and a variable load 
creep test under the worst conditions investigated is about 12%. 

Strain rate has an effect on stress-birefringence and stress-strain curves. 
This is in line with the behavior of viscoelastic materials. Strain rate was 
varied from about 6 to 600 yin./in./min. and the maximum variations in 
stress-birefringence and stress-strain curves at 1600 psi stress level are 20% 


9707 


on birefringence and 27% on strain. Therefore, 
ne = f,(¢,6 
eo = fr(G,€) 
Experimental Observations 
In the above analysis, if one makes p = v = 1 in eqs. (8) and (10) and 
combines them with the quantities in the elastic state, one gets the linear 


viscoelastic law. 
If (n(t)/o) plotted against (€(t)/o) is a straight line, then 


AO eT OAD et Lp 


og og og 
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That is, 
Ne +n,(t) ~ Ale, + «(t)] + Bo (22) 


It may be mentioned that A and B are constants while, n,(¢) and «,(t) de- 
pend upon the history and rate of loading. 
The above relation, eq. (22) is valid for the following ranges investigated : 


0 < ¢< 180 min. 
0 < « < 2,000 psi. 


Equation (22) is valid as a stress-strain-time-birefringence relation for some 
model materials which creep. 
On the basis of experiments conducted by the author on Araldite 502 
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Fig. 3. n(t)/o versus e(t)/o for Araldite 502 at 74°F. 
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Fig. 4. n(t)/o versus «(t)/o for CR-39 at 72°F. 
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FRINGES/p-s:‘i 


E(t) MICROIN/IN. 
o p-s.i 


Fig. 5. n(t)/o versus ¢«(t)/o for Catalin 800 at 19.5°C."5 


with 9% hardener HN 951, it is found that the relation between (n(t)/<c) 
versus (e(t)/c) is slightly nonlinear while, on the basis of experiments on 


CR-39, the author found the relation between (n(t)/c) versus (e€(t)/c)as 
linear. Both the model materials are found to behave as linearly visco- 
elastic in the range tested. The details of the experiments for CR-39 are 
given in earlier publications;'* and, Figures 3 and 4 indicate the rela- 
tionships for both the model materials. In Figure 3 and 4, 0.95 ¢« takes into 
consideration the stiffening effect of the strain gages in measuring the 
strains. 

For Catalin 800, the author has taken the experimental values from the 
literature and obtained 


n(t) = Ale(t) — e(t)|] + Bo 


for 0 < ¢< 10 min. and 261.7 < o < 963.4 psi (see Fig. 5). 

The above laws, eqs. (22) and (23) derived by the author for photovisco- 
elasticity agree in principle with the earlier work of Bayoumi and Frankl” 
on photoplasticity. The present work is done independently by the 
author and some of the results were presented as an IAS (Institute of Aero- 
space Sciences) student paper in April 1962. A recent paper of Javornicky” 
from Czechoslovakia has results which confirm the main observations re- 
ported in this paper. 


The Optical Boltzmann Principle 


It is well known that the Boltzmann Principle holds good for linear 
mechanical viscoelastic materials, under certain conditions; and it was 
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proved earlier in some typical model material that optical creep parallels 
mechanical creep. There is a great amount of literature on mechanical 
linear viscoelasticity; and, if one replaces « by n, one obtains the linear 
photoviscoelastic laws. The author experimentally verified the photo- 
viscoelastic law in the following form for CR-39 and Araldite 502 for arbi- 
trary uniaxial tensile loading 


; l 
n(t) = f [c, creep] (t — 7) wate) dr (24) 
0 dr 


where [C, creep] is the birefringence variation in time following a unit 
stress input, generally termed as the optical creep function. The percentage 
error was found to be less than 2%. The above experiments confirm the 
hypothesis proposed by Dill'! on the basis of electromagnetic theory. 


The Model Materials 


The prime purpose in using CR-39 and fully polymerized Araldite 502 as 
the model materials for the present analysis is the fact that their mechanical 
and optical properties are stable as contrasted to the plasticized, incom- 
pletely cured epoxies, prepared in the laboratory, whose properties vary 
substantially with time.?!_ CR-39, like most other model material, includ- 
ing fully polymerized epoxies, does not flow while, some incompletely 
(and/or plasticized) polymerized epoxies do flow. The great disadvantage 
with the latter is the fact that their properties vary from day to day and 
one is likely to get non-repeatable results. Creep properties are very 
sensitive to the composition and the curing of the specimen and the author 
found discrepancies as high as 10% even between two identical specimens 
from the same sheet. Hence, one is forced to use the same specimen for all 
the tests. This is not possible with a model material which flows, as some 
residual creep strain is left over after every test. Annealing, if resorted to, 
besides removing the residual strain, also changes the composition of the 
material. These are basically the difficulties which the research group at 
New York University”? has run into. 

CR-39, Columbia Resin, was a very popular model material for photo- 
elasticity all over the world, until the discovery of Epoxy resins in 1951, in 
spite of the fact that it creeps. It is linearly viscoelastic up to about 
2,000 psi stress level. 

The other model material with which the author worked is Araldite 502 
(also known as Araldite D) with 9% of hardener HN 951. The monomer 
belongs to the family of epoxy resins manufactured and marketed by Ciba, 
Inc. The author is one of the earliest to investigate the suitability of this 
resin for photoelastic investigations” in 1956. For the present investiga- 
tion, the author used a complicated heating and annealing cycle to sta- 
bilize the optical and mechanical properties as follows: Make the casting 
at room temperature; and, after fifteen days, subject it to annealing at 
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100°C. for five days, and store the specimen at room temperature for nine 
months. 


Conclusions and Possible Extensions 


In conclusion, the present analysis may be extended to the following 
cases: 


(a) Nonlinear creep problems. 

(b) Large (finite) deformation. 

(c) A(o) # Bia) (see eqs. (19) and (20)). 
(d) f(t) ¥ g(t) (see eqs. (19) and (20)). 


Work is in progress to extend the present analysis to biaxial stress states 
and arbitrary loadings. 

The above analysis is applied to predict the time-dependent stress and 
strain distributions in a rectangular beam in pure bending. The results 
will be published somewhere else. 


The author is grateful to Professor Ervin O. Stitz for serving as his thesis advisor during 
the course of this investigation. The author is indebted to the reviewer for valuable 
suggestions. 
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Résumé 


Ce travail constitue une étude préliminaire commencée en septembre 1961 en vue 
d’expliquer la dépendance du temps des distributions de tension et d’élongation dans le 
probléme de structure par une analyse de modéles optiques. La présente analyse se 
limite cependant A |’étude de systémes si la charge augmente de maniére monotone. 
Enfin, la mise en paralléle du fluage optique et du fluage mécanique est la base qui a 
servi A l’élaboration de cette étude. 


Zusammenfassung 


Die vorliegende orientierende Untersuchung wurde im September 1961 mit der Ab- 
sicht begonnen, die zeitabhiingige Spannungs- und Verformungsverteilung bei Struktur- 
problemen durch eine optische Modellanalyse zu klaren. Diese Analyse ist aber auf 
monoton wachsende einachsige Belastung beschriinkt. Schliesslich ist ‘“Optisches 
Kriechen geht parallel zu Mechanischem Kriechen”’ die Grundlage, auf der diese Analyse 
aufgebaut ist. 
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Dynamic Birefringence of Amorphous Polymers 


B. E. READ, Basic Physics Division, National Physical Laboratory, 
Teddington, Middlesex, England 


Synopsis 


The principles underlying the dynamic birefringence method are outlined, giving par- 
ticular emphasis to the use of this technique in investigating the molecular mechanism of 
relaxation in amorphous polymers. Data are presented for amorphous polyacetaldehyde 
({—CH—O—],,) in the frequency range 1-20 cycles/sec. and temperature range 0 to 


CH; 
—25°C. The results are discussed and compared briefly with dielectric data for poly- 


acetaldehyde and also with earlier dynamic birefringence data obtained for polymethyl- 
acrylate. 


Introduction 


The dynamic birefringence technique has recently emerged as a valuable 
new tool for investigating the molecular mechanism of mechanical re- 
laxation in polymers. In this method the real and imaginary components 
of both the complex strain-optical coefficient and the complex stress- 
optical coefficient are determined as a function of frequency. Such data 
may be compared with the corresponding frequency dependence of the 
components of the complex modulus or compliance. Furthermore the 
strain-optical and stress-optical relaxation spectra may be evaluated and 
compared with the mechanical relaxation and retardation spectra. Since 
birefringence provides a measure of molecular orientation, such comparisons 
may yield information concerning the molecular orientational rearrange- 
ments responsible for the mechanical relaxation behavior. 

In this laboratory dynamic birefringence studies have initially been 
made upon amorphous polymers in the relaxation region associated 
with their rubberglass transformation. The behavior in this region is 
generally attributed io the segmental motions of the polymer chain back- 
bone. Our experimental method has been described in some detail! 
and results have been reported for polymethylacrylate (PMA).? In this 
paper data are presented for amorphous polyacetaldehyde ({[—-CH—-O—],). 


CH; 
This polymer was chosen for investigation because the —-CH; side 
groups were expected to have less influence on the relaxation of bire- 
fringence than the —COOCH; side groups in the PMA chain. Further- 
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more, comparable dielectric results have also been obtained for poly- 
acetaldehyde in this laboratory. Owing to the presence of the dipole 
moment within the chain backbone, the dielectric relaxation of this polymer 


also arises from main chain motions. 


Theory 


Phenomenological 


For a polymer sample subjected to a sinusoidally varying tensile strain 
e we arbitrarily assume that the stress / and birefringence An cycles 
lead the prescribed strain cycle by phase angles 6 and 6’, respectively. 
Hence we have (in complex notation), 


e* = e« exp {iwt} 
f* = foexp{i(wt + 8)} (1) 
An* = Ano exp{i(wt + 6’)} 


where ¢€0, fo, and An are the strain, stress, and birefringence amplitudes, 
respectively, w the experimental frequeney in radians/sec. and ¢ is the time 
in seconds. In the present experiment ¢0, fo, Ano, 6 and 6’ are each meas- 
ured as a function of frequency and temperature. From these measure- 
ments we obtain the real and imaginary components of (a) the complex 
modulus (2* = f*/e* = EL’ + iE”), (b) the complex compliance (D* = 
e*/f* = D’ — iD"), (c) the complex strain-optical coefficient (K* = 
An*/«* = K’ + iK”") and (d) the complex stress-optical coefficient (C* = 
An*/f* = C’ — iC”), where for example, 


K’ = |K| coss’, K" = |R| sin 8’ (2) 
C’ = |C| cos(6 — 8’), Cc" = C| sin (6 — 6’) 7 


where |K| = Ano/e and |C| = Ano/fo. 

According to the phenomenological theories, the dependence of EF’ 
and E” on frequency may be expressed generally in terms of well-known 
expressions involving the mechanical relaxation spectrum ¢z(In 7). Also 
the frequency dependence of D’ and D” can be written in terms of the 
retardation spectrum @p(In 7)} We have similarly defined a strain-optical 
relaxation spectrum, ¢x(In 7), and also a stress-optical relaxation spec- 
trum, ¢c(In 7), which relate the components of K* and (*, respectively, to 
frequency. For example, 


+o 2,2 
K’ = . : ox(In 7) ae T 


+o 
Ks f - ox (In 7) aa din +r 


oT” 


where ¢x(In 7) represents the contribution to K’.. — K’» from processes 
having strain-optical relaxation times between In zandInz7+dInr. K’. 
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and K’y are the limiting values of K’ at infinite and zero frequency, re- 
spectively. An analogous definition has been given? for ¢,(In 7). Toa 
first order approximation, each of the above distribution functions may be 
valculated from the experimental data according to equations given.” 
It is to be emphasized that whereas ¢z(In 7) and ¢p(In 7) are always positive, 
¢x(In r) may be either positive or negative in different region of 7 depending 
on whether the mechanism with relaxation time 7 contributes a positive or 
negative amount to K’.. — K's. Hence it is possible for K’ to either 
increase or decrease with frequency in a particular frequency range de- 
pending on whether ¢x(In 7) is positive or negative for the appropriate 
relaxation mechanisms. 

The above phenomenological discussion refers to the dynamic mechanical 
and optical behavior of a polymer over a wide frequency range at con- 
stant temperature. Since our present measurements are limited to about 
1.5 decades of frequency, a method of reduced variables has been tentatively 
proposed for effectively extending this frequency range. This procedure 
involves an initial reduction of both the mechanical and optical data to a 
standard temperature 7’) as described.2, Subsequently the reduced experi- 
mental plots are shifted along the log (frequency) axis until they superpose. 
Master curves are thus obtained of E’7,, E"7,, D’r,, D"r,, K’r, K" x, 
C’,,, and C"7, against log arf where f is the frequency in c/s, and log ar 
is the horizontal shift factor. Since greater precision has so far been at- 
tained with the mechanical data than with the optical data the log a, 
values are usually determined from the mechanical plots and then used to 
construct the optical master curves. 


Molecular 


According to the proposed molecular theory,* the mechanical and optical 
relaxation arises from the diffusion of chain molecules toward a new set of 
most probable configurations subsequent to a deformation of the sample. 
Each real molecular chain is replaced by an idealized chain comprising a 
large number of freely jointed, optically anisotropic, links. The con- 
figurational rearrangements therefore lead to a positive entropy change but 
zero change in internal energy. The molecules are each subdivided into v 
submolecules each of which contains sufficient links so that its end-to-end 
distance conforms to the Gaussian distribution. The configuration of an 
entire chain is specified initially in terms of the Cartesian coordinates of 
one end of each submolecule relative to the other end as origin. 

A basic postulate underlying the theory is that when the specimen 
undergoes an instantaneous deformation, the components of the vectors 
joining the ends of each submolecule initially undergo a change in the 
same ratio as the change in dimensions of the bulk specimen (i.e., an affine 
deformation). If the sample is held at constant external strain the sub- 
molecules return to less extended configurations during the subsequent 
relaxation process. After transforming the original Cartesian coordinates 
to a system of normal coordinates, the variation of the birefringence and 
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stress with time (at constant strain) are calculated in terms of the time 
dependence of each normal coordinate. Under dynamic conditions the 
above configurational changes are superimposed onto the configurational 
changes due to the prescribed sinusoidal strain. ‘The following results are 


obtained, 
4 w?7? é 4 WT 
E’ = 3NkT > ——~., E” = 3NkT P__ 4 
» l + w*r’, 2X 1 + w?r, (4) 
and 
, w?7? E WT 
K’ =K 2, . 2 5 
: 2» 1 + w?r*, . X 1 + w?7?, (5) 


where 7, is the relaxation time characterizing motions of the p’th normal 
mode, N is the number of molecules per unit volume, k is Boltzmann’s 
constant, and 7’ the absolute temperature. K is an optical constant 
proportional to N and a — a,, the polarizability anisotropy of the sta- 
tistical link. To a good approximation the above summations may be 
replaced by integrals if the three longest relaxation times (71 —~ 73) are not 
involved. The theory then predicts that ¢,(In 7) and ¢x(In 7) should be 
proportional to 7~'”*. 

The theory also predicts that the birefringence should be in phase with 
the stress (6 — 6’ = 0) and that |C| should be independent of frequency. It 
follows that C’ should also be independent of frequency and that C” and 
¢@c(In 7) should be zero at all frequencies or relaxation times. These 
predictions arise from the fact that at any instant both the birefringence 
and the stress are determined by the internal configurational states of the 
molecules. The stress is related to these configurations via the con- 
figurational entropy. 

From eqs. (4) and (5) the limiting high frequency values of EF’ and K’ 
are given by FL’. = 3vNkT and K’.. = Kv, respectively. However, 
there is an upper frequency limit at which the theory is valid since con- 
figurational changes within the submolecules are ignored. Furthermore, 
at very high frequencies (or low temperatures) the elastic deformation 
mechanism will involve the distortion of valence angles and separation of 
chains against their forces or attraction. These local distortions will 
result in an internal energy change rather than the entropy change assumed 
in the above orientational mechanism. The release of these distortions 
should also contribute to the relaxation of birefringence, although at 
present no molecular theory has been developed for this mechanism. 
However the dynamic birefringence technique should provide information 
concerning the frequency ranges in which the orientational and distor- 
tional mechanisms, respectively, predominate. 


Experimental 


Our experimental method for measuring Ano, fo, €, 5, and 6’ over a range 
of frequency and temperature has been described previously.' In the 
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present work measurements have been made upon amorphous poly- 
acetaldehyde in the frequency range 1-20 cycles/sec. and temperature 
range 0 to —25°C. The polymer was prepared from acetaldehyde 
monomer at —75°C., using silica gel as catalyst, as described by Weissermel 
and Schmieder.‘ The crude polymer was twice precipitated from acetone 
solution with water, dried under vacuum, dissolved in dry benzene contain- 
ing a little Nonox antioxidant, and then freeze-dried. The molecular 
weight of the polymer, determined from intrinsic viscosity measurements in 
methyl ethyl ketone, was about 1.5 & 10°. The glass transition tempera- 
ture of the polymer, determined dilatometrically, was —30.4 + 1°C. For 
the dynamic birefringence measurements strips of the polymer, approxi- 
mately 2.5” & 0.5” X& 0.1”, were molded at 50°C. 


Results 


Mechanical 


The variation of £’, EH”, and tan 6 with frequency is shown in Figure 1. 
The different temperatures studied are indicated by different symbols as 
illustrated on the E’ plots. This same system of symbols is used in all 
later diagrams. The data shown in Figure 1 were obtained on three dif- 
ferent samples each from the same preparation. In the temperature 
range 0 to — 7.2°C. a +0.75% dynamic strain superimposed on a 16% 
static strain was employed. At lower temperatures a combination of 
+0.17% dynamic strain and 22% static strain was used. L’ is seen to 
increase with increasing frequency at each temperature, while E” exhibits 
a maximum at about 3-4 cycles/sec. at —21.5°C. A maximum in the tan 6 
curves is observed within the experimental frequency range at —10.8, 
— 13.1, and —15°C. 

Reduced master plots of log E’r,, log F'"7,, tan 6, log D’7,, and log D”r, 
against log azf are shown in Figure 2. The reference temperature 7) was 
conveniently chosen as 0°C., the highest temperature of measurement. 
At —18.7°C., and lower temperatures, the shift factors (log ar) derived 
from the formation of the modulus and compliance master curves are 
somewhat larger than those required for the tan 6 master plot. This 
discrepancy has probably arisen from errors in determining the stress 
magnitudes, since both the sample and proof ring were changed at this 


TABLE I 





Temperature, °C. Log ar Temperature, °C. Log ar 
, 


G.0 0 > .6 3.09 
—3.0 AS : 3.45 
—7.! wy 20. 4.01 
—10. Jf 21. 4.36 
—13. 2.37 25. 5.56 


—15. 
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Fig. 1. The variation of £’, FE”, and tan 6 with frequency for polyacetaldehyde. The 
different temperatures studied are indicated by different point symbols as shown on the 
E’ plots. 


temperature. Since the tan 6 values are independent of the stress mag- 
nitude, the shift factors obtained from the tan 6 plots are probably more 
accurate and furthermore they yield a smoother log a7 versus temperature 
curve. Hence the log ay values quoted in Table I are those derived 
from the tan 6 results. 
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Fig. 2. Reduced master curves of E'r,, E"r,, tan 6, D’r, and D"r, against log arf. 


The reference temperature 79 = 0°C. Point symbols as in Figure 1. 


The master curves given in Figure 2 are typical of those observed for the 
rubberglass relaxation of amorphous polymers. The maximum in the 
k” master curve occurs at a higher reduced frequency than the tan 6 
maximum and corresponds to the region where /’ is approaching its limit- 
ing glassy value. A maximum in D” is not observed within the range of 
reduced frequency investigated, but it is estimated to occur at about | 
cycles/sec. at 0°C. 
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Birefringence 


Figure 3 illustrates the determination of the factors Anz,/Anz which are 
required inthe proposed method for reducing the optical data to the 
reference temperature 7. The top diagram shows the variation of 
birefringence with temperature (in the range 0 to —30°C.) for a polyacet- 
aldehyde sample maintained at 19.7% extension. The sample was initially 
stretched at 0°C. and the birefringence allowed to relax to a steady value. 
The agreement between the results obtained during the cooling and 
subsequent heating experiment shows that the measurements were re- 
versible. The small increase in birefringence with decreasing temperature 





| 








fh ccintinsciianienintsi emesis 
-20 -10 
Temperature (c) 
Fig. 3. (Upper), plot of the birefringence against temperature for polyacetaldehyde 
at 19.7% extension during cooling (©) and reheating (A); (lower), the variation of Ano°c./ 
Anr with temperature for polyacetaldehyde. 


is attributed to an increase in ay — a,, the equilibrium optical anisotropy 
of the statistical link. The values of Anoec./Anr shown in the lower dia- 
gram are averages determined from two separate runs. 

The master plots of |K|r,, Cle 6’, and 6 — 6’ against log a7f are given 
in Figure 4. At a reduced frequency of about 10‘ cycles/sec. |Klz, exhibits 
a maximum and the birefringence and strain are in phase (6’ = 0). At lower 
frequencies |K|r, increases with frequency and the birefringence cycle 
leads the strain cycle (6’ positive), whereas at higher frequencies the 
reverse situation exists. |C|, is seen to increase as the frequency de- 
creases but tends toward the static value of 1.1 K 10~—" em.?/dyne at the 
lowest reduced frequencies. This static stress-optical coefficient, indicated 
by the point of the ordinate, was measured for a (rubberlike) specimen at 
room temperature after hanging known weights onto its end. The stress 
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Fig. 4. Master curves of |K|7,, !Cl7,, 6’, and 6 — 6’. Point symbols as in Figure 1. 
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Fig. 5. Master curves of K'r, and K",7,. Point symbols as in Figure 1. 
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Fig. 6. Master curves of C 7, and C"7,. Point symbols as in Figure 1. 


cycle leads the birefringence cycle at all reduced frequencies (6 — 6’ posi- 
tive), and the 6 — 6’ plot passes through a broad maximum in the region 
of 10° cycles /sec. 

Plots of K’;, and K";, against log azf are shown in Figure 5. At low 
frequencies the increase in K’;, with increasing frequency is accompanied 
by a positive K";, peak. In the high frequency region a decrease in K’,, 
with increasing frequency is observed and also a corresponding negative 
K" 7, peak. 

Figure 6 shows the C’;, and C";, curves. The C",7, curve has a single 
maximum which correlates approximately with the inflection point in the 
C7, plot. 

Plots of log ¢z(In 7), log @p(In 7), dx(In 7) and ¢¢(In 7) against log r at 
7’) = 0°C. are shown in Figure 7. These distributions were evaluated from 
the £"7,, D"7,, K"r, and C"7, master curves respectively according to the 
first-order approximation equations given. The dashed line shown in the 
long 7 region of the log ¢z(In 7) versus log 7 plot has a slope of —'/2 as 
predicted by the molecular theory. The slope of the experimental curve 
generally differs from —'/, but tends to this value at the longest 7 inves- 
tigated. A maximum in the ¢g¢ spectrum is observed at r ~ 107° sec. 
while the @p spectrum appears to tend toward a maximum in the long 7 
region. The ¢x spectrum divides into two regions. Processes for which 
r < 10~*? sec. contribute a net negative amount to K’.. — Ky and ¢x 
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Fig. 7. Plots of log ¢,(In 7), log ¢)(In 7), ¢x(In 7), and @-(In 7) against log rat O°C. The 


dashed line shown with the ¢, spectrum has the theoretical slope of —'/2. 


exhibits a minimum at 7 ~ 10~** sec. which correlates with the me- 
chanical ¢, maximum. For 7 > 10~‘*-’sec. the relaxation processes contrib- 
ute a net positive amount to K’.. — K’y and ¢x shows a maximum at 
rt ~ 10-*8sec. The dc spectrum has a single maximum at 7 ~ 10~**’ sec. 

The apparent frequencies of the E”, K”, and C” maxima and the K” 
minimum are plotted against 1/7’ in Figure 8. Points are included on 
the diagram for which a maximum was directly observed within the 
experimental frequency range. Other points were subsequently derived 
from the log a7 values given in Table I. The dashed curve illustrates the 
estimated frequency-temperature location of the D” maximum assuming that 
at 0°C. the peak occurs at 1 cycle/sec. Also plotted in l‘igure 8 are points 
corresponding to the maximum in the dielectric loss factor ¢«”.6 These 
points lie fairly close to those corresponding to the K” maximum. 


Discussion 


The observed frequency dependence of the mechanical and optical 
properties probably arises from the motions of segments within the main 
polymer chain and corresponds to the transformation from the rubberlike 
to glasslike state of polyacetaldehyde. 

The high modulus observed on the high frequency side of the reduced 
modulus curve is typical of values found for polymers in the glassy state 
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and may be attributed to the distertion of local groups within the structure. 
The corresponding distortional birefringence in this region leads to rel- 
atively small positive values for K’;, and C’;,. The initial increase in 
K’,, with decreasing frequency indicates that motions arising from the 
release of these local distortions produce an increase in the average ori- 
entation of chain links. Such a mechanism would be responsible for the 
negative peak in the K”,7, master plot and the maximum in the "7, 
curve and also the associated peaks in the ¢x and ¢x spectra at + ~ 10-58 
sec. at 0°C. 

After passing through a maximum, the subsequent decrease in K’,, 
with decreasing frequency suggests that an additional relaxation mech- 
anism predominates in the low frequency region. It is probable that this 
mechanism involves the motions of relatively long flexible chain segments 
such that at constant strain the chains tend to more coiled-up configura- 
tions as the frequency decreases. This process, which forms the basis of 
the molecular theory outlined earlier, leads to the positive peaks in the 
K",, and C” 7, plots and the corresponding maxima in the ¢x and @¢ spec- 
tra. However the theoretical predictions are not obeyed since the slope of 
the log ¢z(In 7) versus log 7 plot differs from —'/2, C’7, is not independent 
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of frequency, and 6 — 6’, Cr, and ¢¢ are not zero. A probable explanation 
for the failure of these predictions is that the ‘distortional”’ mechanisms, 
which predominate in the high frequency region, overlap the orientational 
processes which predominate at low frequencies. In order to resolve 
satisfactorily the two proposed mechanisms a molecular theory appropriate 
to the high frequency range is required. The correlation noted above 
between the frequency-temperature location of the K” maximum and the 
dielectric «” peak indicates perhaps that the dielectric relaxation in poly- 
acetaldehyde is determined largely by the longer range chain rearrange- 
ments. 

In the case of PMA it was also concluded that overlapping distortional 
and orientational mechanisms were involved in the relaxation region as- 
sociated with the rubber-glass transformation. However in this case 
chain orientation produces a negative birefringence, possibly owing to the 
orientation of —COOCH,; groups perpendicular to the stretching direction. 
Consequently a positive K” peak was observed at high frequencies and a 
negative K” peak in the low-frequency range. Also, for PMA a close 
correlation was found between the locations of the high frequency K” 
maximum and the e” maximum. Hence, unlike polyacetaldehyde, the 
dielectric relaxation of PMA probably results from the short-range dis- 
tortional rearrangements which could perhaps induce, or modify, rotations 
of the dipolar side groups. 

The author thanks Mrs. J. Hawkins for preparing the polyacetaldehyde and for general 
assistance. The work described above has been carried out as part of the research pro- 
gram of the National Physical Laboratory and this paper is published by permission of the 
Director of the Laboratory. 
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Résumé 


On souligne les principes 4 la base de la méthode de biréfringence dynamique en 
attribuant une emphase particuliére A l’usage de cette technique dans la recherche de 
mécanismes moléculaires de relaxation dans les polyméres amorphes. Les résultats sont 
donnés pour le polyacétaldéhyde amorphe ([—CH—O—],) pour des fréquences de 


CH; 
l’ordre de 1 4 20 cycles/seconde et des températures de l’order de 0°C 4 —25°C. 
Ces résultats sont discutés et comparés bri¢vement avec les données de biréfringence 
dynamique obtenues précédemment pour le polyméthacrylate de méthyle. 
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Zusammenfassung 


Die Prinzipien fiir die Methode der dynamischen Doppelbrechung werden unter be- 
sonderer Betonung der Verwendung dieser Methode zur Untersuchung des molekularen 
Relaxationsmechanismus bei amorphen Polymeren dargelegt. Daten fiir amorphen 
Polyacetaldehyd ({[—CH—O—}],,) werden im Frequenzbereich von 1-20 Hertz und im 

CH: 
Temperaturbereich von 0 aus —25°C angegeben. Die Ergebnisse werden diskutiert und 
kurz mit dielektrischen Daten fiir Polyacetaldehyd und auch mit friiheren dynamisclien 
Doppelbrechungsdaten fiir Polymethylacrylat verglichen. 
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Temperature Dependence of Orientation Birefringence 
of Polymers in the Glassy and Rubbery States 


R. D. ANDREWS* and T. J. HAMMACK, * Materials Research Laboratory, 
Department of Civil Engineering, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 


Synopsis 


The temperature dependence of orientation birefringence has been measured for 
several different polymers, by producing samples with a fixed degree of orientation by 
mechanical stretching, and varying the temperature. The following three polymers 
were studied: polystyrene, poly(methy] acrylate), and poly(methyl methacrylate )—the 
last in both atactic and stereospecific forms. These all have negative orientation bire- 
fringence at room temperature. Measurements were carried out in a range of tempera- 
tures below the glass transition temperature, except for poly(methy! acrylate) which was 
crosslinked. The curves obtained for the different polymers showed considerable dif- 
ferences in form; they were reversible in regard to temperature cycling, and showed no 
appreciable time effects. Polystyrene showed little change in its orientation birefrin- 
gence with temperature. The other polymers showed appreciable changes with tem- 
perature, and an actual change of sign of the birefringence from negative to positive with 
increasing temperature. The crossover temperature was near the glass transition 
temperature for the various stereo forms of PMMA, and measurement of this crossover 
temperature might therefore provide a fairly accurate method for determining the stereo 
composition of a PMMA sample. However, the crossover temperature is considerably 
above 7’, in the case of PMA. We believe that these curves reflect thermal motions of 
the molecular groups with polarizability anisotropy (such as the phenyl] and ester side 
groups) whose orientation gives rise to the birefringence. The small change in the case 
of polystyrene is undoubtedly a consequence of the rigidity of the phenyl ring and the 
steric hindrance of its immediate environment. These results may also reflect chain 
backbone motions, and other effects, in an indirect way. The theory and experimental 
data for a complete interpretation of these results is not yet available. However, some 
of the evident consequences of the present findings are discussed. 


Birefringence studies of solid polymers have generally involved measure- 
ments of the stress-optical coefficient, the strain-optical coefficient, or ori- 
entation birefringence at some standard temperature (usually room tem- 
perature) as an index of the degree of orientation. In the present paper 
we describe some preliminary results from a somewhat different type of 
experiment, in which samples with a fixed degree of orientation are pro- 
duced, and the effect of temperature alone on the orientation birefringence 
is measured. This study of the temperature variable has provided results 
of unusual interest. 


* Present address: Textile Division, Department of Mechanical Engineering, M.I.T. 
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The polymers examined in the present study were: 


(1) Polystyrene—linear 
atactic 

(2) Poly(methy] acrylate)—crosslinked 
atactic 

‘“. (3) Poly(methyl methacrylate)—linear 

atactic 
isotactic 
syndiotactic 
stereo-block 


These polymers were amorphous in all cases; the stereospecific PMMA 
samples were quenched into the amorphous glassy state in the original 
compression molding process. The three atactic materials were polymerized 
in place as cast sheets. All of these polymers have negative orientation 
birefringence at room temperature. 


Samples and Experimental Procedure 


The original polymer sheets were 1/i—'/s in. thick and a few inches 
square; small strips were cut from these (by machining, or with a razor 
blade in the case of PMA) which were perhaps */s in. wide and 3 in. long. 
These strips of the glassy polymers were suspended in a temperature box 
under moderate load, and: the temperature was slowly raised to a point 
slightly above the glass transition temperature. The sample was allowed 
to creep to a suitable extent, to give the degree of orientation desired, and 
the box was then opened and the sample quickly cooled by circulation of 
room-temperature air, while remaining under load. The sample could then 
be removed if desired. The procedure was slightly different for the poly- 
(methy] acrylate), which was in the rubbery state at room temperature and 
apparently crosslinked; in this case the sample was loaded and allowed 
to reach elastic equilibrium in the rubbery state, and this fixed load was 
then kept on the sample at all times as the temperature was varied. 

Temperature changes involved the use of two different air thermostats, 
one operating above room temperature with electric heating coils of the 
conventional sort, and the other operating below room temperature, down 
to liquid nitrogen temperature, in which liquid nitrogen was added to a 
double-cylinder brass cooling jacket surrounding the sample. The con- 
struction of the latter thermostat will be described elsewhere. The un- 
crosslinked polymers could not be heated above the glass transition tem- 
perature in the oriented state, since this would have led to significant dis- 
orientation. Birefringence was measured by use of a Sénarmont com- 
pensator, which has been described in a previvus publication.! Birefringence 
was not calculated as such; the data are presented simply as compensator 
extinction angle values (which are proportional to the birefringence) as a 
function of temperature. The aim in the present measurements was 
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merely to establish the general form of the curves obtained in this type of 
experiment. 


Experimental Results 


Polystyrene 


The results obtained for a sample of cast polystyrene sheet are shown in 
Figure 1. The temperature was varied between liquid nitrogen tempera- 
ture (—196°C.) and +50°C., and the curve is essentially flat, showing an 
almost complete lack of temperature dependence of the orientation bire- 
fringence. Actually there appears to be a slight minimum in the curve, 
at about 0°C., which might justify some further investigation. 


‘ 
g 
3S 
5 


-100 -50 0 100 ISO 200 
TEMPERATURE , °C 


Fig. 1. Temperature dependence of orientation birefringence: polystyrene cast sheet. 


Poly(methyl acrylate) 


An analogous curve for poly(methy! acrylate) is shown in Figure 2. 


This sample appeared to have been crosslinked (presumably by chain trans- 
fer) during the polymerization, and therefore the temperature range 
covered could be extended into the rubbery as well as the glassy state, and the 
curve shown extends from liquid nitrogen up to about +115°C., approxi- 
mately 110° above the glass transition temperature of +4°C. As noted 
previously, these measurements were carried out by hanging a dead load 
on the sample in the rubbery state and changing temperature under this 
constant load. The experiment therefore corresponded strictly neither to 
constant length nor constant elongation, because of combined effects of 
thermal expansion and the temperature dependence of modulus in the rub- 
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Fig. 2. Temperature dependence of orientation birefringence: poly(methyl acrylate). 
bery state (kinetic theory of rubber elasticity); however, the experiment 
‘an still be regarded as being carried out at an approximately constant de- 
gree of orientation, since these perturbing effects are small. Carrying out 
the experiment by holding the sample at a fixed length by rigid clamping 
would actually be less satisfactory: Stein, Krimm, and Tobolsky? carried 
out a measurement in this way on o sample of crosslinked poly(ethy] acry- 
iate), and they found that when the sample was cooled below the glass 
transition, thermal contraction introduced a distortional stress in the 
sample; this introduced a corresponding distortional birefringence effect 
of opposite sign, which caused the experimental curve to turn around and 
increase at the lowest temperatures. | 

As seen from the figure, there is ¢ relatively small change of orientation 
birefringence with temperature in the glassy state, although the effect is 
larger than was observed with polystyrene: the negative orientation bire- 
fringence decreases in magnitude with increasing temperature, at a slightly 
increasing rate as the glass transition is approached. As temperature in- 
creases above 7’, (which is indicated in the figure by a vertical dashed 
line) there is an abrupt rise in the curve, followed by a change in slope in the 
neighborhood of +40°, leading to a perfectly linear region of the curve. 
An extrapolation of the final linear region suggests that the birefringence 
will cross the zero axis at approximately +130°C. and become positive in 
sign. The sample unfortunately broke before this temperature could be 
reached. It is of interest that Stein, Krimm, and Tobolsky observe no 
abrupt change around the glass transition in their birefringence-tempera- 
ture curve for crosslinked poly (ethyl! acrylate). 
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Poly(methyl methacryiate) 


Somewhat more detailed measurements were made on PMMA, including 
investigation of the effects of stereospecific structure. Although the type 
of experiment which we are carrying out here seems to be a relatively new 
one, some fragmentary data of this sort on PMMA have previously been 
reported. Kolsky and Shearman,’ for example, have published a curve 
showing the effects of temperature cycling on the birefringence of a sample 
of oriented PMMA, in the glassy state. A _birefringence-temperature 
curve has also been published by Trapeznikova and Zhurina‘ which is 
reproduced as Figure 3. They made their measurements on an oriented 
sample prepared from a solvent-cast film (presumably atactic polymer), 
and they report their data in the same manner as our own—i.e., as Sénar- 
mont compensator extinction angle vs. temperature. This curve is not in 
very good agreement with our own curves for PMMA, as will be pointed 
out later. Other data on Plexiglas II commercial cast sheet have been 
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Fig. 3. Temperature dependence of orientation birefringence: poly(methyl methacry- 
late) data of Trapeznikova and Zhurina.‘ 


reported in a thesis by McLoughlin.’ He oriented strip samples by load- 
ing at temperatures above the glass transition (7’, = 100°C., roughly), 
followed by cooling through the glass transition region under dead load. 
He observed a positive birefringence on initial loading, and the curves of 
birefringence vs. temperature passed smoothly through zero to increasingly 
negative values as the temperature progressively decreased. His data 
covered only a limited temperature range, however. 

Data of our own on PMMA are shown in Figure 4, which allow a compari- 
son of curves for atactic, isotactic and syndiotactic polymers. These dif- 
ferent polymers all give curves of the same general form. The curves show 
an appreciable change of orientation birefringence with temperature in the 
glassy state, and even at liquid nitrogen temperature the curves have not 
flattened to a final constant value. This argues for the existence of a sig- 
nificant degree of molecular motion in the glassy state, in contrast to the 
concept of a glassy polymer as a totally frozen structure; this will be dis- 
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Fig. 4. Temperature dependence of orientation birefringence: poly(methyl methacry- 
late), different stereo forms. 


cussed in more detail in a later section. As temperature increases, the 
curves rise at an increasing rate and finally cross the zero axis, where the 
birefringence shows a change of sign from negative to positive, in agreement 
with the observations of McLoughlin. It will be noted that Trapeznikova 
and Zhurina draw their curve so as to approach the zero axis asymptotic- 
ally as temperature increases, although their data points would equally well 
allow the curve to be drawn through the zero axis into the positive region. 
If the curve were drawn in this way, their crossover point would be about 
+50°C., which would be considerably below that (+110°) which we ob- 
serve for our own atactic sample; this is probably due to residual solvent 
in their film which would lower the glass transition temperature of their 
polymer. Finally, we see no evidence in our own curves of the abrupt 
drop which they show at the low-temperature end of their curve, at about 
—100°C.; this must be some artifact in their measurements. 

The crossover points (or temperatures) for the different stereo forms are 
indicated in Figure 4 by vertical dashed lines above the zero axis; these 
values seem to correspond roughly with the glass transition temperatures 
determined for these same polymers from volume-temperature measure- 
ments, as follows: 


Polymer Crossover temp. (°C.) (°C. 





Isotactic +31 +48 
Atactic +110 +100 
Syndiotactic +135 +110 
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However, there is no reason to think that the crossover temperature and 
glass transition temperature should coincide, and the lack of generality of 
any such principle, even for acrylic polymers, is evident from the curve for 
PMA in Figure 2, where the crossover point lies more than 100° above 
T,. Two curves are shown for both the isotactic and atactic polymers 
in Figure 4; these were obtained from samples with different degrees of 
orientation, and they show that such curves differ essentially only by a ver- 
tical scale factor, and seem to have the same crossover point (ef. particu- 
larly the isotactic curves). This suggests that the crossover temperature, 
which can be easily measured, may provide a convenient method of deter- 
mining the average stereo composition of a polymer suck as PMMA with 
good accuracy. It is interesting that the crossover temperature for our 
atactic polymer (which is believed from other measurements to contain about 
75% syndiotactic structure) lies about 75% of the way between the values 
for the isotactic and syndiotactic polymers; this would suggest that the 
relation between crossover temperature and stereo composition might be 
almost linear. 

The birefringence vs. temperature curve for a stereo-block PMMA (which 
is a block copolymer of isotactic and syndiotactic sequences, having a dila- 
tometric 7’, value of +82°C.) is shown in Figure 5. This curve is similar 
in form to the other PMMA samples, and appears to have a crossover tem- 
perature of about +97°C. However, the curve seems to show a peculiar 
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Fig. 5. Temperature dependence of orientation birefringence: poly(methyl methacry- 
late), stereo-block. 
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. 


inflection at about. —50°C. which was not seen in any of the curves in Fig- 
ure 4; this requires further verification, but is presented here as being of 
general interest, and because it may have some relation to the stereo- 
block structure. 

Experiments with temperature cycling indicated that the curves ob- 
tained in all these experiments were reversible: the same curve was ob- 
tained with both increasing and decreasing temperature. Also, no time- 
delay effects were observed: birefringence changes followed temperature 
changes instantly, as nearly as could be determined. 


Discussion 


In interpreting these results, as in all cases of birefringence phenomena 
in polymers, we may postulate the following general principle: Birefrin- 
gence arises from the orientation of molecular groups having polarizability ani- 
sotropy. Some degree of polarizability anisotropy is associated with the hy- 
drocarbon backbone chain of the polymers investigated here. However, the 
anisotropy is predominantly associated with the phenyl and ester side groups. 
This is illustrated in a schematic way in Figure 6, where longer and shorter 
arrows have been drawn beside these groups to indicate the approximate 
magnitudes and directions of their principal polarizabilities. 


SIDE GROUP POLARIZABILITIES ROTATIONS 
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Fig. 6. Principal polarizabilities and rotational motions for phenyl and ester side groups. 


The reversibility of the birefringence-temperature curves in regard to 
temperature cycling, and the effectively instantaneous ‘change of birefring- 
ence with temperature change (no time-delay effects being observable) 
strongly suggests that the birefringence values which we are observing are 
thermal equilibrium values, which are associated with the thermal energy 
level of the material. More specifically, it would appear that we are ob- 
serving the effects of thermal agitation (vibration and rotation) on the state 
of orientation of the side groups which were produced by extending the chain 
backbone during the original orientation process. 

In the case of polystyrene, the polarizability of the phenyl] ring is approxi- 
mately twice as great in the plane of the ring as in the direction perpendicu- 
lar to the plane, and the ring is a rigid, planar structure which has no inter- 
nal degrees of freedom of motion; the only motion possible is a rotation 
around the bond connecting it to the chain backbone, as illustrated in Fig- 
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ure 6. The lack of temperature dependence of the birefringence in oriented 
polystyrene seems to illustrate not only the rigidity of the ring, which is 
unaffected by temperature change, but also a lack of rotation of the ring 
about the connecting bond to the chain. This suggests that rotational 
motion of the ring must be highly hindered sterically; and previous in- 
vestigations of the stress-optical effect in glassy polystyrene® have indicated 
that this is indeed the case. Signer’ has also pointed out that for poly- 
styrene to have a negative orientation birefringence, the plane of the phenyl! 
ring must lie not in the plane defined by the oriented chain backbone (which 
would produce an approximately zero birefringence contribution), but in 
a plane roughly perpendicular to it. 

The greater temperature dependence of the orientation birefringence in 
the acrylic polymers seems to reflect the greater mobility of the nonrigid 
ester side group. The rotation of the end methyl! group around the ad- 
jacent C—O bond axis (see Fig. 6) can have an effect not only on the ani- 
sotropy of the ester group, but may also facilitate rotation of the ester group 
as a whole around the bond connecting it to the chain backbone. In anal- 
ogy with the behavior of the phenyl ring, the —COO— part of the ester 
group is a rigid planar structure, and if this plane were to lie in the plane 
of the oriented chain backbone, the birefringence effect would be positive; 
this would become negative as the plane rotated into a position perpendic- 
ular to the backbone plane. It would therefore be possible to explain the 
decrease of negative birefringence with temperature, and even the transi- 
tion to positive values, on the basis of thermal motions of the ester group 
alone. This is essentially the explanation offered by Kolsky and Shearman® 
to explain their own data, and we are in agreement with this interpretation. 
Tsvetkov and Verkhotina® have estimated, using chain models and bond 
polarizabilities, that to obtain a positive birefringence for oriented PMMA 
requires a practically free rotation of the ester group. A somewhat dif- 
ferent explanation is offered by Trapeznikova and Zhurina.‘ They pro- 
pose that the observed birefringence in PMMA is a summation of two con- 
tributions—a negative contribution from the ester group, and a positive 
contribution from the alpha-methyl group on the chain. They interpret 
the positive birefringence observed at high temperatures as resulting from 
the disorientation of the negative component; they note that the positive 
birefringence could not be “frozen in” even with very rapid cooling. Peu- 
kert® also was surprised to find that it was not possible to freeze in the posi- 
tive birefringence first produced when PMMA was hot-stretched. He 
gives a slightly different version of the two-contribution theory, in which 
he postulates a positive component due to molecular orientation and a 
negative component due to stress. 

Assigning a significant positive contribution to the alpha-methyl group 
seems unreasonable, since the polarizability anisotropy of this group is 
almost zero; its major effect is probably simply as a source of steric hin- 
drance to motion of the ester group. The effect of the chain backbone re- 
quires more extended comment. The polymers studied here have the 
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major part of their polarizability anisotropy concentrated in the side groups, 
and although the oriented chain backbone should contribute a small posi- 
tive effect, this is less important than the effect of the side groups. It 
would be interesting to carry out this same type of measurement on poly- 
mers (such as natural rubber and silicone) where the predominant contribu- 
tion comes from the chain backbone. In addition, it should be noted that 
backbone motions will of course affect the orientation of side groups in- 
directly, and this may. be making a contribution to the form of our present 
curves. The abrupt change in birefringence observed in PMA just above 
the glass transition (Fig. 2) might be explainable as resulting from the sud- 
denly increased freedom of motion of the chain backbone. 

Such effects as a change in equilibrium chain configuration with tempera- 
ture (due to a change in statistical distribution between trans and gauche 
forms) would also be expected to show up in measurements of this sort, 
though this type of effect should be observable particularly in the rubbery 
state above 7',, rather than in the glassy state. Stein, Krimm, and To- 
bolsky? interpret the change in birefringence with temperature which they 
observe in crosslinked poly (ethyl acrylate) simply as a change in size of the 
“statistical segment” with temperature. This would be one way of look- 
ing at the effects of chain configurational energy. 

Data of the present type clearly have important implications in connec- 
tion with the interpretation of the variation of stress-optical and strain- 
optical coefficients with temperature. And in a broader sense, if birefrin- 
gence behavior is studied as a function of both temperature and time (or 
frequency), in a way analogous to static and dynamic rheological or vis- 
coelastic properties of polymers (as has been recommended in another pub- 
lication”), it seems evident that time-temperature superposition proce- 
dures, for example, cannot be successfully carried out for birefringence data 
without an explicit allowance being made for this intrinsic temperature 
dependence of birefringence, in addition to the effect of temperature on the 
rates of the rheo-optical processes. 

The interpretation of these curves as reflecting thermally activated mo- 
lecular motions suggests that it would be valuable to compare data of this 
sort with other types of physical data which also give indications of mo- 
lecular motions, such as that from dynamic mechanical, dielectric, and nu- 
clear magnetic resonance measurements. An interesting and also slightly 
puzzling point, for example, is that the birefringence-temperature curves 
for PMMA show no evidence of the 6-peak or dispersion region below the 
glass transition which is seen clearly in dynamic mechanical and dielectric 
data, and‘is interpreted in these latter cases as resulting from the rotation 
of the ester side group. 

Still another valuable use for data of this type would be in providing a 
clearer understanding of the glass transition phenomenon. The glass 
transition is interpreted as a change in the state of motion of the mole- 
cules, and these orientation birefringence curves measured through the 
glass transition region, when properly interpreted, should provide consid- 
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erable information on the exact nature of the changes in molecular motions 
which take place. 

In conclusion, it is clear that these birefringence-temperature curves 
show a great deal of individual character in the case of different polymers, 
and it therefore seems likely that considerable information of value could 
be deduced from them. Unfortunately, the theory needed for this pur- 
pose is not available at the present time, and must be developed. It is 
planned to obtain additional data of this type covering a wider range of 
polymers and experimental conditions, to provide clearer indications of the 
general nature of the theory which is needed, and thus facilitate its rapid 
and successful development. 
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Résumé 


La dépendance vis-a-vis de la température de la biréfringence d’orientation a été 
mesurée pour différents polyméres, en faisant des échantillons possédant un degré 
d’orientation déterminé par étirage mécanique et variation de la température. On a 
étudié les trois polyméres suivants: le polystyréne, le polyacrylate de méthyle, et le 
polyméthacrylate de méthyle—ce dernier dans ses formes atactiques et stériospécifiques. 
Tous ces polyméres ont une biréfrirgence d’orientation négative 4 température ambiante. 
Des mesures ont été effectuées dans un domaine de température en dessous de la tempéra- 
ture de transition vitreuse, excepté pour le polyacrylate de méthyle, qui était ponté. 
Les courbes obtenues pour les différents polyméres montrent des différences considérables 
de forme: elles sont réversibles quant 4 la variation cyclique de la température et ne 
montrent pas d’effets de durée appréciables. Le polystyrene montre un léger change- 
ment de biréfringerce d’orientation avec la température. Les autres polyméres mon- 
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trent des changements appréciables avec la températures et un changement de signe de la 
biréfringence, du négatif au positif par une augmentation de la température. La tem- 
pérature de changement de signe est proche de la température de transition vitreuse 
pour les différentes stéréoformes du PM MA, et la mesure de cette température de change- 
ment de signe peut donc constituer une méthode excellente pour déterminer la composi- 
tion stéréo des échantillons de PMMA. Néanmoins la température de changement de 
signe est considérablement plus élevée que 7’, dans le cas du PMMA. Nous croyons que 
ces courbes reflétent les mouvements thermiques des groupements moléculaires dus 4 
l’anisotropie de polarisabilité (comme les groupes latéraux phényles et esters), et dont 
l’orientation engendre une augmentation de la biréfringence. Le faible changement 
dans le cas du polystyrene, est indubitablement la conséquence de la rigidité du cycle 
phényle, et de l’empéchement stérique de son environnement immédiat. Ces résultats 
peuvent ainsi refléter des mouvements de la chaine principale et d'autres effets, d’une 
La théorie et les données expérimentales ne permettent pas encore de 


fagon indirecte. 
Néanmoins certaines conséquences évidentes des 


donner une interpretation complete. 
résultats actuels sont discutés. 


Zusammenfassung 


Die Temperaturabhangigkeit der Orientierungsdoppelbrechung wurde fiir verschiedene 
Polymere durch Herstellung von Proben mit einem fixen Orientierungsgrad durch 
mechenische Dehnung und Variierung der Temperatur bestimmt. Folgende drei Poly- 
mere wurden untersucht: Polystyrol, Polymethylacrylat und Polymethylmethacrylat— 
letzteres sowohl in ataktischer als auch stereospezifischer Form. Alle Proben zeigen bei 
Raumtemperatur negative Orientierungsdoppelbrechung. Die Messungen wurden in 
einem Temperaturbereich unterhalb der Glasumwandlungstemperatur aus gefiihrt, 
ausgenommen Polymethylacrylat, das vernetzt war. Die fiir die verschiedenen Poly- 
meren erhaltenen Kurven zeigen betriichtliche Formunterschiede, sie waren in bezug auf 
den Temperaturzyklus umkehrbar und zeigten keine merkliche Zeitabhingigkeit. Poly- 
styrol zeigte kaum eine Anderung seiner Orientierungsdoppelbrechung mit der Tem- 
peratur. Die anderen Polymeren zeigten merkliche Anderungen mit der Temperatur, 
und sogar einen Vorzeichenwechsel der Doppelbrechung von negativ zu positiv mit 
zunehmender Temperatur. Die Ubergangstemperatur lag fiir die verschiedenen Stereo- 
formen von PMMA in der Nihe der Glasumwandlungstemperatur; Messungen dieser 
Ubergangstemperatur kiénnten daher eine ziemlich genaue Methode zur Bestimmung der 
Stereozusammensetzung einer PMMA-Probe bilden. Bei PMA liegt jedoch die Uber- 
gangstemperatur betrachtlich iiber 7',. Wir glauben, dass diese Kurven die thermische 
Bewegung der Molekiilgruppen mit Polarisationsanisotropie (wie die Phenyl- und Ester- 
seitengruppen) wiedergeben, deren Orientierung die Doppelbrechung verursacht. Die 
kleine Anderung bei Polystyrol ist zweifellos eine Folge der Starrheit des Phenylringes 
und der sterischen Hinderung seiner unmittelbaren Umgebung. Diese Ergebnisse kénn- 
ten indirekt auch die Bewegung der Hauptkette und andere Effekte erkennen lassen. 
Die Theorie und die experimentellen Daten fiir ein vollstindige Erklirung dieser Ergeb- 
nisse sind zwar noch nicht greifbar, einige offensichtliche Folgerungen aus den vorlie- 


genden Ergebnissen kénnen aber diskutiert werden. 
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Synopsis 


Assuming that the crystalline polymer consists of the crystalline and the non-crystal- 
line phase, a mechanical model for the crystalline polymer is constructed. We introduced 
into this model such complicated circumstances as some parts of the crystalline phase will be 
largely deformed under a stress concentration and other parts not so much. The tempera- 
ture dependence of complex moduli of the crystalline polymers with various degrees of 
crystallinity is illustrated by using the model. On the basis of the model, we have derived 
an equation for calculating the strain-optical coefficient of the crystalline polymer from 
those of the crystalline and the non-crystalline phase. The equation gives a relation 
among the rheo-optical quantities, the viscoelastic quantities and the parameters repre- 
senting the fine structure and the crystallinity of the polymer. It is concluded from this 
equation that the strain optical coefficients cannot be expressed by such a simple additive 
relation as usually employed. 


1. Introduction 


One method of description of the rheo-optical properties of polymers has 
been recently presented by Stein, Onogi, and Keedy.! They introduced 
two distribution functions A(r) and B(r) to express the birefringence of 
polymers in the similar way as in the viscoelasticity. The two distribution 
functions for amorphous polymers have a possibility to be interpreted by 
the molecular theory. It is, however, difficult to interpret the distribution 
functions for the crystalline polymers by the molecular thecry, because 
there exists a heterogeneity of molecular aggregation in the crystalline tex- 
tures. Asa result of this circumstance, some parts in the crystalline region 
will be largely deformed under a stress concentration and other parts not 
so much. 

In this report, a model method developed by us to interpret the visco- 
elastic properties of the crystalline polymers and the polymer blends is 
applied to express the relationship among the rheo-optical properties, the 
viscoelastic properties and the texture of the crystallized polymer. In this 
case, the texture of the crystalline polymer is simplified to be expressed by 
the degree of crystallinity and one parameter representing the mixing state. 

113 
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2. Equivalent Model Method for Viscoelastic Properties of Two-phase 
System 


The left side of Figure 1 is a schematic representation of the two-phase 
system, in which the particles 2 (as phase 2) are dispersed in the medium 1 
(as phase 1). The phase 1 and the phase 2 are perfectly bonded to each 
other. 

When such a sample as the left side of Figure 1 is subjected to tension, 
some of the lines of force in the sample pass only through the phase 1 and 
the others through both the phases in series. 

Taking the above circumstance into consideration, the sample may be 
symbolized in the simplest case as shown in the right-hand side. A and ¢ 
in the figure are parameters relating to the mixing state and the composition 


ie 


Fig. 1. Schematic representations of the mixing state of the two-phase system (the left 
side) and their equivalent models (the right side). The upper figure shows a homogen- 
eously dispersed system and the lower figure a heterogeneously dispersed system. 


(b) 


Fig. 2. Two mechanical models for calculating the modulus of the equivalent model in 
Figure J. Model (a) is adopted.in our method. 
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Fig. 3. E’ and E” versus temperature for the interpolymer of polybutadiene (PBD) 
with styrene (ST) and acrylonitrile (AN) monomer at 138 cycles/sec.: (©) observed 
values for the volume fraction of PBD, V2 = 0.23; (@) observed values for the volume 
fraction of PBD, V2 = 0.44; solid lines for white circles are calculated based upon the 
equivalent model with \ = 0.600, ¢ = 0.380, and V2 = 0.23; solid lines for black circles 


are calculated based upon the equivalent model with \ = 0.775, ¢ = 0.568, and V2 = 
0.44; dotted lines are the observed ones for the component polymers. 


of the sample. When these simplifications are allowed for the actual two- 
phase system, the models on the right-hand side of Figure 1 can be said to 
be the equivalent models for them. ‘Two simple mechanical models are 
considered for calculating the modulus of the equivalent model in the right 
hand side. Figure 2 shows these two models. In the model (a) in Figure 
2, the strain of the element of the phase 2 is the same as that of the right- 
element neighboring it. In the model (6), the strain of the phase 2 is 
different from those of the other elements. It has been proved by the 
model experiments that the model (a) is more suitable to calculate the 
modulus of the equivalent model.? Therefore, the model (a) will be 
adopted for calculating the modulus of the actual system. 

The gross complex modulus /£* of the system, in which the particles 2 
are dispersed in the medium 1, is expressed by use of its equivalent model 
of Figure 1 as follows: 


4 = 


ii g Seog el 1 
E — r»)Ei + ABR t . (1) 


where EZ} and E3 are the complex moduli of the medium and the dispersed 
particles respectively, \ and ¢ are the parameters representing the mixing 
state, and they are in the relation that A¢ is equal to the volume fraction of 
the particles. When the particles are spheres, \ and g depend on the vol- 
ume fraction of phase 2, V2, as follows: 

h = (2 + 3V2)/5 
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These equations have been derived from the condition that Kerner’s 
equation® and eq. (1) should hold at the same time. Kerner’s equation re- 
lates the modulus of the two-phase system to the elastic constants of the 
continuous phase and the dispersed phase, assuming that the dispersed 
phase is in the shape of spherical particles and that there is perfect adhesion 
between the two phases. The equation is derived on the basis of the elastic 
theory. If the dispersed phase 2 takes a more extended shape, such as sheet 
or rod, the value of \ increases and that of ¢ decreases in comparison with 
the predicted values by eq. (2) or (3), that is, the tendency of the series- 
connection is intensified. 

Figure 3 shows 2’ and £” versus temperature at 138 cycles/sec. for the 
system in which the polybutadiene particles (phase 2) are dispersed in the 
styrene-acrylonitrile copolymer phase (phase 1), being graft-polymerized in 
their intersurface. The fu!! lines. in the figure are calculated by 
using the values listed in Table I for the parameters of the equivalent 
model of Figure 1. Satisfactory agreement between the observed and the 
‘alculated curves will be seen, except a slight deviation on the dispersion 
region. This tendency is duly ascribed to the unconttollable factors in 
polymerization. The simple equivalent model contains only one param- 
eter representing a mixing state, if the composition of the system is given. 
If necessary, eq. (1) can be extended to the generalized form to express 


a more complex behavior.*4 


TABLE I 


Volume ratio 





of 
(ST-AN)/PBD V.~PBD) d ¢ Curves 
77/23 (white circle) 0.23 0.600 0.380 Full line 
(0.54) 
56/44 (black circle) 0.44 0.775 0.568 Full line 


(0.66) 


* Values when the phase 2 is assumed to be dispersed in the spherical shape. 


3. Application of Equivalent Model Method to Crystalline Polymer 


It is assumed that the texture of the crystalline polymer is composed of 
the crystalline and the non-crystalline region. The non-crystalline region 
concerned with the primary dispersion is considered to exist separately from 
the lamella phase, being squeezed out from the lamella phase, as illustrated 
schematically in Figure 4. After the model of Figure 4, the lamella phase 
has a strong tendency to form a continuous crystalline phase, the support 
of which is also obtained from the mechanism of crystal growth. The left 
side of Figure 5 is a schematic representation of the crystalline texture, 
corresponding to the local structure given in Figure 4. Owing to these 
circumstances of the crystalline texture, the actual crystalline polymer will 
be reasonably represented by the simplified mechanical model, in which the 
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Center of 
* Spherulite 
Fig. 4. Model of the local structure of spherulite. Region A is the non-crystal- 
line region, region A’ is the defective region, and region C is the crystalline region with per- 
fect regularity 


Fig. 5. Schematic representations of the texture of the crystalline polymers with 
crystallinities of (a) about 90% and (b) about 50%, and their equivalent models. Black 
area expresses the non-crystalline region. 


non-crystalline region is dispersed in the crystalline phase as shown in the 
right hand side of Figure 5. In the same way as in the section 2, setting the 
non-crystalline region as phase 2 and the crystalline region as phase 1, the 
mechanical model can be constructed as shown in the left hand side of 


Figure 2. 
It should be noted that the viscoelastic behavior of the crystalline poly- 
mer cannot be interpreted unless the properties of the non-crystalline re- 
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Fig. 6. E’ and E” versus temperature for the model specimens constructed from a PVC 
film (curve 1) as the substitute for the crystalline phase and a PVC/NBR = 42/58 film 
(curve 7) as that of the amorphous phase. The models used for calculating the complex 
moduli are indicated in the figure. x is the volume fraction of PVC and corresponds to 
the crystallinity. PVC and NBR are the abbreviations of polyvinyl chloride and nitrile- 
butadiene rubber, respectively. 


gion can be clarified. Because of the difficulty of its clarification, we have 
inevitably adopted the model experiment. 

Figure 6 shows the viscoelastic behaviors of the model specimens, which 
are constructed by bonding together a polyvinylchloride (abbreviated as 
PVC) film as phase | to a film of the homogeneously mixed system of poly- 
vinylehloride and nitrile-butadiene ribber (NBR) as phase 2. The con- 
structions of the model specimens are shown in the figure. The PVC/NBR 
film is a molecularly mixed system of both polymers with its volume ratio of 
42/58, and gives a transparent film, showing a very broad shape of disper- 
sion. Curve 1 in Figure 6 corresponds to the dispersion curve of the PVC 
film, and curve 7 to that of the PVC/NBR film. The value of x cited in 
Figure 6 is the volume fraction of the phase 1 and corresponds to the degree 
of crystallinity, when this model specimen is considered to represent the 
crystalline polymer. In case of x taking 80-90% (curves 2 and 3), the 
crystalline dispersion becomes remarkable and this tendency accords with 
characteristics of linear polyethylene. In cases of x = 50-60% (curves 4 
and 5), the increase of intensity of the non-crystalline dispersion is seen 
in contrast with the decrease of the crystalline dispersion. These features 
are actually seen in the dispersion curves of polypropylene or branched 
polyethylene. Polyester, polyamide, and others with crystallinity of 
about 50% correspond to curve 6. Equivalent models employed here are 
selected not only to satisfy the composition of the system but also to satisfy 
the S-shape law,® which holds generally in the relation of log modulus vs. 
composition for the two-phase system. Figure 7 represents this relation. 
This relation comes from the fact that the dispersed phase and the medium 
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Oo O02 O04 O06 O8 10 
PVC/NBR(42/58) X pvc 
Fig. 7. Dependence of dynamic modulus at 40°C. upon the composition for the two- 
phase model specimen constructed by a PVC/NBR = 42/58 film and a PVC film in the 
form such as indicated in the right side of the figure. Model number corresponds to those 
in Figure 6. 
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Fig. 8. Model for calculating the strain-optical coefficient of the crystalline polymer. 


phase change their roles inversely in the intermediate compositions. The 
parallel and the series relation correspond to this limiting relation, respec- 
tively. From the results mentioned above, it may be said that the mechani- 
cal model in the left side of Figure 2 is the most simple mechanical model] of 
crystalline polymer in which the fine structure of crystalline polymer is 
taken into consideration. 


4. Application of Equivalent Mechanical Model to Dynamic Rheo-Optical 
Properties of Crystalline Polymer 


The viscoelastic properties of the crystalline polymers have a possibility 
to be interpreted by the equivalent model method, as mentioned before. 
The similar model method may be applied to the rheo-optical properties of 
polymers. In this case, the complex strain-optical coefficient of the crystal- 
line polymer will be expressed in terms of the complex strain-optical 
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coefficients of the crystalline and the non-crystalline region, the crystallin- 
ity and the texture parameter A or ¢g of the crystallized polymer. 

When a dynamic stress S* is applied to the model in Figure 8, the dy- 
namic strain y* is given as follows: 


i= - 
: U1 — dv) +B ‘J Ec | @) 


4 


where 2 and £4 are the complex moduli of the crystalline and the non- 
crystalline region, \ and ¢ are parameters mentioned above. Let the dy- 
namie strain of the upper element in the model be y; and the strain of the 
two lower elements ys, the stress S* satisfies the following equations. 
™ * 7y* ~ 
S* = "Eo (5) 
-y * , 7* > 
S* = y2{(1 — Ec + AEG} (6) 
From eqs. (4), (5), and (6), we obtain 
: (1 — A)Eo + AEA 
a ta... =. Eh ee ae aoe 
(1 —A+Ay)Ee + M1 — g)E4 


»s 
" Ec 


(8) 


We C= d+ OES + MI = QS 


Since the retardation of an element is proportional to the strain of the 
element and to the strain-optical coefficient corrected by the magnitude 
of the element, the retardation Aj of the upper element is given by 

Ai = (1 — ¢)Keri (9) 
The retardation A} of the lower elements is the sum of the quantities of 
both the elements connected in parallel, that is 


A> = (1 — \)eKers + AeK Ave (10) 


In eqs. (9) and (10), Ke and K are the strain-optical coefficients of the 
crystalline and the non-crystalline region. The first term of the right side 
of eq. (10) corresponds to the lower left element of Figure 8 and the second 
corresponds to the right element. The total retardation A* of the model is 
assumed to be the sum of the retardations of all elements. 


A* = Ai + Ag (11) 

From eqs. (9) and (10), the total retardation A* is 
A* = (1 — g)Kevi + {(1 — \eKe + AeKalye (12) 
Since the magnitude of the model is a unity, the retardation of the model 


is its birefringence. The strain-optical coefficient K* of the model is de- 
fined as the birefringence divided by the strain. 


K* = A*/y* (13) 
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From eqs. (7), (8), (12), and (13), we obtain the following expression for the 
strain-optical coefficient A*, 
lyhe 


5 s(K,—Ke 14 
(l1—A+Ag)Ec + 1 — v)E, (Aa dian 


K* = Ko + 


By rearranging this equation, we obtain 


_ Ke — K* _ heEc 
Ke —Ki, (1l-At+A—) Lot Al — OE: 


k* = (15) 


Bi icmpsteee 


- a, 4 * 16 
(1 = A)Ko + NE 4 ( » 


where V4 = Ag and #* is the non-dimensional complex number defined by 
eq. (15). 

Equation (15) represents the relationship between the viscoelastic proper- 
ties and the rheo-optical properties of the crystalline polymers. If the 
viscoelastic and the rheo-optical properties of the crystalline and the non- 
crystalline region of the real system are given respectively, the rheo-optical 
properties of the real system will be determined by these equations. 

It is noteworthy that the non-dimensional complex number R* is not 
equal to the volume fraction V4, but it depends on the moduli, the volume 
fraction of the non-crystalline region and the mixing state parameter. 
Strictly speaking, the frequently used eq. (17) is not correct. Instead of 
V, in eq. (17), R* must be used as expressed in eq. (18), which is derived 
from eq. (15). 


K* = Ko (1 — Va) + Vika (17) 
K* = Ke (1 — R*) + R*Ky (18) 


Characteristic properties of the crystalline region has a possibility to be 
evaluated by the direct measurement on a single crystal sample.’ If the 
properties of the non-crystalline region could be separated from those of the 
crystalline texture by use of an appropriately oriented material, the rela- 
tionship among the crystalline texture, the viscoelastic properties and the 
rheo-optical properties of the bulk crystallized polymer will be described by 
this method. 
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M. TAKAYANAGI, S. UEMURA, AND S. MINAMI 


Résumé 


On a construit un modéle mécanique pour un polymére cristallin en admettant qu’il 
se compose d’une phase cristalline et d’une phase non-cristalline. Nous avons introduit 
dans ce modéle des conditions compliquées au point que certaines parties de la phase 
cristalline seront fortement déformées par des tensions de concentration et d’autres ne 
le seront pas. En utilisant le modéle, on illustre la dépendance de la température vis-a-vis 
des modules complexes des polyméres cristallins pour différents degrés de cristallinité. 
Sur la base du modéle nous avons dérivé une équation permettant de calculer le coeffi- 
cient optique de tension du polymére cristallin a départ de ceux de la phase cristalline et 
non-cristalline. Cette équation fournit une relation entre les grandeurs rhéo-optiques, 
viscoélastiques et les paramétres représentant la structure fine ainsi que la cristallinité 
du polymére. De cette équation il est conclu que les coefficients optiques de tension ne 
peuvent étre exprimés par une relation additive simple comme cela est fait habituelle- 


ment. 


Zusammenfassung 


Unter der Annahme eines Aufbaus des kristallinen Polymeren aus der kristallinen und 
der nicht-kristallinen Phase wird ein mechanisches Modell fiir das kristalline Polymere 
konstruiert. In das Modell wurden so komplexe Verhiiltnisse eingebaut, dass z. B. 
gewisee Teile der kristallinen Phase bei einer Spannungskonzentration weitgehend de- 
formiert wurden und andere Teile nicht. Die Temperaturabhangigkeit des komplexen 
Moduls von kristallinen Polymeren mit verschiedenem Kristallinitatsgrad wird anhand 
des Modells erliutert. Auf der Grundlage des Modells wurde eine Gleichung zur 
Berechnung des spannungs-optischen Koeffizienten des kristallinen Polymeren aus den- 
jenigen der kristallinen und nicht-kristallinen Phase abgeleitet. Die Gleichung liefert 
eine Beziehung zwischen den rheo-optischen Gréssen, den visco-elastischen Gréssen 
und den Feinstruktur- und Kristallinititsparametern des Polymeren. Die Gleichung 
fiihrt zu dem Schulss, dass die spannungs-optischen Koeffizienten night durch die iiblich- 
erweise verwendete, einfache additive Beziehung dargestellt werden kénnen. 
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Synopsis 


The complex strain-optical coefficient, complex stress-optical coefficient, complex 
modulus of elasticity and the corresponding loss tangent of phase angles between the 
stress, strain and birefringence have been measured for several high polymers in the fre- 
quency range from 0.001 to 10 cycles/sec. Frequency and temperature dependences of 
the dynamic optical properties of vulcanized Hevea rubber are very similar to those of 
the dynamic mechanical properties. On the other hand, the dynamic optical properties 
of polyethylene and polypropylene show remarkable dispersions in the above frequency 
range. The dynamic optical behavior of these three materials are very typical and are 
quite different from each other. The differences have been explained by considering dif- 
ferent frequency dependences of three kinds of deformation mechanism, namely, orienta- 
tion of molecules in amorphous region, deformation of spherulites, and orientation of 


crystals. 


Introduction 


Introduction of dynamic birefringence into the high polymer study has 
recently been achieved by Stein, Onogi, and Keedy,'~* and the experiments 
performed with several high polymers based upon the theoretical founda- 
tion developed by them have revealed a sequence of novel information. 
One of the most important results obtained by them would be that the 
strain-optical coefficient, that is—the birefringence per unit strain, is de- 
pendent upon the frequency. This suggests that the dispersion as ob- 
served in dynamic modulus might be found with respect to the strain-optical 
coefficient. The effect of dispersion of modulus is verified by the measure- 
ment of the imaginary part of modulus, loss tangent or the phase difference 
between stress and strain by any means. However, the measurement of 
such quantity with respect to birefringence has not so far been achieved yet. 

The main purpose of this study is to measure the phase difference be- 
tween birefringence and strain, as well as the phase difference between 
birefringence and stress, and thereby to enrich our knowledge on the dis- 
persion effect, which has been confined to the conception of viscoelastic 
system. 

In analogy to the complex modulus in the viscoelastic system, which is 
given by 
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we define the following complex quantities with respect to the strain-optical 
and stress-optical coefficient as follows: 












A 
K* = — = K’'+ 1K” 
7 













A 
M* = = M’ — iM” 





A 





where S is the stress, 7 the strain, A the birefringence, and the star denotes 
the complex quantity. 

Like the loss tangent in the viscoelastic system, that is 
E”’ 
KE’ 








= tan 6 







we define the following quantities with respect to the strain-optical and 
stress-optical coefficient, respectively, 
x” 


| ed 









= tana 






and 







As 6 is the phase angle between stress and strain, so is a the phase angle 
between birefringence and strain, and 6 the phase angle between birefrin- 
gence and stress. As can be deduced from the expression of S/y, the stress 
is always ahead of the strain in the phase by 6, but the birefringence can be 
ahead of or behind the strain in phase, aecording as K”’ is positive or nega- 
tive. 

The experiment was performed with the specimens of crosslinked rubber, 
polyethylene, polypropylene, nylon 6, and polyethylene terephthalate. 
Kh’, E'', K'’, K", M’, M”, tan 6, tan @ and tan 6 of these specimens were 
measured at a certain frequency (1.50 cycles/sec.), and it was displayed that 
there are several different types of dynamic birefringence in reference to the 
relation among the vectors of stress, strain, and birefringence. Further, the 
frequency dependence of these quantities was measured within the fre- 
quency range from 10-2 to 10. Some of the materials manifested the dis- 
tinct dispersion or indication of dispersion of dynamic birefringence in this 
frequency range. The change in birefringence would arise from the line-up 
of the molecules in the amorphous domains, orientation of crystallites 
(spherulites) and deformation of crystallites (spherulites). The result of 
the experiment was discussed in this light. 
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Apparatus 


The phase difference was measured by the hysteresis loop method with the 
aid of a cathode ray oscilloscope, as well as by the direct method with the aid 
of an electromagnetic oscillograph. The construction of the apparatus and 
the block diagram are sketched in Figures 1 and 2. The periodical ex- 
tension of the sample was effected by the sinusoidal movement of the lower 


Fig. 1. Sketch of apparatus: (7) motor; (2) PIV; (3) worm gear; (4) cam; (4) strain 
gage, (6) spring; (7) sample; (8) sample chamber; (9) thermometer; (10) strain gage; 
(11) extension device; (12) worm gear; (13) motor; (14) lamp house; (14) lens and iris; 
(16) filter; (17) polarizer; (18) compensator; (19) analyzer; (20) viewer; (21) photo- 
multiplier and, (22) optical bench. 


clamp which is connected to an eccentric cam device rotating with a constant 
speed. The upper clamp is connected with a strain gage, by means of 
which the stress is changed into an electric signal. The movement of the 
lower clamp is also changed into an electric signal by means of another 
strain gage. The phase difference which is occasionally produced by the 
electric circuit itself was compensated by a compensating circuit. 

The optical system for the measurement of birefringence consists princi- 
pally of a light source (a high pressure mercury vapor lamp), filter, lens sys- 
tem, polarizer and analyzer which are crossed and placed before and behind 
the sample chamber, respectively, and a photomultiplier tube. The plane 
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Fig. 2. Block diagram: (A) strain gage; (B) sample; (C) cross nicols; (D) strain 
gage; (£) strain meter; (/) electric source; (G@) photomultiplier tube; (H) high 
voltage source; (J) ammeter; (J) strain meter; (K) electric source; (1) phase com- 
pensator; (7) cathode ray oscilloscope; and (NV) electromagnetic oscillograph. 


of polarization of the polarizer and analyzer is tilted 45° to the direction of 
stretching. A compensator is supplemented so that the measurement can 
be performed in the region where the relation between retardation and the 
current generated by the photomultiplier tube is practically linear. The 
monochromatic light of \ = 546 my was used. 

The sample chamber is furnished with two windows on the opposite 
sides to permit the light to pass through the sample. The chamber is 
sufficiently insulated, and the temperature in the chamber can be controlled 
between —40° and 200°C. 

IXvery time before starting the experiments, an elastic spring is employed 
in place of the sample, and the circuit is so adjusted that the linear Lissa- 
jous’ figure is obtained. Similar test was done also before starting the 
measurement of birefringence. In this case the area of the slit was changed 
sinusoidally in the same phase as that of the displacement of the lower 
clamp. 

Figures 3 and 4 show the Lissajous’ figure after phase-matching. Figure 
3 refers to the same phase vibration of force and deformation, and Figure 4 
to the same phase vibration of light intensity and deformation. 

The sample was put between the upper and lower clamps, and being 
slightly stretched, and left in the chamber for a certain period of time to 
permit the sample to relax sufficiently before starting the measurement. 

The amplitude of vibration was from 1 to 5% of the length of sample. 
The lengths of samples were from 6 to 8 cm. The frequency of vibration 
can be changed from 1/200 to 10 cycles/sec. 


Specimens 


(1) Rubber. Three parts of crosslinking agent (tertiary butyl-perbenzo- 
ate) and 0.5 part of antioxidant (2.6-ditertiary butyleresol) were mixed thor- 
oughly with 100 parts of pale crepe, and the blend was heated at 140°C. for 
30 min. The product was formed into a sheet having a thickness of 
0.45 mm. 
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Deformation 


Fig. 3. Lissajous’ figure of the same phase vibration of force and deformation 
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Fig. 4. Lissajous’ figure of the same phase vibration of light intensity and deformation. 


(2) Low-density polyethylene. Commercial pellets (Sumikathene 
I'204) were cast into a sheet having a thickness of 0.2 mm. by the inflation 
method. 

(3) Polypropylene. Commercial pellets of isotactic polypropylene 
(Avisan 1041) were extruded from a linear slit and formed into a sheet 
having a thickness of 0.35 mm. 

(4) Nylon 6. Commercial pellets were cast into a sheet in the same 
manner as before. The thickness of the sheet was 0.073 mm. 

(5) Polyester. Polyethylene terephthalate was cast into a sheet hav- 
ing a thickness of 0.073 mm. in the same manner as the above. 


Results 


A. Rubber 
The Lissajous’ figures of stress-strain, birefringence-strain and_bire- 
fringence-stress at a frequency of 1.50 cycles/sec. at 30°C. are given in 
Figures 5a, 5b, and 5c, respectively. As can be deduced from the direction 
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Fig. 5. Lissajous’ figures of (a) stress-strain; (b) birefringence-strain; ‘and (c) stress- 
birefringence of rubber. Arrows indicate the direction of scanning. 


LOGV 


Fig. 6. £’, B’’, and tan 6 versus log v curves of rubber. 
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LOG Vv 
Vig. 7. K’, E’’ and tan @ versus log » curves of rubber. 


Fig. 8. M’, M’’, and tan 6 versus log » curves of rubber. 


of scanning shown in l‘igure 5a, the phase of stress is ahead of the phase of 
strain. Similarly, Figure 5b leads to the conclusion that the strain is be- 
hind the birefringence. On the other hand, the phase difference between 
stress and birefringence is very small, as shown in Figure 5c, but careful 
observation reveals that the stress is ahead of the birefringence. 
Subsequently to the measurements at a constant frequency, the de- 
pendence of birefringence upon frequency will be treated below. Tligure 
6 shows the values of HL’, £’’, and tan 6 of rubber at 30°C. as a function of 
frequency. Within the region of frequencies covered by the experiment 
these quantities at 30°C. have only slight tendencies to increases with in- 
creasing frequency, but on the whole the curves are monotonous. Similar 
results are obtained for K’, K’’, and tan aas well as for 1/7’, 4/’’, and tan 8, 


as shown in Figures 7 and 8. 
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The analogy between modulus and strain-optical coefficient is seen also in 
the temperature dependence of these quantities. At temperatures below 
—20°C., the glass transition takes place, and the indication of dispersion 
comes in sight, that is—H’, 2”, and tan 6 measured at a frequency of 1.5 
cycles/sec. increases rapidly with decreasing temperature, as displayed in 


LOG (TAN 6) 
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Fig. 9. Log EZ’, log £”, and log tan 6 of rubber as a function of temperature. Frequency 
is 1.5 cycles/sec. 
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Fig. 10. Log K’, log K”, and log tan a of rubber as a function of temperature. Frequency 
is 1.5 cycles/sec. 


Figure 9. The same thing happens for strain-optical coefficient, K’, K”, and 
tan a, as shown in Figure 10. In a word, the dynamic birefringence of 
rubber as a function of frequency and temperature is quite similar to the 
dynamic modulus as a function of frequency and temperature. 
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Retardation 


Deformation Deformation 
(a) (b) 


Retardation 
(c) 
Fig. 11. Lissajous’ figure of (a) stress-strain; (b) birefringence-strain; and (c) stress- 
birefringence of low-density polyethylene. Arrows indicate the direction of scanning. 


B. Low-Density Polyethylene 


Considerably different results have been obtained with the specimen of 
low-density polyethylene. The Lissajous’ figures of stress-strain, bire- 
fringence-strain, and birefringence-stress at a frequency of 1.50 cycles/sec. 
at 30°C. are given in Figures lla, 11b, and llc. As can be expected, the 
stress is ahead of the strain (I*ig. 11a), but, contrary to the case of rubber, 
the retardation is behind the strain (lig. 116). In the case of rubber, the 
stress and birefringence were almost in the same phase, but in the case of 
polyethylene the stress is far ahead of the birefringence, as can be seen in 
Figure lle. In other words, 6 is the greatest and a is the smallest. 

With respect to the frequency dependence also remarkable difference from 
that of rubber is observed. The /’, FE”, and tan 6 versus frequency curves 
are shown in Figure 12. The tan 6 versus frequency curve seems to have a 
flat peak around the frequency of 0.1 cycles/sec. The relaxation spectrum 
measured by Faucher‘ with the low-density polyethylene also shows a flat 
peak at the time scale of 10sec. This may be regarded wishfully in the same 
light with the peak of tan 6 curve. The dynamic measurements of bire- 
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fringence offer the curves having more distinct character as can be seen in 
Figure 13, which shows K’, K’’ and tan a as a function of frequency. The 
most important point is that the K’’ versus frequency and tan a@ versus fre- 
quency curve have a distinct peak at about 0.1 cycles/sec. Stein, Onogi, and 
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Vig. 12. 2’, £’’, and tan 6 versus log v curves of low-density polyethylene. 
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Fig. 13. A’, A’, and tan @ versus log v curves of low-density polyethylene. 


Keedy** suggested preliminarily by their measurement of K’ that the dis- 
persion might happen at this frequency, but the measurement of K’’ and 
tan a as a function of frequency as shown in Figure 13 seems to authenticate 
their opinion. The tan 8 versus frequency curve shown in [igure 14 also 
suggests the dispersion as happening at the frequency of about 0.1 cycles/ 


sec. 
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Fig. 14. M’, M’’, and tan 6 versus log » curves of low-density polyethylene. 


C. Polypropylene 


The strain-optical coefficient of polypropylene is characteristic in the 
point that it is negative-signed within the region of small strain, as is the 
case in the present experiment. The observation of Lissajous’ figures 


shows that 6 is the greatest and 8 is the smallest. In this meaning, the 
relation is the same as that of rubber, but 6 is not so small as in the case of 


rubber. 

In Figure 15 the relation between 6, a, and @ is displayed schematically 
together with the vectors of stress, strain, and birefringence, and the values 
of tan 6, tan a, and tan 8 are given for the reference. 


RUBBER 


S 82a>B 
tans 1.68-10°7 
tana |,62-" 
tanB 6.50-10°4 


POLYETHYLENE 


Ss B>8>a a 
tan8 19.4-10 
tana 6.50-" 
tang 26.2: " 


POLYPROPYLENE 


8>a>B 

ton’ 12.8-10% 
tana 9.0-" 
tanB 3.76-" 


Fig. 15. Vectors of stress, strain, and birefringence, and the relation between angles of 
phase difference, 5, a, and 6. The values of loss tangents are given 
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The frequency dependence of 2’, ’’ and tan 6 is shown in Figure 16. 
The tan 6 versus frequency curve displays a tendency to increase both 
towards the higher frequencies and towards the lower frequencies. This 
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Fig. 16. 2’, E’’, and tan 6 versus log » curves of polypropylene. 
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Fig. 17. K’, K’’, and tan 8 versus log » curves of polypropylene. 


seems to imply that the dispersions would take place in the regions of both 
the higher and lower frequencies than those covered by this figure. Figure 
17 presents the K’, K’’, and tan a versus frequency curves. K’ and K”’ 
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Fig. 18. Log K’ versus log » curves of nylon 6, polyethylene terephthalate, rubber, 
polyethylene, and polypropylene. 


have negative values. The shape of the curve of tan a versus frequency is 
very similar to that of the tan 6 versus frequency curve. 


D. Nylon 6 and Polyester 


Figure 18 presents the plots of log K’ versus log v of nylon 6 and poly- 
ethylene terephthalate, together with those of rubber, polyethylene, and 
polypropylene. In each case the dependence of K’ upon frequency is not 


noticeable in the region of frequencies covered by this experiment. This is 
the case also with the frequency dependence of FL’. 


Discussion 


With reference to the frequency dependence of birefringence and loss 
tangent three types of response of high polymers to vibration are dis- 
tinguished. 

(1) As can be seen with rubber, K’ is positive and dependent only little 
upon the frequency. The stress and birefringence are almost the same in 
phase, and tan 8 is extremely small. 

(2) A distinct dispersion is observed within the limit of frequencies cov- 
ered by this experiment. Both the tan a versus frequency curve and 
tan 8 versus frequency curve show a peak almost at the same frequency. 
The birefringence is far behind the strain. This is the case with low-density 
polyethylene. 

(3) As represented by polypropy!ene K’ is negative, and the absolute 
value increases with the increase in frequency. The phase of strain is be- 
hind those of stress and birefringence. 

The first type of response is characteristic of the rubber state. The speci- 
men of rubber employed in this study behaves as a material having the ideal 
rubber elasticity at the ordinary temperatures, and birefringence arises 
mainly from the orientation of molecules to the direction of stretching, and 
should be almost independent of the frequency at lower frequencies covered 
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by this experiment. Since the deformation is mostly elastic, not only the 
phase difference between stress and strain, but also the phase difference be- 







tween stress and birefringence is very small. 

The response of the second type involves the deformation characterized 
by the relaxation time comparable with the time scale of the present experi- 
ment; that is 10~! to 10? see. The distinct dispersion which was revealed 
very clearly by the dynamic birefringence measurement is related to this 
mechanism. The relaxation time of orieritation of crystallites is naturally 
shorter by far than that of the orientation of molecules, and is postulated 
to be the order of 10-100 sec. Therefore, the dispersion as found in this 
study with polyethylene may possibly be ascribed to the rotation of 










spherulites. 

The value of K’ decreases with increasing frequency, and this happens the 
more distinctly, the greater is the crystallinity and the greater the size of 
crystallites. The rotation of the spherulites can take place only slowly and, 
therefore, when the frequencyNincreases, the motion of spherulites can not 
follow the rapid vibration, and consequently K’ decreases. In a word, the 
orientation of molecules and orientation of spherulites are the response of 
polyethylene to the external force. 

The interpretation of the third type of response as is found in polypro- 
pylene is very difficult at present. According to Sasaguri, Hoshino, and 
Stein,’ the birefringence which would arise from the deformation of spheru- 
lites is negative-signed. Since the relaxation time of the deformation of 
spherulites is assumed to be longer by far than the time scale of this experi- 
ment, the birefringence should be influenced only little by the change in 
frequency. The K’ versus frequency curve of polypropylene would be the 
arithmetic sum of the negative contribution due to the crystal deformation 
and the positive contribution arising from the orientation of molecules and 




















that of crystallites. 
We postulate that K’ is composed of three types of contribution, that is 


K’ = K', + K'p + Ky 


where K’, is the contribution from the orientation of molecules in amor- 
phous domains, K’y the contribution from the deformation of spherulites 
and Ko’ the contribution from the orientation of crystallites. 

In the first type of response (rubber) K’ is almost equal to K’,4. In the 
second type of response (polyethylene), three terms of the right-hand side of 
the equation have significance, but in the region of frequencies covered by 
the experiment K’, and K’o are predominant, and in the region of lower 
frequencies K’y would surpass the K’,. In the third type of response 
(polypropylene), all the three terms are significant. However, the fact that 
K’ is negative suggests the K p to be most. predominant. 
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Résumé 


Le coefficient optique complexé de déformation, le coefficient optique complexe 
de tension, le module complexe d’élasticité, la tangente correspondante a l’angle de 
déphasage entre la tension et V’élongation, et la biréfringence ont été étudiés pour 
plusieurs hauts polyméres a des fréquences variant de 0.001 A 10 cycles/seconde. 
La dépendance des propriétés optiques dynamiques de la fréquence et de la température 
sont trés semblables 4 celles des propriétés mécaniques dynamiques. D’autre part, les 
propriétés optiques optiques dynamiques montrent des dispersions remarquables dans 
le domaine de fréquence rapporté. Le comportement optique dynamique de ces trois 
matériaux est trés caractérisé et tout a fait différent l'un de l’autre. On explique ces 
différences en considérant que la dépendance de la fréquence varie pour chacun des 
trois mécanismes de déformation envisagés a savoir l’orientation des molécules dans la 
région amorphe, la déformation des sphérulites et l’orientation des cristaux. 


Zusammenfassung 


An einigen Hochpolymeren wurden im Frequenzbereich von 0,001 bis 10 Hertz der 
komplexe dehnungs-optische Koeffizient, der komplexe spannungs-optische Koeffizient, 
der komplexe Elastizitiitsmodul sowie der entsprechende Verlusttangens des Phasen- 
winkels zwischen Dehnung, Spannung und Doppelbrechung gemessen. Frequenz- und 
Temperaturabhiingigkeit der dynamisch-optischen Kigenschaften von vulkanisiertem 
Heveakautschuk sind derjenigen der dynamisch-mechanischen Eigenschaften sehr 
iihnlich. Andrerseits zeigen die dynamisch-optischen Eigenschaften von Polyiithylen 
und Polypropylen im obigen Frequenzbereich eine bemerkenswerte Dispersion. Das 
dynamisch-optische Verhalten dieser drei Stoffe ist recht typisch und in jedem Fall ganz 
verschiedenartig. Die Unterschiede wurden durch Annahme einer verschiedenen Fre- 
quenzabhiingigkeit der drei Arten von Deformationsmechanismen, niimlich Molekiilo- 
rientierung im amorphen Bereich, Sphiirolithdeformation und Kristallorientierung, 


erkliirt. 
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Dynamic Birefringence of Polyolefins* 


K. SASAGURI} and R. 8. STEIN, Polymer Research Institute, 
University of Massachusetts, Amherst, Massachusetts 


Synopsis 


The use of the dynamic birefringence method was extended to the study of the rela- 
tionship between structural changes and mechanical properties for polyethylene, atactic 
and isotactic polypropylene, and polybutene-1. The phase angle between birefringence 
and strain was measured which permitted the resolution of the strain-optical coefficient 
into real and imaginary parts. The real and imaginary parts of the dynamic relaxation 
modulus were also measured and compared with the optical data. Results for the 
crystalline polymers were interpreted in terms of the two-stage mechanism of deforma- 
tion, involving spherulite deformation followed by internal crystalline reorientation. 
For the atactic polymer, the viscoelastic flow of amorphous regions contributes to the 
birefringence. The temperature dependence of birefringence and stress is similar for 
this mechanism but not for the crystalline contribution. 


Introduction 


In previous papers, a linear phenomenological theory for time-dependent 
birefringence was presented and dynamic birefringence data for several 
polymers were reported.'~> The birefringence values for polyethylene, 
polypropylene, and polybutene-1 were found to be time-dependent, both in 
experiments in relaxation at constant length, at constant strain rate, and 
dynamic experiments where the sample is subjected to vibrational strain at 
varying frequency. All of these experiments were shown to be inter- 
related through a strain-optical coefficient relaxation spectrum. This spec- 
trum characterizes the distribution of optical response mechanisms of the 
crystalline polymers which is dependent upon the structure and morphology 
of the polymer as affected by thermal and mechanical treatment and by the 
temperature of the experiment. The response has been shown to vary with 
aging of polypropylene and polybutene-1 samples as influenced by sec- 
ondary crystallization and crystal structure transitions. This time-de- 
pendence of birefringence was related to the characteristic time required 
for crystal orientation following deformation of the polymer. 

A theory of crystal orientation involving spherulite deformation followed 
by rearrangement of crystals and recrystallization within the deformed 

* Supported by grants from the Air Force Office of Scientific Research, the Army 
Research Office (Durham) and the Petroleum Research Fund, and by a contract with 
the Office of Naval Research. 

T On leave from the Ashai Chemical Company, Japan. 
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spherulites has been proposed.**7 Direct microscopic observation of the 
spherulite deformation has been made and related to changes in the low- 
angle light scattering pattern.? A measurement of the rate of this spheru- 
lite deformation process by the laser light scattering technique has been 
reported.’ An attempt has been made to relate the predictions of this 
theory to changes in crystal orientation functions measured by x-ray dif- 
fraction.’ It is useful to interpret the dynamic birefringence in terms of this 









mechanism. 

In this paper, further experiments on the dynamic birefringence of poly- 
olefins in periodic strain, relaxation at constant length, and constant rate 
of strain experiments are reported. Dynamic mechanical properties are 
studied simultaneously in periodic strain experiments and their relationship 
to optical properties is discussed. Phase relationships between stress and 
strain and birefringence have been measured and are in agreement with 
predictions of the theory and experimental measurements reported by 
Yamada, Hayashi, Onogi, and Horio.'° 
















Polyethylene Dynamic Birefringence 





The time dependent birefringence of polyethylene films was studied at 
room temperature in vibrational experiments. Strain-optical coefficient 
for low density quenched, low density annealed, and high density quenched 
films of polyethylene are shown as a function of frequency in Figure 1. 
Structural properties of these polymers are listed in Table I, where density 
was measured by density gradient column, crystallinity by x-ray diffrac- 
tion, applying Hermanns’ method'*:'* and crystallite size from half-width 












of (110) diffraction. 


TABLE I 














Crystal- Cryst. size, 


Density linity, % A. 




















Low density quenched polyethylene 0.911 44.4 78 
Low density annealed polyethylene 0.919 47.4 89 
High density quenched polyethylene 0.946 66.9 118 








The strain-optical coefficient decreases with increasing frequency. This 
is related to the time dependence of crystal reorientation within the de- 
formed polyethylene spherulites. The spherulite deformation theory indi- 
cates that the birefringence change produced by affinely deforming a poly- 
ethylene spherulite without any internal crystalline rearrangement would 
be negative. The fact that a positive birefringence is always observed 
would indicate that there is sufficient internal crystal mobility within these 
polyethylene spherulites to allow sufficient change in the crystal orientation 
with deformation from a b-axis radial to a c-axis radial state to always pro- 
duce positive birefringence. That is, if the parameter y of the spherulite 
deformation theory which characterizes internal mobility of crystals*® is 
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sufficiently great, the predicted birefringence-elongation curve will always 
be positive. This high mobility of polyethylene crystals is associated with 
the probable ease with which the crystals may deform by means of disloca- 
tion migration mechanisms. There has been experimental indication for 
this in the electron microscopic studies of Geil, demonstrating deformation 
of polyethylene single crystals.' 

With increase in degree of crystallinity, crystal size, and density, the 
crystal dispersion is shifted to lower frequencies indicating a greater re- 
laxation time for the crystal motion. There is also a greater decrease in 
birefringence with increasing frequency under these conditions. 

By measurement of the phase angle between strain and birefringence, it 
was possible to resolve the dynamic strain-optical coefficient into its real 
and imaginary parts, K’ and K”. This was accomplished for the larger 
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Fig. 1. The variation of the dynamic strain-optical coellicient for 3 types of poly- 
ethylene with frequency. 


phase differences (with tan a > 0.2), by directly comparing the tracings of 
strain and birefringence obtained with the two channel Sanborn recorder 
of the apparatus which has been described.? lor smaller phase differences, 
a Lissajous’ figure technique was used similar to that employed by Ya- 
mada, et al."° The strain signal, amplified by the D.C. amplifier of a Ballan- 
tine Model 643 VT VM, and the birefringence signal, amplified by a General 
Radio Model D.C. Electrometer Amplifier, were placed on the X and Y de- 
flecting plates of a Tektronix Model 502 Oscilloscope. With no phase dif- 
ference, a 45° straight line trace is obtained, while with a finite phase dif- 
ference, an elliptical figure results. The phase balancing of the two circuits 
was checked by vibrating a quartz wedge in the optical path where a 
birefringence signal certain to be in phase with the strain is obtained. This 
gave the expected 45° straight line. 
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The phase angle a between birefringence and strain was determined in 
the conventional fashion from the Lissajous’ figure by the equation 









sin q = iA” Ay) (1) 










where Ao is the maximum amplitude of the oscilloscope trace in a particular 
direction (say the x direction) and A* is the intercept of the Lissajous’ figure 
with the coordinate axis in this direction. 

The real and imaginary parts of the complex strain-optical coefficient K* 
(ratio of maximum birefringence to maximum strain) may be obtained by 











K’ = K* cosa (2) 






and 









(3) 


= K* sin a 











0. 
FREQUENCY 


0.01 





CPS 


The variation of the dynamic strain-optical coefficient of a medium density 
quenched polyethylene with frequency at 20°C. 









A plot of the variation of the real and imaginary parts of the strain- 
optical coefficient with frequency is given in Figure 2 for a quenched 
medium-density polyethylene film at 20°C. It is seen that the imaginary 
component exhibits a maximum in the region where the real component is 
decreasing with frequency as predicted by the phenomenological theory. 
For comparison, the real and imaginary components of the complex 
dynamic modulus were determined on the same sample with the same ap- 
paratus over the same frequency range. A stress measuring device con- 
sisting of a Schaevitz TDC-4M miniature dynamometer is placed in series 
with the upper clamp. This was excited by a 10 K.C. carrier signal which 
was amplified by a Ballantine Model 643 amplifier demodulated with a 
solid state diode circuit and compared with the strain signal using the 
Lissajous’ figure technique on the oscilloscope. Phase balance was checked 
by determining whether a 0° phase difference was obtained when a spring 
The real and imaginary modulus was 
















was connected between the clamps. 
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calculated from the ratio of the magnitude of the stress and strain and the 
phase difference was obtained by conventional techniques. 

The real and imaginary moduli are plotted in Figure 3 for the same 
sample as in Figure 2. A maximum is seen in 2” at about the same fre- 
quency as for K” but the 2” maximum is lower. It is suggested that the 
K” maximum is associated with the time for crystal orientation change 
within the spherulite involving both the motion of crystals within the 
amorphous matrix as well as changes within the crystals themselves involv- 
ing plastic flow, defect migration, and recrystallization. These changes 
apparently affect the birefringence more then they do the stress, resulting 
in the higher K” maximum. 
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Fig. 3. Dynamic modulus and loss modulus-frequency plots for low density poly- 
ethylene quenched film at 20°C. 


It is also seen that while K” falls off at higher frequencies, ” does not 
and appears to be increasing toward a second high frequency maximum. 
This second maximum is likely to be associated with motion in amorphous 
or disordered regions. The strain-optical coefficient for amorphous regions 
is known to be less than that associated with crystal orientation. 

These experiments may be interpreted in terms of the phenomenological 
theory of dynamic birefringence?’ in terms of which 
: Ag(In 7)w?r? + Bo (In 7) 

—© 1 + w?r? 


& f [Ao(In T) = Bo (In T) lwr 
sg 1 + w*r? 


K* = K’' + iK’ 
K(w) = |K*| = [(K’)? + (K")?]” 


din +r (4) 


dint 


tan a = K°/K’ 
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The determination of the variation of Bo(In 7) — A (In 7) with In 7 has 
been discussed. This may be approximated by 


Ao(In 7) — Bo(In 7) = (21K eons al) im 
= w=1/r 


By this technique, it was shown that for polyethylene, Bo(In +) > Ao(In 7) 
over the entire range of measurements. Apo (In 7) is the strain-optical co- 
efficient associated with the elastic response while Bo(In 7) is that associated 
with viscous processes. In the case of polyethylene, it was concluded that 
the positive viscous birefringence contributed by the c-axis orientation of 
crystals upon stretching is greater than the elastic contribution arising from 
amorphous chain orientation and spherulite deformation. It is conceivable 
that at sufficiently high frequencies or low temperatures, the negative 
elastic contribution for spherulite deformation may dominate, but there 
is no indication of this over the present range of measurements. Since By 
> Ao, both K’ and K” are positive, giving rise to a positive tan a and over 
the accessible frequencies. If at sufficiently high frequencies or low tem- 
peratures, Ay is negative and Ayw*r? > Bo, then A’ will reverse sign and tan 


a will become negative. 
Dynamic Birefringence of Polybutene 


The optical behavior of polybutene-1 was reported previously.*!4 The 
birefringence of polybutene-1 film becomes negative on stretching due to 
b-axis orientation before necking. The deformation of the spherulites con- 
taining stable modification-I crystals occurs elastically with little reor- 
ganization leading to negative total birefringence. Dynamic birefringence 
of polybutene-1 has been shown to be independent of frequency over the 
range of 0.01-10 cycles/sec. and is negative.‘ 
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Fig. 4. The variation of dynamic modulus and loss modulus of polybutene-1 with 
frequency at 20°C. 
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FREQUENCY OR 1I/RELAXATION TIME 


Vig. 5. A postulated variation of (Ao — Bo) with linear 7 and of K’, K”, and tan a@ 
with log (1/w) for aged isotactic polybutene-1 and polypropylene. 

In Figure 4 dynamic modulus /’ and loss modulus /” for polybutene-1 
measured at 20°C. are plotted against frequency. Increases in both the 
dynamic modulus and the loss modulus with frequency are observed. 
These changes may be related to the so-called 6-transition which is located 
in the neighborhood of — 10°C." and associated with motion of amorphous 
regions. The mechanical loss occurring in these regions of negligible 
change in dynamic birefringence indicates again that the amorphous 
transitions affect modulus much more than birefringence. 

For this polymer there are no mechanisms having an appreciable strain- 
optical coefficient for which there is a relaxation time corresponding to 
(1/w) over the range of frequencies of the measurement. The measured 
negative strain-optical coefficient is for the elastic mechanism of spherulite 
deformation which has a relaxation time rz < 0.01 sec. for aged polybutene- 
1. The relaxation time for internal crystal rearrangement, 7, > 10* sec.: 
Consequently, the spectrum of (By) — Ay) might appear as in Vigure 5 to- 
gether with the accompanying values of K’, K", and tan a. This figure will 
also be referred to in discussing the optical properties of the other polymers. 


Dynamic Birefringence of Isotactic Polypropylene 


The variation of birefringence of isotactic polypropylene during exten- 
sion, relaxation at constant length, retraction, and during vibrational ex- 
periments has been reported.‘ These results are consistent with a two- 
stage deformation process involving first, the elastic deformation of the 
spherulites contributing a negative birefringence, followed by a viscoelastic 
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Fig. 6. The variation of strain-optical coefficient or isotactic polypropylene with 
frequency at various temperatures. 
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Fig. 7. Frequency dependency of loss modulus of isotactic polypropylene at various 


temperatures. 


rearrangement of the crystals within the spherulites making a positive 
birefringence contribution. Relaxation at constant length involves an in- 
crease in birefringence with time which has been associated with an increas- 
ing positive contribution resulting from the viscoelastic intra-spherulitic 
rearrangement. The birefringence initially increases on retraction as a re- 
sult of the decreasing negative elastic contribution, and then decreases as 
a result of the decreasing positive contribution from internal crystal orienta- 
tion. The strain-optical coefficient in vibrational experiments is negative 
since the negative elastic deformation contribution predominates at the 
small strains involved. It becomes less negative with decreasing frequency, 
however, because of the increasing positive contribution from the intra- 
spherulitic rearrangement which is only seen at longer times. This fre- 
quency dependence was found to be greater for aged samples and at lower 
temperatures where the internal mobility of the crystals is less. 

This temperature dependence is seen in the strain-optical coefficient data 
in Figure 6 where the strain-optical coefficient of an aged isotactic poly- 
propylene sample is plotted against frequency at a number of temperatures. 
The coefficient is seen to be negative at all temperatures and is more tem- 
perature dependent at the lower temperatures. The curves are not super- 
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posable. The temperature dependence is not uniform but is greater in the 
interval of 20°-40°C. indicating the presence of an optical loss maximum 
in this region. This is also seen in the room temperature measurements of 
Yamada, et al.!° which exhibit a negative K’ which decreases with increas- 
ing frequency and a K” maximum at about 1 cycle/sec. The strain-optical 
coefficient appears to be again changing greatly with temperature at tem- 
peratures in the region of 100—120°C. or greater. 

The mechanical dynamic loss modulus /” for this polymer is plotted in 
‘igure 7 for this range of temperature. The loss maxima shift toward higher 
frequencies with increasing temperature. Loss maxima are seen in the 
region of 1 cycle/sec. at 10°C. and at about 90°C. The lower temperature 
maximum has been associated with the onset of micro-Brownian motion of 
amorphous chains, '*~'® and has been shown to decrease in intensity with de- 
creasing atactic content.’ It is likely that the contribution of crystal orien- 
tation to the strain-optical coefficient is much greater than that of amor- 
phous orientation in this region, so that the effect of this amorphous transi- 
tion on the optical properties is believed to take place through the effect of 
these amorphous regions on the ability of the crystals to orient. That is, the 
effective viscosity of the amorphous regions for crystal orientational motion 
is changing in this region. 

The optical behavior of isotactic polypropylene may be understood in 
terms of Figure 5. The relaxation time for crystal rearrangement is shorter 
for isotactic polypropylene than for polybutene-1 so that the frequency at 
which the optical dispersion for this process occurs is within the range of 
normal measurements. Consequently, the increase in negative K’ with in- 
creasing frequency arises from this higher relaxation time mechanism and 
leads to a negative tan a. It is conceivable that at larger relaxation times 
and higher temperatures, K’ will reverse in sign and become positive and 
tan a will go to — ~ and then reverse to + © as a goes through 90° as will 
be seen for atactic polypropylene. Thus far this has not been observed for 
isotactic polypropylene. 

The higher-temperature mechanical loss maximum is associated with the 
onset of motion within crystals” and is found to increase in intensity with 
decreasing atactic content. It seems possible that such motion enables 
crystal orientation to occur through plastic deformation involving motion 
of chains within the crystal with respect to each other. 


Dynamic Birefringence of Atactic Polypropylene 


Atactic polypropylene samples obtained by elution fractionation were 
used, which were considered to contain 3—5% soluble low molecular weight 
crystalline species. If this atactic film is stretched and held at constant 
length, birefringence decays rapidly to a positive value at a low tempera- 
ture and to a negative value at a higher temperature as shown in igure 8. 
At a much higher temperature, birefringence is initially negative and de- 
cays to a much more negative value. 
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Fig. 8. Birefringence relaxation and recovery curves for atactic polypropylene. 
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Fig. 9. Temperature dependency of equilibrium birefringence of atactic polypropylene. 


If the load is removed after the birefringence has leveled off, the bire- 
fringence decreases instantaneously to a negative value and then increases 
rapidly towards zero. Similar results were observed by Tsvetkov and his 
coworkers for atactic polypropylene fraction.'* He attributed this negative 
birefringence to a small amount of isotactic species in the atactic poly- 
propylene, which cannot be revealed by x-ray. This might be reasonable 
also in our case because small amounts of isotactic polypropylene in atactic 
polypropylene can be detected by infrared spectra.'® 

Thus the initial positive birefringence is caused by amorphous chain 
orientation, and birefringence decays with the relaxation of amorphous 
chains with time. In addition to this amorphous relaxation, orientation of 
microcrystallites formed by small amounts of isotactic polypropylene 
species probably leads to negative equilibrium birefringence. 

This equilibrium birefringence was studied further by heating and cooling 
the sample at a constant length and measuring birefringence changes. 
This equilibrium birefringence is decreased to a negative value by heating 
the sample at constant length to 60°C. On cooling, birefringence of the 
sample is observed to be constant (lig. 9). This result indicates that this 
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Fig. 10. Master curves of dynamic birefringence K’, K”, and tan a@ of atactic 
polypropylene plotted against frequency. 
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Fig.11. Master curves of dynamic modulus and loss modulus plotted against frequency 


negative birefringence arises from an irreversible process and may be due 
to microcrystallite orientation in the amorphous medium. 

The dynamic birefringence of atactic polypropylene was studied at vari- 
ous temperatures with simultaneous measurements of dynamic stress. 
The curves obtained were superposed by shifting along the In 7 axis by 
adding factors In a, at each temperature to obtain the master curves which 
are shown in Figures 10 and 11. 

In Figure 10 the dynamic-optical coefficients K’ and K” and the tangent 
of the phase angle a are plotted against reduced frequency, where the 
reference temperature is 30°. The imaginary part of the dynamic-optical 
coefficient K” is always negative and decreases with increasing frequency. 
The real part of the dynamic-optical coefficient K’ is initially negative at 
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Fig. 12. The calculated reduced strain-optical spectrum for mostly atactic 
polypropylene. Reference temperature 30°C. 
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Fig. 13. Shift factors versus (1/7') for birefringence and stress of atactic polypropylene. 


low frequency or at high temperature and increases rapidly with frequency. 
The birefringence is observed to change sign at a frequency which is re- 
ferred to as the critical frequency. 

The relaxation curve at 0°C. shows one optical dispersion which may be 
the glass transition of the amorphous chains. The tangent of the phase 
angle a is initially negative at low frequency or at high temperature and in- 
creases to infinity at the critical frequency. Above this critical frequency, 
tan a changes sign and increases from negative infinity. This behavior may 
be understood in terms of the optical relaxation time distribution for this 
polymer shown in Figure 12 which has been calculated from the experimen- 
tal data. The broad optical relaxation time distribution is that associated 
with the viscoelastic flow of the amorphous chains and is similar to the 
mechanical relaxation spectrum for linear amorphous polymers.” It has 
been shown?! that the birefringence and stress should be proportioned in the 


DYNAMIC BIREFRINGENCE OF POLYOLEFINS 151 


rubbery region for such viscoelastic flow. The strain-optical coefficient is 
positive for this orientation of amorphous polypropylene chains. 

The reversal in sign of birefringence at long times or low frequencies is a 
consequence of a long relaxation time mechanism with a negative strain- 
optical coefficient which, as indicated before, is associated with the negative 
orientation of microcrystalline structure contributing negatively to By but 
not appreciably to Ao. When the real part of the strain-optical coefficient 
changes sign and goes through zero, it is apparent from eq. (5) that tan a 
will pass through infinity. This “maximum” in tan a does not signify any 
molecular loss mechanism at this frequency but is merely an accident of the 
difference in sign of birefringence of competing mechanisms. 

The temperature dependence of both the birefringence and the stress of 
atactic polypropylene may be expressed in terms of the variation of the 
shift factors a; required to produce superposition in Figures 10 and 11 with 
temperature. A plot of log a7 against (1/7') is given in l’igure 13 for stress 
and for birefringence. The variation is similar for both indicating that the 
same molecular mechanism is responsible for the temperature dependence 
of both phenomena. This is the approach to a glass transition at the higher 
frequencies and lower temperatures associated with the loss of motion of the 
main chains. 
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Résumé 


L’usage de la méthode de biréfringence dynamique a été étendu 4 |’étude de la relation 
entre les changements de structure et les propriétés mécaniques pour le polyéthyléne, 
le polypropylene atactique et isotactique et le polybuténe-1. L’angle de phase entre la 
biréfringence et |’élongation est mesuré, ce qui permet la résolution du coefficient optique 
de |’élongation en parties réelles et imaginaires. Les parties réelles et imaginaires du 
module de relaxation dynamique ont été mesurées et comparées avec les données op- 
tiques. Les résultats pour les polyméres cristallins sont interprétés en termes d’un 
mécanisme de déformation en deux étapes, comprenant une déformation du sphérulite 
suivi d’une réorientation cristalline interne. Pour le polymére atactique |’écoulement 
viscoélastique des régions amorphes contribue a la birtiringence. La dépendance vis-a- 
vis de la température de la biréfringence et de la tension est semblable pour ce méca- 
nisme, mais pas en ce qui concerne la contribution cristalline. 


Zusammenfassung 


Die Verwendung der dynamischen Doppelbrechungsmethode wurde auf die Unter- 
suchung der Beziehung zwischen Strukturiinderungen und mechanischen Eigenschaften 
fiir Polyathylen, ataktisches und isotaktisches Polypropylen und Polybuten-1 ausgede- 
hut. Der Phasenwinkel zwischen Doppelbrechung und Verformung, der die Trennung 
des spannuags-optischen Koeffizienten in Real- und Imaginiirteil gestattet wurde gemes- 
sen. Auch der Real- und Imaginirteil des dynamischen Relaxationsmoduls wurde 
bestimmt und mit den optischen Daten verglichen. Die Ergebnisse fiir die kristallinen 
Polymeren wurden auf Grundlage des Zweistufenmechanismus der Deformation mit 
Spharolithdeformation und nachfolgender innerer kristalliner Reorientierung gedeutet. 
Beim ataktischen Polymeren triigt das viscoelastische Fliessen der amorphen Bereiche 
zur Doppelbrechung bei. Die Temperaturabhingigkeit der Doppelbrechung und Span- 
nung ist bei diesen Mechanismen jihnlich, nicht aber beim kristallinen Anteil. 
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Rheo-Optical Properties of Polymers. X. Relaxation 


D. G. LEGRAND and W. R. HAAF, General Electric Research Laboratory, 
Schenectady, New York 


Synopsis 


Light scattering relaxation in thin strained polyethylene films has been determined. 
Results for swollen and unswollen films are presented. Concurrent studies of stress re- 
laxation indicate that the mechanism by which relaxation takes place is associated with 
changes in the average fluctuation in refractive index of the material. The ratio of the 
change in the fluctuation in refractive index to stress appears to be independent of 


time. 


INTRODUCTION 


In previous papers, we have presented data on time dependent studies 
of light scattering from thin polymer films of polyethylene.':? Due to the 
necessity of working at large deformations, the results were difficult to 
interpret. Due to many modifications in the original experimental appa- 
ratus, as well as the use of samples prepared under very carefully controlled 
conditions, which will be described in detail elsewhere,* we have been able 
to obtain more reproducible data which can be discussed in terms of current 
theories and which elucidates some of the mechanisms by which relaxation 
takes place in polyethylene. 


Experimental 


In order to perform relaxation experiments, we use the straining system 
shown in Figure 1. The sample § is clamped on two blocks D and F, 
which are caused to move apart via a screw A, which has a left-handed 
thread at one end and a right-handed one at the other. A Heller thyratron 
DC motor is connected to the screw A. By connecting the motor control 
to a timed circuit, the motor can be turned on and off for any period of time 
between 1 msec. to 10 min. This unit can be placed in the light scattering 
apparatus either with the stretching axis vertical or horizontal to the 
scattering plane. A displacement transducer measures the rate of separa- 
tion of the clamps. Since the sample is deformed from both ends, the 
center of the sample stays in the same position at all times. The rate of 
separation of the clamps was measured by the displacement transducer T. 

Prior to usage, each sample was examined for microscopic defects and 
orientation in a polariscope. In order to eliminate surface scattering, a 
drop of silicone oi] was placed on each surface; thin glass microscope cover 
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Fig. 1. Straining system. 


slides were then placed on top of the silicone oil. With a judicious applica- 
tion of silicone oil, it was found that the presence of the cover slides did not 
affect the results since reruns on the same samples showed very little change 
in observed intensity. For samples swollen in various solvents the same 
procedure was used. 


Results and Discussion 


For thin polymer films, most of the scattering occurs due to orientation 
and density fluctuations. Stein and Wilson have considered this in detail 
for static scattering from unstrained samples, and using the assumption of 
random orientation fluctuations, arrive at the following expression for 
vertically polarized incident and scattered light 


ge = 
V,=K " fro =: r°dr 


uw 
in f f(r) [" v(r) + | = reart (1) 
45 a? hr 


where K is an angular independent proportionality constant, a is the mean 
polarizability per unit volume of the medium, 7? is the mean-squared differ- 
ence between the principal polarizabilities of the volume elements (assum- 
ing cylindrical symmetry), h = (47/X) sin (6/2), \ is the wavelength of the 
light in the medium, @ is the scattering angle, and y(r) is a correlation func- 
tion for fluctuations in the density. For the horizontal component of 
scattered light using vertically polarized incident, they obtain 


H, = Kile f f(r) E v(r) + | = rart (2) 


In order to justify the use of the random orientation fluctuation approxi- 
mation, we show the change in scattering at 6’ = 6.6° and 10° as a function 
of rotation of crossed polars for an unstrained sample in Figure 2. The 
angle ¢ characterized the angle between the polarizer and scattering plane. 
For a system where this approximation is not valid, a large variation be- 
tween 0° and 45° should be seen. The data indicates about a 10% change. 
Data on other types of polyethylene show changes of approximately 60%. 
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Fig. 2. Change in Rayleigh ratio (Hy) at @ = 10° as a function of angle £. 
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Fig. 3. Change in Rayleigh ratio (V) of @ = 10° as a function of log time at 2 = 0 and 
2 = 90°. 


Stein has proposed an extension of parts of this theory to explain time 
dependent phenomena.’ However, he limits his treatment to only the 
terms pertaining to orientational fluctuations. 

In Figure 3 we present data for the change in the Rayleigh ratio as a 
function of log time at 6 = 10° for V, scattering at 2 = 0 (stretching axis 
perpendicular to scattering plane) and 2 = 90° (parallel to scattering plane). 
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It is observed that over the time period from 1 sec. to 100 sec. considerable 
relaxation takes place. 

In igure 4, we present data for H, scattering at Q = 90° and 6 = 10°. 
The obvious difference between /y, and Jy, shows that Jy, is associated 
Further, since Jy, is changing 














with much shorter time effects than Jy,,. 
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Fig. 4. Change in Rayleigh ratio (H,) of @ = 10° asa function of log time at 2 = 90°. 





while 7, is constant, then the changes are primarily.a result of density 
rather than orientational fluctuations during this period of time. Using 
eq. (1), we then have 





sin hr 


ly, = mf y(r) dr (3) 
hr 





For an oriented sample with only density fluctuations, eq. (3) is rewritten as 
ly, = 7? S y(t) cos k(r-s)dr (4) 


Following the original theory of Stein and Norris,® this equation can be 






evaluated to give 







Dq2ar2 
Ty, (Q = 0, 6) = 7? (Q = 0) exp — or aul a) (5) 
2 J 2 a? ~ 4 2 T 
Ty, (Q = 90, 0) = 7° (Q = 90) exp — Pa ial — cos 0)? + 2 sin’ 





i) 


where y(r) has been replaced by an empirical Gaussian function proposed 
by Stein and Norris 






a? 


y(t) = exp — (= + 5 + )t (7) 





and Q is the angle characterizing the relative orientation between the 
stretching and scattering direction. The parameters a and ¢ specify charac- 
teristic correlation distances in the sample perpendicular and parallel to the 
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Fig. 5. Change in correlation distance a during relaxation at constant strain as a function 
of log time. 


stretching direction. In order to evaluate a, /y, data is compared at two 
és for 2 = O, and a is given by 
l Ty, (Q = 0, 6;) I/, 
n 
i= —] — ‘ 
rv/2| cos 0 — cos 8; 


To evaluate c, Jy, (Q = 90) at different 6’s are used simultaneously with the 
data from a, 
A? Ty, (Q = 90, 6) 


] 
"Ty, Q = 90, 6) 


/ 
9 


| l 
+ a*[(1 — cos 6)? — (1 — cos @)?] | — 


; ; (9 
sin?é, — sin?6, 
l‘inally, in order to evaluate 4? (Q = 0) and 7? (Q = 90) use is made of eqs. 
(5) and (6), 

R(6', 2 = 0) Ao! 
rFy, (Q = 0) 
R(O’, 2 = 90) rot 


7 (Q = 90) = (11 
A, )= Fy, @ = 90) , 


n? (Q = 0) = (10) 


quations (8-11) were used to evaluate the changes in a, c, and 7? as a 
function of time in our experiments. Using data taken from several runs, 
the change in these parameters is shown in Figures 5-7 as a function of log 
time. It is of interest to show all the results from the start of an experiment 
totheend. This is shown in Figure 8 on a linear time basis. These results 
were obtained at a strain rate of ~55%/sec. and at the strain .03. The 
following conclusions can be drawn. During the straining a decreases 
while 7? (Q = 0) increases and simultaneously c increases while 7? (Q = 90) 
decreases, thus accounting for the behavior of the Rayleigh ratio. How- 
ever, during the subsequent relaxation, a@ continues to decrease while 7? 
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Fig. 7. Change in n? at 2 = 0 and 2 = 90 during relaxation at constant strain as a func- 


tion of log time. 


(2 = 0) decreases and similarly ¢ increases while 7? increases. It would 
thus appear that the major cause for the change during a relaxation experi- 
ment is connected with the n? term in eqs. (8) and (9). Also the final state 
of the sample after relaxation gives approximately the same results as 


noted by Stein and Norris, i.e., a general decrease in the scattering level for 
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Fig. 8. Change in Rayleigh ratio (Vy) at 2 = 0 and 2 = 90 from beginning to end of 
experiment on linear time scale. 
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Fig. 9. Effect of swelling with various solvents on Rayleigh ratios during relaxation 
experiments. 


Iy, at 2 = 0,2 = 90.6 The magnitudes of the changes are also in accord 
with the findings of Stein and Norris.6 Several explanations can be put 
forth to explain the change in 72, which we have found. 

It may be due to: (1) voids in the sample; (2) fluctuations in refractive 
index between crystalline and amorphous material; (3) melting of crystal- 
line material during stretching with recrystallization and orientation taking 
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Fig. 10., Stress relaxation curves for swollen and unswollen samples. 


place during relaxation; (4) rearrangement of crystals within the spherulite ; 
or (4) differences in orientation of aggregates of crystals within the spheru- 
lite. 

In order to eliminate some of these possibilities, samples were swollen in 
various solvents varying in refractive index from 1.46 to 1.56. Since it 
has been previously shown that the solvent does not enter the crystalline 
areas, the swelling must take place in the amorphous material.’ A general 
decrease in scattering of 10-20% was observed. In Figure 9, we compare 
relaxation data for samples at the same scattering angle, but swollen with 
different solvents (as indicated). It must be pointed out that depending 
upon the angular region which is observed, the effect of swelling varies. It 
is found that as noted by Stein et al., that the lower angles are affected more 
by swelling.’ Thus, in the angular region discussed here, the correlation 
size parameter appears to decrease on swelling. During the subsequent 
straining, a increases slightly and then decreases. This is due in part to the 
fact that the sample deswells. 

A separate series of experiments were done to measure the stress relaxa- 
tion in unswollen and swollen samples. The results of these experiments 
are shown in Figure 10. From the data, it is obvious that although the 
same strain rates and final strain were approximately equivalent, the stress 
isnot the same. This effect is caused by nonequivalent swelling ratios and 
different degrees of plasticization of the swelling liquid. However, in 
general, it appears that swelling does cause faster relaxation. Another 
effect noted was an increase of stress around 100 sec. which appeared to be 
caused by deswelling of the samples, since, as the solvent is removed, the 
sample contracts. Thus, if the rate of contraction is faster than the stress 
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relaxation, then the stress will increase. Finally, a comparison of 7? and 
stress as a function of time for unswollen samples shows them to behave in a 
similar manner, and in fact, the ratio of these quantities is independent of 
time, similar to the stress optical coefficient in rubbers. 


Conclusions 


It was suggested previously that when a film of polyethylene was de- 
formed uniaxially, that the crystallites would both rotate and translate.” 
Evidence has been previously offered which indicated ~ orientational 
relaxation times of crystallites are approximately .032 sec.* Our current 
data suggests that after the orientation of the crystals has stopped, that a 
subsequent translational motion of both crystalline and amorphous material 
takes place from approximately 1 sec. to more than 100 sec. This motion 
is in some way related on a molecular scale to the simultaneous relaxation of 
stress which occurs within the material as indicated by the constancy of the 
Rayleigh scattering optical coefficient in this region of time. Finally, the 
data shows that long-time relaxation in this material occurs over very large 
distances and is related to the mobility of the crystalline and amorphous 


material. 
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Résumé 


On a déterminé la relaxation de la lumiére diffusée dans des films minces de poly- 
éihyléne sous tension. Des résultats pour des films gonflés et nongonflés ont été pré- 
sentés. Des études simultanées de relaxation de tension montrent que le mécanisme par 
leque! a lieu la relaxation est associé & des changements dans la fluctuation moyenne de 
lindice de réfraction du produit. Le rapport entre le changement de la fluctuation de 
lindice de réfraction et la tension semble étre indépendant du temps. 


Zusammenfassung 


Die Lichstreuungsrelaxation in diinnen gedelinten Polyiithylenfilmen wurde bestimmt. 
Ergebnisse an gequollenen und ungequollenen Filmen werden vorgelegt. Eine gleich- 
zeitig durchgefiihrte Untersuchung der Spannungsrelaxation zeigt, dass der Relaxa- 
tionsmechanismus mit einer Anderung der mittleren Brechungsindexschwankung des 
Materials verknupft ist. Das Verhiiltnis der Anderung der Brechungsindexschwankung 


zur Spannung scheint von der Zeit unabhiingig zu sein. 
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Morphological and Rheological Studies of 
Polyethylene by Light Scattering 


ROBERT 8. MOORE and SHIRO MATSUOKA, Bell Telephone 
Laboratories, Incorporated, Murray Hill, New Jersey 


Synopsis 


The technique of photographic light scattering was employed to study films of medium 
and high molecular weight linear polyethylene. Results on the high molecular weight 
polymer show that secondary crystallization is severely restricted, that oxidation induces 
further crystallization, and that annealing may lead to tilted chain lamellae. Results of 
experiments conducted on films of the medium molecular weight polymer while subjected 
to tensile stress indicate that spherulite deformation occurs at very small strains in agree- 
ment with proposed theories. 


Introduction 


The technique of photographic light scattering'~* has been shown 
to be an excellent method for observing differences in the. structure of 
polyethylene spherulites. Modification of this technique has permitted 
us to obtain light scattering patterns from films subjected to tensile stress 
as well as scattering patterns from undeformed films. During these experi- 
ments it was found that scattering was not readily observed from films of 
very high molecular weight linear polyethylene whereas the usual scattering 
patterns were observed from films of lower molecular weight. Because 
this meant that crystallization was restricted in the high molecular weight 
polyethylene further studies were carried out by the technique of photo- 
graphic light scattering. This paper presents the results of these studies as 
well as an interpretation of the light scattering measurements and the 
mechanical behavior of a sample of lower molecular weight polyethylene 
subjected to tensile stress. 


Experimental 


Some of the characteristics of the materials used in these experiments are 
listed in Table I. As can be seen from the table both the medium molecular 
weight linear polyethylene, MMWL (Philips Petroleum Co.) and the high 
molecular weight linear polyethylene, H\IWL (Allied Chemical Co.) were 
extremely linear with fewer than three methyl groups/1,000 carbon atoms. 
The values of the density and x-ray degree of crystallinity of the HMWL 
polyethylene correspond to the values observed in a branched polyethylene 

163 






















RK. S. MOORE AND 8S. MATSUOKA 








TABLE I 
X-ray Weight- 
q Melting average 
Density," crystal- point, molecular 
Polymer g./ce. linity °C. CH;/100 C» weight 
DYNH (Union Carbide 
Co.) 0.932 65 112.0 2.4 1x 10° 
Marlex 6000-60 
(Philips Pet. Co.) 0.960 sO 137 .5 <0.3 1 X 10 
ACX Lot 13 
(Allied Chem. Co.) 0.933 67 138.5 <0.3 > 108 


* As determined at room temperature. 
» Methyl] groups/100 carbon atoms.® 





of medium molecular weight with about 24 methyl groups/1,000 carbon 





atoms. 

In these experiments a pulsed optical maser with a 90° oriented ruby 
crystal was used as the source of coherent, collimated, plane-polarized light 
(6943 A.) (see Fig. 1). The emitted energy ranged from 0.2 to 0.5 joules 
and the pulse duration was about 0.5 msec. The beam from the 5 by 0.6 
cm. ruby rod was passed through the circular aperture to remove uncolli- 
mated fluorescent radiation and onto the vertically mounted sample. 
The scattered radiation was transmitted to the photographic film (ASA 
speed 400) via the analyzer whose plane of polarization was kept either 
parallel or perpendicular to the horizontal plane of polarization of the 
incident beam. The sample to film distance ranged from 12-39 cm. for 
these experiments. Those samples subjected to tensile stress were held in 
the tensile creep apparatus shown schematically in Figure 1. This 
consisted of a platform with two vertical rods to which were mounted 
two horizontal cross bars, the upper of which was clamped. The lower 
bar was mounted on two linear movement bearings by means of which 
the bar moved with slight friction on the vertical rods. As shown in 
ligure 1 a weight could be attached to vary the sample stress. The 
sample, of nominal dimensions 7.6 by 0.76 by 0.005-0.02 em., was clamped 
at the upper and lower ends between the bars. The lower bar was held on a 
















Fig. 1. Light scattering apparatus showing the light source, the defining aperture, the 
tensile creep apparatus, and the camera bellows. 
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platform which could be rapidly lowered to start the experiment. Not 
shown in Figure 1 is a direct current differential transducer by means of 
which sample elongation was recorded on a strip chart recorder. Vor 
static experiments the creep apparatus was replaced by a holding frame 
to which were attached samples meited between cover slides to reduce sur- 
face scattering. Tor both types of experiments samples of the MMWL 


polyethylene were prepared from films molded at 160°C., while samples of 
the HMWL polyethylene were prepared from films molded at 190°C. 
between sheets of polyethylene terephthalate to facilitate release. Samples 
for the static experiments were melted between cover slides at 260°C. for 
20 see. and quenched to room temperature on a metal block. These 
samples were subsequently annealed under various conditions as described 


in the next section. 


Results 


Shown in Figure 2 are the scattering patterns from a quenched MMWL 
polyethylene sample and the scattering patterns from a quenched H\IWL 
polyethylene sample after it had been annealed at 125°C. for 90 min. For 
the MMWL sample the pattern observed with parallel polarizers is elon- 


(d) 


Parallel Polarizers Crossed Polarizers 


Fig. 2. Scattering patterns from (a,b) quenched MMWL polyethylene and (c,@) annealed 
HMWL polyethylene. 
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Fig. 3. Dilatometric measurements of the time dependence of the specific volume of 
two polyethylenes at 125°C. The arrows indicate the calculated time for spherulite 


impingement. 


gated in the direction of polarization, while the pattern observed with 
crossed polarizers resembles a clover leaf with leaves at 45° to the hori- 
zontal and vertical axes, as observed previously.2, The distence to one of 
the maxima is related to the average size of the spherulites,? the radius of 
the average spherulite being about 5.0 u as determined from this scattering 
pattern. 

During the primary stage of crystallization of spherulitic polymers the 
lamellar fibrils branch out from the center until they impinge on fibrils from 
other growing spherulites.’?. After this the secondary stage of crystallization 
predominates and the degree of crystallinity continues to increase.’ The 
latter phenomenon has been attributed to crystallization of the material 
between fibrils’ and to thickening and perfection of the fibrils themselves. °:!° 
In polyethylene the c or chain axis and the a axis are tangential while the b 
axis lies in the direction of the radial fibrils. 

It has been shown? that the development of the clover leaf pattern 
(see Fig. 2) occurs at about the same time that the pattern under parallel 
polarizer conditions changes from a circular pattern to one elongated in the 
direction of polarization. These changes occur during the secondary 
crystallization process when most of the chains are incorporated into fibrils 
with the chains lying tangential. The scattering observed under crossed 
polarizer conditions depends entirely on the existence of a difference be- 
tween the tangential and radial refractive indices within the spherulites, 
and no observable pattern exists unless sufficient anisotropy is present. 





POLYETHYLENE 167 


Experimentally, a large increase in the intensity of scattering was observed 
subsequent to spherulite impingement* and was interpreted as arising from 
secondary crystallization of the material in the interfibrillar regions.? 
Thickening and perfection of the fibrils would also lead to enhanced scatter- 
ing. 

The lower pictures shown in ligure 2 are for the HMWL polyethylene 
under similar conditions. It can be seen that essentially no secondary 
crystallization has occurred as evidenced by the absence of a clover leaf 
scattering pattern even though the sample had been annealed at 125°C. for 
90 min. A low intensity scattering pattern was obtained from this sample 
with the effective exposure about an order of magnitude higher indicating 
that the secondary crystallization was severely restricted. 
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Fig. 4. Comparison of the tensile creep compliance at various stresses with data de- 
rived from dynamic mechanical measurements. Stress in units of 10° dynes/cm.’: 
(@), 0.40; (O), 0.94; (®), 1.2. ( ), derived data; (---), derived data displaced ver- 
tically; (—), times at which scattering patterns in Figure 8 were taken. 


In Figure 3 are shown plots of the specific volume against the logarithm 
of time obtained from dilatometric measurements at 125°C. for two poly- 
mers comparable to the HMWL and MMWL polyethylenes. The arrows 
indicate the time at which spherulite impingement was calculated to occur. 
The time of impingement was determined from knowledge of the radial 
growth rate and the average spherulite size obtained by microscopy.’ 
These dilatometric measurements indicate a marked difference in the 
behavior of HMWL and MMWL polyethylene during crystallization. 
The difference is more pronounced after impingement has occurred when 
the specific volume of the MMWL polyethylene changes rapidly, but the 
specific volume of the HMWL polyethylene changes much more slowly 
and to a much lesser extent. Other evidence for restricted crystallization 
has been obtained from studies of the oxidative annealing,® the permeabil- 
ity,® and the dynamic mechanical properties of HMWL polyethylene!! 
and will be referred to in a latter section of this paper dealing with static 
light scattering experiments. 
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Since interpretation of scattering from samples of HMWL polyethylene 
subjected to tensile stress was complicated by the lack of secondary crystal- 
lization, MMWL polyethylene was chosen from simultaneous creep and 
light scattering studies. The results of the creep measurements were con- 
verted to logarithmic plots of the tensile creep compliance D against time 
in seconds t as shown in Figure 4 for the MMWL polyethylene. The lower 
solid line in Figure 4 was calculated with the method of Hopkins and 
Hamming"? from dynamic mechanical data obtained on a modified Maxwell 
Dynamic Tester.°:"' These results were shifted vertically to more easily 
compare the shape of the curve with the results obtained from the creep 
measurements. The creep compliance is observed to increase gradually 
from about 10>" to 10“? em.?/dyne over about four decades at room tem- 
perature. The samples subjected to creep were more compl ant than those 
subjected to dynamic stressing, due to the higher strains involved in the 
creep experiments. 

As can be seen in Figure 4 the compliances are not coincident at long 
times and high strains which indicates a departure from linear viscoelastic 
behavior. Though the nature of the two experiments is somewhat. dif- 
ferent because simple tension occurs in the creep experiments and simul- 
taneous tension and compression occur in the flexing of the rotating rod 
in the Maxwell Dynamic Tester,'* the effect of this difference on the com- 
pliance and modulus is expected to be quite small. 

The nonlinear behavior is more likely to arise from slippage and micro- 
yielding within spherulites. (This is not to be confused with the macro- 
scopic phenomenon commonly called necking.) Micro-yielding would arise 
from the difference in the response of the fibrils perpendicular and parallel 
to the direction of stress in a manner similar to that observed by Horio’ in 
‘ast films of polyethylene. The behavior of the load versus elongation 
curves for this cast polyethylene indicated that a large increase in stress 
occurred for a sample stretched in the a direction, the direction of casting, 
whereas the stress remained almost constant as the sample yielded when 
strained at right angles to the direction of casting, the 6 direction.'4 The 
stress at small strain required to produce this behavior" was about 10° 
dynes/em.? However, the actual stress required to cause yielding in our 
experiments would be below this value because there are fewer fibrils per 
unit volume and because the required stress level decreases with increasing 
time and strain. It can be seen from Figure 4 that under a stress of 1.2 X 
10° dynes/em.? the compliance reaches a value of about 10~*° em.?/dyne 
in 102° sec.; whereas at a stress level of 4 X 10’ dynes/cm.’, a value well 
below that believed necessary to cause micro-yielding at low strains, a 
compliance of 10~** em.? dyne is reached only after 10?7 sec. Thus, the 
results of our creep experiments, which encompass the required stress 
level, probably reflect micro-yielding. While micro-yielding would be 
likely to cause spherulite deformation in a sample under the required 
stress such deformation could occur at small strains for other reasons as 


discussed below. 
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Fig. 5. Temperature dependence of the tensile loss modulus at one cycle per second for 
four polyethylenes. 


Results of dynamic mechanical experiments on these polymers are pre- 
sented in Figure 5 in which the loss modulus at one cycle per second as 
determined on a modified Maxwell Dynamic Tester is plotted as a function 
of temperature. Two maxima are observed: the a transition which occurs 
from 30 to 80°C., and the 8 transition which occurs at about —30°C. 
The height of the 8 peak is observed to increase with an increase in the 
amorphous content of the polymer, indicating that the transition is sensitive 
to the amorphous regions. (The lower peak height of the HMWL poly- 
ethylene, due to restricted segmental motion is discussed in a later section.) 
Further evidence that the 8 transition is influenced by the amorphous 
regions is obtained from the temperature dependence of the a7 shift factors 
used to superpose dynamic measurements obtained at different tempera- 
tures.'' Determination of the coefficients of the WLI equation in the 
usual manner’ resulted in a value for c,° of 17.5 for a reference temperature 
of —27.5°C. This indicates that the mechanical behavior of polyethylene 
in the vicinity of the 8 peak is similar to the behavior of an amorphous 
polymer in the vicinity of its glass transition. The experiments of Nielsen" 
indicate that the crystalline regions also influence the 8 transition, and it is 
concluded that this transition involves motion of crystalline and amorphous 
regions coupled in series.'* Though the height of the 8 peak has sometimes 
been related to the degree of branching it can be seen from ligure 5 that 
annealing of the MMWL polyethylene causes a drastic reduction in the 
height of the 8 peak. 

The height of the a peak increases with increasing crystallinity and the 
maximum shifts to higher temperatures. This shift of the maximum with 
temperature is associated with an increase in the degree and perfection of 
the crystals, making motion more difficult until the temperature is in- 
creased. More important however is the fact that the a peak is due to a 
primary transition within the crystalline lamellae. The temperature 
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Fig. 6. Temperature dependence of the a dimension of the unit cell and of the de- 
gree of crystallinity for MMWL polyethylene. (Initially annealed at 125°C. for 5 


1 ys. ) 


of the maximum in the @ peak for each of the polyethylenes is comparable 
with the temperature at which the thermal expansion coefficient of the 
crystals undergoes a change as observed by the change in the a dimension 
with temperature as shown in Figure 6. The a dimension and the degree 
of crystallinity are plotted here as a function of temperature. The change 
in the slope of the plot of the a dimension versus temperature occurs 
at about 100°C., which is comparable with the maximum in the a peak, as 
has also been observed by Takayanagi.'® It can be seen that there is 
little change in the degree of crystallinity until well beyond the a peak of 
Figure 5, and the peak cannot be due to partial melting of the crystals. 
Rather, the peak must be due to a primary transition within the crystalline 
lamellae which has been transmitted throughout the sample. 

One would not expect to observe the effect of the intralamellar transition 
unless mechanical continuity existed from one spherulite to another. 
By this it is meant that the polymer is made up of spherulites which have 
impinged upon their neighbors with little or no amorphous material in the 
regions between the spherulites. This would result in the transmission of 
any change within the lamellae from one spherulite to another, and so on 
throughout the sample. In contrast to this, an amorphous polymer loaded 
with rigid spheres would not be expected to reveal a transition within the 
spheres since the effect of the transition would be taken up by the amor- 
phous material because the spheres would lack mechanical continuity. 

The fact that two peaks are observed over a comparatively narrow tem- 
perature range —30 to +20°C. for some branched polyethylenes) indicates 
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Fig. 7. Schematic diagram of a cross section of a spherulitic polymer under crossed 
polarizers subjected to vertical stress. The a and @ lines represent parallel viscoelastic 
response mechanisms discussed in the text. 


that there are two viscoelastic response mechanisms acting in parallel. 
A fringed micelle or series model would give one broad peak (unless the 
responses occurred at widely different temperatures) since the strains would 
be additive from the crystalline and amorphous portions. Moreover, a 
series model would show the compliance of the amorphous regions by a low- 
frequency storage modulus characteristic of a rubbery material, about 10° 
dynes/em.”. Because a low frequency storage modulus of 10°-10'° dynes 

cm.” 1s observed even for 60% crystalline polyethylene'! the response 
mechanisms must be coupled in parallel in order for the high modulus of 
the crystalline fibrils to be observed. 

Any morphological model would have to be consistent with all of the 
above observations. The model proposed is shown in Figure 7 and is 
similar to one proposed by Takayanagi.”” The diagram in Figure 7 
represents a schematic cross section of a spherulitic polymer under crossed 
polarizers with the sample stressed in the vertical direction. The two 
viscoelastic response mechanisms are represented by the vertical solid lines 
a, and the vertical dashed lines 8. The a lines represent the mechanical 
continuity from one spherulite to another through the vertical crystalline 
fibrils. ‘These crystalline fibrils would transmit the effect of the primary 
intralamellar transition throughout the system and would give rise to a 
high storage modulus. The amorphous regions and the crystalline regions 
at right angles to the vertical fibrils, designated by the dashed lines in 
Figure 7, would reflect the large compliance of the amorphous regions 
and are to be associated with the 8 mechanism which must include motions 
of both regions in series. Because of the low compliance of the fibrils and 
because of the difference in the viscoelastic responses parallel and _ per- 
pendicular to the direction of stress, spherulite deformation would have to 
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occur at very small strains for the model to be correct.’ Experimental 
determination of spherulite deformation is discussed in the next section. 
Photographs of the light scattering patterns taken at the times indicated 
in Figure 4 are shown in Figure 8 for the sample of MMWL polyethylene. 
Since the direction of sample elongation was vertical one would expect a 
horizontally elongated scattering pattern. In fact, Stein and co-workers 
observed an abrupt change from an elongated scattering pattern to two 
horizontal streaks at about 200% elongation which indicated that complete 
disruption of the spherulitie structure had occurred.? Two effects can be 
seen in Figure 8 as the amount of strain is increased. ‘The patterns 
gradually become more diffuse and the clover leaves gradually compress 
toward the horizontal axis, until at about 6% strain the effect is quite 
marked. (The effects are more easily seen if the pattern at zero strain is 
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Fig. 8. Scattering patterns for MMWL polyethylene subjected to a tensile stress of 1.2 < 
10° dynes/cm.*. The per cent strain is indicated. 


compared with each of the other patterns.) The first indication of de- 
formed spherulites is detected at about 3% strain, though some deforma- 
tion probably occurs before this. Thus, as would be expected from the 
theory, spherulite deformation occurs at very small strains. 

The picture on the lower right in Figure 8 was taken shortly after the 
sample was released from stress. As can be seen the clover leaves are not 
as compressed and the picture is less diffuse. Pictures taken at later times 
show essentially complete recovery of the spherulites despite the fact that 
the strain has not returned to zero. This might arise from a process 
in which strained crystals unfold and reform within the spherulites during 
sample elongation, the rate of recrystallization not being fast enough to 
keep up with the increasing strain. Upon release of stress the reforming 
process would dominate and recovery of the light scattering pattern 
would occur. Residual strain would then be taken up in a nonrandom dis- 
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Fig. 9. Scattering patterns (a) from annealed MMWL polyethylene and (6) from 
HMWL polyethylene after further annealing. 


tribution of amorphous material between fibrils. ‘The amorphous material 
would be rather insensitive to light scattering, however. 

In order to further investigate the restricted crystallization in the HMWL 
polymer the sample of Figure 2 was further annealed at 130°C. for about 
three hours. ‘The scattering patterns from this sample are compared in 
Figure 9 with the patterns from the sample of M\MWL polyethylene after 
annealing at 125°C. for 90 min. in the same air oven. The distribution of 
scattering in the upper pictures of igure 9 for MMWL polyethylene is 
narrower than the distribution for the quenched MMWL polyethylene, as 
shown in Figure 2, indicating that the spherulites are more uniform in size 
and that they are slightly larger. ‘The lower pictures in Figure 9, the one 
on the right being at a slightly higher output of the maser, show that a 
clover leaf pattern has developed in this sample of H\IWL polyethylene 
upon annealing, though the weak intensity of the pattern indicates that the 
extent of secondary crystallization is still small. 

In order to determine the effect of oxidation, a sample of HMWL poly- 
ethylene was melted between cover slips, quenched, and subjected to an 
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Fig. 10. Scattering patterns from (a) a quenched and (6) oxidized sample of HMWL 
polyethylene. 


oxygen atmosphere for 17 hr. at 113°C. At this temperature little re- 
crystallization would be expected to occur. Any effect would then be due 
mainly to crystallization of previously uncrystallized material. Scattering 
patterns from this sample are shown in Figure 10. The top picture of the 
quenched sample shows no scattering characteristic of secondary crystalliza- 
tion. The lower photographs were taken after oxidation, the picture on the 
left being of the scattering from the center of the sample while the picture 
on the right is of scattering from the periphery of the sample. The edge 
of the sample became whitened indicating that oxidation had occurred. 
l'rom the striking difference in the two patterns it is clear that oxidation 
greatly enhanced the secondary crystallization, probably through scission 
of polymer chains imbedded in more than one lamella"! 

The sensitivity of HM WL polyethylene to oxidative annealing has been 
observed by Winslow and co-workers.®:*'! A large increase in density 
occurred with a small amount of oxidation in the HMWL polyethylene 
whereas much more oxidation was required to achieve the same density 
and degree of crystallinity in lower molecular weight linear polyethylenes.?! 
Branching and crosslinking were found to impede the ordering process 
which proceeded rapidly, however, in the HMWL polyethylene with an 
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Fig. 11. Scattering pattern from a quenched (a) and vacuum-annealed (6, c) sam- 
ple of HMWL polyethylene; scattering pattern from a sample of polytetrafluoroethylene 


film (d). 


accompanying increase in the x-ray degree of crystallinity® from 67-95%. 
The similarity in total oxygen uptake was attributed to crystallization 
upon chain scission.’ Evidence for the orderly packing of chains upon 
oxidation was observed by Winslow and co-workers in electron photomicro- 
graphs of the replicated fracture surface of a sample of oxidized HMWL 
polyethylene. Long banded structures were visible with striations per- 
pendicular to the length of the bands.?!_ These banded structures are 
similar to those observed in low molecular weight polyethylene by Ander- 
son?” in which the chains are believed to be parallel to the striations. 
Oxidation might be expected to lead to similar bands in the case of HMWL 
polyethylene. 

To observe the effect of annealing without oxidation a sample of HMWL 
polyethylene was melted between cover slips and quenched to room 
temperature. It was then heated in a vacuum oven preflushed with nitro- 
gen at 130°C. for 24 hr. and allowed to come to room temperature under 
vacuum. Shown in Figure 11 are the seattering patterns obtained with 
crossed polarizers for this sample. No scattering characteristic of second- 
ary crystallization can be observed for the quenched sample. Scattering 
under the same conditions for the annealed sample are shown in the picture 
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at the upper right and at a higher intensity of the maser beam but at a 
further distance from the sample to film in the picture on the lower left. 
Both of these pictures show a cross rather than a clover leaf pattern. This 
is similar to the scattering pattern observed from films of polytetrafluoro- 
ethylene as shown on the lower right in Figure 11. Stein and co-workers 
have shown that this cross type scattering can occur if the scattering is from 
aggregates of crystals stacked at an angle to the direction of maximum 
polarizability, believed to be the case in polytetrafluoroethylene, or else 
it can occur if the scattering is from lamellae in which the chains are tilted 
at an oblique angle to the surface.2* The scattering in the annealed 
samples of H\{WL polyethylene could arise from tilting of the folded chains 
within the fibrils in a manner similar to that observed by Balta Calleja, 
Bassett, and Keller in solution-grown single crystals of polyethylenes 
of lower molecular weight.*4 This tilting might arise in the HMWL poly- 
ethylene in the thickening of lamellae having chains in common. 

The effect of restricted crystallization has also been observed in the 
mechanical properties of the HMWL polyethylene®:" in the behavior of the 
loss modulus £”’ at one cycle per second as a function of temperature (see 
Vig. 5). The height of the 8 peak increases with an increase in the amor- 
phous content for branched polyethylene and for annealed and unannealed 
MMWL polyethylene. However, the peak height of the HMWL poly- 
ethylene is less than that of the unannealed MMWL polyethylene and less 
than that of the branched polyethylene even though the latter material 
and the HMWL polyethylene have the same density and x-ray degree of 
erystallinity.'' Chains tied to more than one lamella would be quite 
restricted in their motion and would be expected to cause a reduction in 
the height of the 8 peak of the H\IWL polyethylene. 


Summary 


Through the use of photographic light scattering from thin films of high 
molecular weight linear polyethylene we have shown that secondary crystal- 
lization is severely restricted in this polymer, that oxidation induces further 
crystallization, and that annealing may lead to tilted chain lamellae. The 
restricted crystallization and the restricted motion observed in dynamic 
mechanical experiments probably have a common origin in the occurrence 
of lamellae with chains in common. Light scattering experiments con- 
ducted on films of medium molecular weight linear polyethylene subjected 
to tensile stress indicated that spherulite deformation occurred at very 
small strains, as predicted by Matsuoka’ for a system with two quite 


different viscoelastic responses acting in parallel. 
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Résumé 


La technique de la diffusion lumineuse a été employée pour |’étude de films de poly- 
éthyléne linéaire de poids moléculaire moyen et élevé. Les résultats de l'étude du poly- 
mére de poids moléculaire élevé montrent que la cristallisation secondaire est fortement 
empéchée, que l’oxydation induit une cristallisation ultérieurement et que le recuit peut 
conduire 4 des lamelles enchevétrées. Les résultats expérimentaux obtenus pour des 
films de polymére de poids moléculaire moyen indiquent qu’il y a dans le polymére, 
soumis 4 des tensions d’élongation, une déformation des sphérulites pour de tiés faibles 
tensions ce qui est en accord avec ies théories proposées. 


Zusammenfassung 


Das photographische Lichtstreuungsverfahren wurde zur Untersuchung von Filmen 
aus linearem Polyiithylen mittleren und hohen Molekulargewichts herangezogen. Die 
Ergebnisse an hochmolekularen Polymeren zeigen, dass die Sekundiirkristallisation stark 
eingeschriinkt, ist, dass Oxydation eine weitere Kristallisation induziert und dass Tem- 
perung zu gekippten Kettenlamellen fiihren kann. Versuchsergebnisse an Filmen aus 
Polymeren mittleren Molekulargewichts unter Zugspannung zeigen, dass, in Uberein- 
stimmung mit der Theorie, eine Sphirolithdeformation be sehr kleiner Verformung 
eintritt. 
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Studies of Rates of Spherulite Deformation 
by Low-Angle Light Scattering* 
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Synopsis 


The rates of spherulite deformation accompanying the elongation and contraction of 
polyethylene and the elongation of polybutene-1 samples have been observed using a low- 
angle light scattering technique. The light scattering pattern is obtained using a ruby 
pulsed laser which is flashed a predetermined time after the extension of the sample. 
It is found that spherulite deformation occurs in times less than 40 msec. which is less 
than the time required for the birefringence increase following deformation. Upon re- 
traction, the light scattering pattern becomes symmetrical before the sample length has 
returned to its initial value. With quenched polybutene-1, the rapid change in shape of 
the light scattering pattern is accompanied by a slower change in intensity occurring over 
times comparable with that for the birefringence change. This is associated with changes 
in the crystal structure, orientation, and density within the spherulite following deforma- 
tion. 


In previous publications,':* it was postulated that the initial orientation 
of crystals accompanying the deformation of a crystalline polymer is a 
consequence of the deformation of the spherulitic or crystalline super- 
structure within which they reside. A mechanism for such a deformation 
was proposed! which accounted for the change in birefringence and crystal 
orientation functions upon stretching the polymer. The crystal orienta- 
tion was found to be time dependent giving rise to time dependent bire- 
fringence changes that have been studied in experiments involving relaxa- 
tion at constant length or periodic straining. This paper represents an 
attempt to directly demonstrate this spherulite deformation and to establish 
whether the time dependence of crystal orientation is related to a time 
dependence of spherulite deformation. 

The direct microscopic observation of spherulite deformation and its 
time dependence is difficult. A photomicrograph illustrating the changes 
seen between crossed polars upon deforming a polyethylene spherulite is 
shown in Figure 1. The birefringence signs (determined from interference 

* Supported in part by contract with the Office of Naval Research and in part by 
Grants from the Air Force Office of Scientific Research and the Petroleum Research 
Fund. 

ft On leave from the General Electric Research Laboratory, Schenectady, New York. 
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Fig. 1. Spherulites in (a) unstretched and (b) 6% stretched polyethylene. Polar direc- 
tions are at +45° to the stretching direction which is indicated by an arrow. The sketch 
under the pictures indicates the sign of birefringence, where a (+ ) sign indicates that the 
greatest refractive index is for light polarized in the stretching direction. (=) indicates a 
more negative bir fringence than (— ). 


colors using a first-order red plate) are shown in the sketch below, where 
(+) signifies that the highest refractive index in that location is in the 
stretching direction and (—) the perpendicular direction. or polymers of 
practical interest, the spherulite size is usually too small to be easily 
resolvable by a light microscope, especially under conditions of samples 
being subjected to rapid strain. Consequently, low-angle photographic 
light scattering techniques were used. 

The relationship between the low-angle light-scattering pattern and 
spherulite size and anisotropy has been deseribed.*~* The deformation of 
the scattering pattern accompanying the stretching of the sample has been 
observed.* The initial change in the scattering pattern (at elongations up 
to 25-50% for low density polyethylene) may be associated with spherulite 
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Fig. 2. A photograph of the laser light-scattering deformation apparatus: (a), laser 
head; (6) adjustable pinhole for elimination of stray light: (c) sample held in deforma- 
tion apparatus; (d) photographic film holder and bellows; (¢) microswitch which closes 
upon completion of sample deformation; (f) traveling microscope for measuring sample 
strain, and (g), electronic timing circuit. } 


deformation while the changes at higher elongation may be associated 
with the reorganization of the crystal structure to fibrillar morphology.* 
The use of a pu'sed laser as a light source for these scattering measurements 
has been described.* This permits the obtaining of light scattering photo- 
graphs in less thar a millisecond. In this work, the laser technique was 
used for a sample subjected to rapid deformation and then relaxation at 
constant length. 

The sample in the form of a film about 1 in. X '/: in. X 5 mil. was held 
between clamps in a lever operated manual deformation apparatus (fab- 
ricated from a drill press attachment for a '/.-in. power home-workshop 
drill) which permitted stretching the sample at elongations up to 20% 
in 80 msee. (Fig. 2.) A microswitch was tripped by the deformation ap- 
paratus which triggered an electronic timing circuit which controlled the 
flashing of a small ruby crystal pulsed laser (Maser Optics, Boston, Mass., 
Model 600). The time interval between completion of stretching and 
flashing of the laser was controllable from 40 msec. upward. An elongation, 
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(a) ) 





(c) (d) 


Fig. 3. H, light scattering photographs for an undeformed polyethylene sample and for 
samples at varying time intervals after having rapidly been stretched 11-12%: (a) 
unstretched; (b) 0.042 sec.; (c) 0.082 sec.; (d)40 min. The polarization directions are 
indicated by arrows at +45° and the stretching direction is vertical. 


using a new sample, must be carried out for each time interval for the short 
time range experiments. <A polarized laser beam is obtained through use of 
a 90°-cut ruby crystal and a polaroid analyzer after sample is adjustable to 
either the Vy or Hy polarization conditions. * 

The sample studied was a Monsanto experimental medium density 
polyethylene similar to that used in our previous work.*4 Films were 
press molded at 140°C. at 630 psi pressure. The heat was turned off and 
the sample was allowed to cool in the press at its natural cooling rate. 

Hy scattering patterns for an unstretched sample and samples stretched 
11-12%, where the picture is taken at varying time intervals following de- 
formation, are shown in Figure 3. The polarization directions are +45° 
and the stretching direction is vertical. It is seen that upon stretching, 
there is a deformation of the general shape of the pattern from circular to 
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(b) 

Fig. 4. H, light scattering photographs for an undeformed polyethylene sample and for 
samples at varying time intervals after having rapidly been stretched at 15.2%. The 
polarization directions are vertical and horizontal and the stretching direction is vertical. 
(a) Unstretched; (b) 15.2% eat 0.051 sec; (c) 15.2% eat 20 min. 


elliptical accompanied by a decrease and then an increase in the intensity 
of the pattern. These changes are associated with (a) a deformation of the 
spherulite from spherical to ellipsoidal shape, and (b) a change in the anisot- 
ropy of the spherulite associated with reorientation of crystals within the 
spherulites. It is suggested that these two processes do not occur at 
equal rates. The deformation is indicated in the first picture obtained at 
0.042 sec. after deformation while further changes are evident in the later 
picture at 0.080 sec. No appreciable changes are seen between 0.080 sec. 
and 4 min. The resulting pattern is extended in a direction perpendicular 
to stretching direction as has been proposed.’ 

Dynamic birefringence measurements show that the birefringence 
relaxation time is of the order of one second for this material at room tem- 
perature. Thus the time for spherulite deformation is less than that 
required for the birefringence change following stretching. The deforma- 
tion of polyethylene appears to be a two stage process with the spherulite 
deformation occurring at very short times followed by an orientation 
process occurring at longer times, 
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Fig. 5. H, light scattering pattern for an unstretched polyethylene sample showing 
the relationship between A and 8; and the path of the radial densitometer scan (- - -) 
The intensity maxima lie along the lines at +45°. 


One obtains photographs in which the changes are less apparent upon 
quick visual examination but which are more subject to quantitative 
interpretation, by having the stretching and polarization directions co- 
incide (Fig. 4). The spherulites are again seen to deform very rapidly from 
spherical to ellipsoidal, with the ellipsoid extended in the stretching direc- 
tion—with the scattering pattern extended perpendicular to the stretching 
direction. Radial scans (lig. 5) using a microdensitometer reveal that the 
are length B (in the vertical direction) between the intensity peaks in- 
creases from 90° in the unstretched case to 104° at 0,051 sec., and 109° at 
20 min. The are length A (in the horizontal direction) decreases from 90° 
to 76° to 71° over the same time interval. The uncertainty in the measure- 
ments is about +3%. The reduction in intensity at short times followed 
by an increase in intensity at longer times is also apparent. 

Observations were made on the contraction of deformed polyethylene 
spherulites as the stretched polymer was allowed to retract freely. The 
sample was first rapidly stretched to 31.9% elongation and held at fixed 
length for one hour. Upon cutting the sample free and observing the 
changes in scattering pattern with time, the deformed pattern was found to 
approach that of unstretched polyethylene (Fig. 6). This result appears to 
be in agreement with that reported recently by Moore and Matsuoka.° 
The change in elongation of the sample on retraction was followed as a 
function of time, by observing the movement of bench marks on the sample 
using a traveling microscope (Fig. 7). The recovery of the deformed 
spherulite appears to be associated with the high rate of recovery portion of 
the curve. : 

It is noted, on comparing Figures 6 and 7, that the sample is still con- 
siderably strained when the spherulite pattern has returned to its unde- 
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(b) (ec) 


Fig. 6. H, light scattering photographs for a polyethylene sample stretched 31.9% 
and for the sample at varying time intervals after the start of free retraction: (a) before 
retraction; (b) 4sec.; (c)10 min. The polarization directions are vertical and horizontal 
and the stretching direction is vertical. 


formed shape. The mechanism by which this residual strain has occurred 
is not understood, and therefore further experiments are planned in an 
attempt to resolve this point. 

Isotactic polybutene-1 has two crystalline modifications, one of which is 
stable at room temperature. The other form is that occurring at the time of 
erystallization.""'! This latter form is very unstable at room temperature. 
It spontaneously changes over to the stable form. 

This transition rate is affected by such processes as aging, temperature, 
and mechanical stretching. '? 

Polybutene-1 samples were prepared by press molding between aluminum 
sheets at 156°C. at 400 psi pressure. The samples were quenched from the 
melt by waving them in air upon removal from the hot press. 

The time dependence of low angie light scattering from quenched poly- 
butene-1 was observed by the photographic technique. Hy light scattering 
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patterns for an unstretched quenched polybutene-1 sample (1 hr. after 
preparation) and for the sample held at 10% elongation after stretching 
are shown in Figure 8. After 0.041 sec. the spherulites are already de- 
formed from spherical to ellipsoidal, with the ellipsoid extended in the 
direction of stretch, as was the case with polyethylene. The later pictures 
(at 1 hr. and 10 hr.) show a marked increase in the intensity of scattering 
while the pattern substantially retains its deformed shape. The bire- 
fringence of this stretched sample was found to be positive and to increase 
by a factor of about 2 over the time period concurrent with the increase in 
scattered intensity. This is in contrast to the observation’? that the bire- 
fringence becomes negative upon stretching stable polybutene-1 which has 
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Fig. 7. The time-rate of retraction for the polyethylene sample for which light scattering 
photographs are shown in Figure 6. 


been aged in the undeformed state. The shift from negative to positive 
values of the birefringence is not a result of the change in birefringence 
accompanying the transition, but rather results from a change in the orienta- 
tion mechanism with aging. Evidence for this was obtained in the follow- 
ing manner. !? 

The mode of crystal orientation accompanying the stretching of a 
partly aged sample was studied by measuring the azimuthal dependence of 
x-ray diffraction. At high elongation, where the birefringence is positive, 
both the (110) and (300) reflections have maximum intensity at the equator, 
indicating that both the a and b axes orient perpendicular to the stretching 
direction. 

At low elongations, where the birefringence is negative, the (110) 
reflection is more intense at the meridian, while the (300) reflection is more 
intense at the equator. Thus, the change from negative to positive 
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(a) (b) 





(c) (d) 


Fig. 8. H, light scattering photographs for an unstretched freshly quenched polybu- 
tene-1 sample and for the sample at varying time intervals after being rapidly stretched 
10%. The polarization directions are vertical and horizontal and the stretching direc- 
tion is vertical. (a)0% esec.; (b) 10% eat 0.041 sec.; (c) 10% eat 1 hr.; (d) 10% e 
at 10 hr. 


birefringence on stretching an aged sample is a result of the shift from a b 
axis to c axis orientation mechanism. 

The crystalline contribution to the birefringence of a spherulite sample 
may be described by the equation! 


™ 3 cos?e — 1 . ( 
4A = A(@) N(0@) sine do (1) 
0=0 Z 


where A(@) is the spherulite birefringence per unit volume (difference 
between radial and tangential refractive index) contributed by those 
crystals in the volume element of width de lying on a spherulite radius 
making an angle @ to the stretching direction. In the unstretched state, 
A(@) is negative and independent of @. If the spherulite deformation is a 
completely affine process (as is approximated in an aged polybutene-1 
sample at low elongation), A(@) remains negative on stretching and the 








188 P. ERHARDT, K. SASAGURI, AND R. S. STEIN 





(a) (5) 


Fig. 9. H, photographic light scattering patterns for unstretched quenched polybutene-1: 
(a) fresh; (b) aged. The polarization directions are vertical and horizontal. 


negative sample birefringence is a consequence of the change in the distribu- 
tion function N(@) sin @ de (representing the number of spherulite radii in 
the volume element N(@) sin @ de) with elongation. 

For a freshly quenched polybutene-1 sample, A(@) changes with elonga- 
tion as a result of the change in orientation of the crystals within the 
spherulite from a state where their c-axis is tangential to one where their c- 
axis is radial accompanying the transformation of their structure from that 
of modification II to modification I. This transformation is induced by 
the internal strain of the spherulite and results in a change of A(e) from a 
negative to a positive value simultaneously with the change in N(@) pro- 
duced by stretching. This results in A becoming positive with elongation 
of such a sample. 

The transformation from modification II to I is time dependent and will 
not be complete at the time of completion of stretching. Hence, if the 
sample is held in this stretched state, A(@) continues to become more 
positive while N(¢) remains constant, resulting in A becoming more positive 
as observed. 

An increase in the seattered intensity was also observed in the spon- 
taneous transition of unstable to stable polybutene-1, in the unstretched 
ease (lig. 9). The transition in the unstretched state occurs at longer 
times after sample preparation, indicating that mechanical deformation 
enhances the transformation from unstable to stable polybutene-1, as has 
been observed from x-ray diffraction measurements. Densitometer scans 
of the photographic light scattering patterns of l’igure 9 show a general rise 
in the overall level of intensity rather than an increase of amplitude of 
the intensity peaks. This indicates, that, in the unstretched state, the 
increase in the intensity of scattering is due to density rather than orienta- 


tion fluctuations. 

This would appear to be in agreement with the reported observation" 
that the conversion of unstable to stable polybutene-1 occurs with a change 
of crystalline form of lower density to one of higher density. Assuming the 
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volume of the unstretched sample to remain constant, such a redistribution 
of density within the sample would give rise to an increase in density 
fluctuations as has been observed. 
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Résumé 


Les vitesses de déformation des sphérulites, qui accompagne I’élongation et la con- 
traction du polyéthyléne et l’élongation d’échantillons de polybuténe-1, ont été ob- 
servées au moyen d’une technique de diffusion lumineuse 4 petit angle. L’image de ia 
lumiere diffusée est obtenue au moyen de la lumitre provenart d’un laser au rubis 
envoyée pendant un temps trés court et prédeterminé apres |’extension de |’échantillon. 
On a trouvé que la déformation du sphérulite a lieu dans un temps inférieur 4 40 msec., 
ce qui est inférieure au temps requis pour que la biréfringence puisse faire augmenter la 
déformation subséquente. Lors de la rétraction, l’image de la lumiére diffusée devient 
symétrique avant que la longueur de !’échantillon ne soit revenue 4 sa valeur initiale. 
Avec le polybuténe-1 refroidi brutalement, le changement rapide de la forme de l'image 
de la lumitre diffusée s’accompagne d’un plus faible changement de |’intensité qui a lieu 
dans des temps comparables a celui necessaire pour le changement de la biréfringence. 
Cela s’accompagne de changements dans la structure du cristal, l’orientation et la 
densité du sphérulite aprés déformation. 


Zusammenfassung 


Die Geschwindigkeit der Sphirolithdeformation bei der Elongation und Kontraktion 
von Polyiithylen und der Elongation von Polybuten-1-Proben wurde mit einer Klein- 
winkel-Lichtstreuungsmethode gemessen. Das Lichtstreuungsdiagramm wird mit 
einem Rubinlaser erhalten, der zu einer vorbestimmten Zeit nach der Dehnung der 
Probe geziindet. wird. Die Spharolithdeformation erfolgt ia Zeiten kiirzer als 40 msec, 
was weniger als die fiir die Doppelbrechungszunahme nach der Deformation notwendige, 
Zeit ist. Beim Zusammenziehen wird das Lichstreuungsdiagvamm vor der Riickkehr der 
Probenlinge auf ihren urspriinglichen Wert symmetrisch. Bei abgeschrecktem Poly- 
buten-1 wird die rasche Gestaltsiinderung des Lichtstreuungsdiagramms von ciner 
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langsameren Intensitétsinderung begleitet, die sich iiber vergleichbare Zeiten wie die 
Doppelbrechungsinderung erstreckt. Sie steht in Beziehung zu Anderungeu der Kristall- 
struktur, Orientierung und Dichte innerhalb des Sphirolithes nach der Deformation. 
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